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NUMERICAL METHOD AND FORTRAN PROGRAM FOR THE SOLUTION OF AN 
AXISYMMETRIC ELECTROSTATIC COLLECTOR DESIGN PROBLEM 

by Oliver W. Reese 
Lewis Research Center 

SUMMARY 

This report describes the determination of the steady- state flow of electrons in an 
axisymmetric spherical collector under a variety of boundary conditions. The electron 
trajectory equations of motion are solved alternately with Poisson's equation for the 
potential field until convergence is achieved. 

The bounding surfaces are a sphere below, a spherical cone above, and, optionally, 
a thin axial spike suspended from the apex of the cone. Negative decelerating potentials 
are prescribed on the cone- spike surface and on segmented zones of the spherical sur- 
face. Polar angles of the cone may vary from 0° to 160°. 

The principal value of this report is that it pushes design research into the high- 
current- density range, where space charge effects are strong. Initial attempts met 
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with convergence difficulties for current densities above 50 A/cm . 

This report describes the method used to surmount these difficulties and handle 
values as high as 100 A/cm . This was achieved in two ways: first, by using variable 
mesh spacing to ensure accuracy in regions where the solution was changing rapidly; 
and second, by using a fast, accurate subprogram for obtaining the potential field. This 
subprogram uses a Wect, noniterative numerical solution to Poisson's equation, the 
details of which are given in NASA TN D-6438. 

The present report includes a description of the mathematical model, a discussion 
of numerical techniques, results from two typical runs, and the FORTRAN computer 
programs. 


INTRODUCTION 

Satellite communication systems have been proposed in which the satellite can re- 
ceive, amplify, and rebroadcast microwave signals. The high efficiency of the micro- 



wave amplifiers to be used is of prime importance to the lifetime and reliability of such 
a system. 

Microwave amplifiers, as a class, convert the energy of bunched electron beams 
into radiofrequency (rf) energy. Since this process is only partially efficient, the re- 
maining beam energy is dissipated as heat in the electron collector. However, collec- 
tors can be designed and maintained at proper potentials to slow the electrons. The 
electron energy is then spent doing work against the collector potential and is recover- 
able as useful electrical power. 

Reference 1 discusses how the results of the program described in this report are 
used in collector design. Reference 2 describes a computer program which solves 
Poisson's equation for a given space charge and which is required as a subroutine in the 
iterative solution of the current problem. 

It is important to realize that there are two characteristics of electron beam trajec- 
tory patterns that facilitate efficient collector design: 

(1) Good, uniform, vertical and horizontal spreading (i.e. , a sufficient rise away 

from the entry hole, and a pushing away from the collector axis) 

(2) Good, uniform separation permitting insertion of electrodes to collect electrons 

near the peak of their trajectories, where their energy is close to minimiun 
Both these characteristics imply no backstreaming, crossing, or impingement of beams 
on the underside of the electrodes. 

The original computer program dealt with negative decelerating potentials pre- 
scribed only on segmented zones of the sphere. But this approach was ineffective be- 
cause of the great distance between these surface potentials and the incoming stream of 
electrons. 

A second stage of the design was to prescribe a decelerating potential on a cone of 
varying polar angle. This slowed down the incoming stream considerably, and in most 
cases the beam spread was sufficient. 

In other instances where even more beam spread was desirable, a third and final 
design modification was made. Here a negative potential was prescribed on a thin spike 
suspended from the cone. This had the effect of spreading the flow even more as it 
entered, and allowing the other decelerating potentials to come into play more effec- 
tively. 

The final FORTRAN program reported herein allows for the presence or absence of 
an axial spike of varying length, as well as for the choice of decelerating surface poten- 
tials on zonal segments of the sphere and on the cone- spike surface. The cone angle 
varies from 0° to 160°. Initial conditions also include variable velocity and injection 
positions of entering electron groups. 

It was assumed that, prior to entry into the collector, the electron beam had been 
sorted into graded velocity groups, with the outermost being the slowest. This is phys- 
ically more realistic than assuming a constant velocity only. The beam was then 


2 



injected through a small entry hole in the collector where th'e negative polar axis pierces 
the spherical surface. The potential distribution on the surface of the sphere was as- 
sumed to be a function of the polar angle only. This report describes the numerical 
solution to a problem used in the design of such a collector under varying conditions of 
surface potential, current density, and collector geometry. 

A typical configuration is shown in figure 1. Notice that the electron groups reach 
minimum velocity at their peaks and start falling back. If electrodes of the proper 
potential are inserted to collect the electrons at these minimum velocity points, the 
efficiency of the collector will then be maximized. 



Electron 

trajectories 


Figure 1. - Typical collector geometry (spherical cone). 


The numerical technique and computer program described in reference 2 were used 
to determine the potential field inside the collector as required. The results were then 
used as part of the design criteria of the collector as described in reference 1. 

In situations where current density is low, less than 10 A/cm , space charge effects 
were found to be negligible. But with higher current densities and associated higher 
space charge, certain problems arise near the entry hole that require special techniques 
to ensure numerical stability. This report describes these techniques. Current densi- 
ties as high as 100 A/cm were dealt with successfully. Initially, calculations were suc- 
cessful only with current densities half as large, before numerical difficulties precluded 
further computation. 

A complete description of the mathematical model is given, along with a discussion 
of the assumptions made. The technique used in determining the space charge distribu- 
tion is then described. Appendixes A to G contain the results of two typical computer 
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runs, a complete listing of the FORTRAN computer program, flow charts, magnetic 
field calculations, the derivation of the equations of motion, and the mathematical 
symbols. 


STATE/VIENT OF THE PROBLEM 

The problem is to determine the steady- state potential field and the associated 
steady- state electron beam trajectory patterns within a spherical collector (1) for any 
specified collector geometry (cone angle; collector radius; entry hole radius; and length 
of axial spike, if any), (2) for negative potential distribution on the surfaces of the col- 
lector, and (3) for magnetic field distribution and electron beam initial entry conditions 
(current, cinrrent density, perveance, injection angles, positioning, and velocity). 

This involves the determination at every mesh point in the sphere of the potential V 
that satisfies Poisson's equation 

V^v = f (V) 

where 

r,2y _ 3^V , 2 3V ^ 1 ^ cot 9 ev , 1 

0p2 p dp ^2 gg2 -2 36 -2 g.j^2^ g^2 

The bars refer to dimensioned values. All mathematical symbols are defined in appen- 
dix G. For convenience, these expressions and all others following are normalized by 
the substitution 

p (Ic) 

R 


(la) 


(lb) 


We may now write 


f(V) (Id) 

"0 

The right-hand term in equation (Id) is related to space charge resulting from the 
beam trajectories. These trajectories are obtained from the solution to the differential 
equation of motion (eq. (3)). Note that the determination of V is an iterative process. 
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in that the right-hand side (or source term), which depends on the space charge, 
changes when the internal potential V changes. 

The source term consists of the dielectric constant £q and Pg(V) where 

PAY) =-jI| (2) 

H 

The velocity term u is easily calculated since it is simply a function of certain initial 
conditions and the potentials V and . 

The current density term J is somewhat more complicated. In general, it is cal- 
culated from the current flow through annuli formed by adjacent trajectories of the vari- 
ous electron energy groups. This flow is shown later in figure 4. These trajectories 
are calculated from the differential equation of motion 

A ^ pM ^ p + P sin^e (p^ ^ _7L A 9V _ 0 ^ _ p' sin 6 cos 6 

gg2 ^2 ^2j^\R9p 36/ ^2 

(3) 

P p2^2p2\30 dp/ 


where 

p normalized radius vector 

9 polar angle 

<p azimuthal angle 

R radius of sphere 

Vjjj magnetic field 

Details of this iterative process are described in the next section. 


METHOD OF SOLUTION 

First, the steady- state flow of electrons with a given velocity distribution was de- 
termined within the collector. This was accomplished by taking the following four steps: 
(1) With the given collector geometry and a specified set of initial conditions, a 
first approximation to the potential field inside the collector was found by assuming zero 
space charge and solving Laplace's equation 
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v^v = 0 


(4) 


or 


d^v ^ 2 dV ^ 1 d^V ^ cot g 5V _ Q 
gp2 p 3p p2 gg2 p2 36 

(2) An arbitrary number of equally spaced electron trajectories were then calculated 
from the equation of motion (eq. (3)). This equation is derived in appendix F. In gen- 
eral, nine electron groups were solved for in this problem. 

(3) From these trajectory patterns, the space charge was then calculated at each 
mesh point, and the results were used as source terms in Poisson's equation which was 
solved next. 


V^V=— (6) 

"0 

Details of how the space charge was calculated are given in a later section. 

(4) Equation (3) was solved again to obtain new trajectories based on the new poten- 
tial field obtained from equation (6). This iterative process was then continued until 
convergence was achieved. The entire process is shown in flow chart form in figure 2. 

Prior to solving the problem, the following four sets of normalized initial conditions 
must be specified. (The quantities listed in the first three sets are shown in figure 3 
for the case of three electron classes. ) 

(1) Collector geometry: radius of collector R, initial polar angle 0 q, and 

radius of entry hole a 

(2) Collector boundary values of electric potential: on a spherical surface 

(Vi, V 2 , • • . , Vj) (In general, the potential on the spherical surface 
decreases linearly to zero (Vj = 0) as the polar angle increases from 0 q 
to some prescribed 6 ^. ); on a conical surface if 0 q 4- 0(Vj); and along the 
negative polar axial spike if 6q 4- O(Vj^) (optional) 

(3) Electron beam entry conditions: type of flow (Brillouin or confined; con- 

fined flow means there is no initial rotational velocity of the beam, or 
(p =0), values of beam energy classes v^, entry positions of electron 
classes in hole at injection 6^^ (6^ is the supplement of the polar angle 
0), and entry angles of^. In general, i = 1, 2, . . ., 9. 
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Given collector qeomet 
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(4) Electron beam parameters: magnetic flux density Bg, magnetic field 
potential electric potential Vg, current Ig, and current density 


Auxiliary Equations 

We now discuss the auxiliary equations which are required in the solution of equa- 
tion (3). 

The first two are the magnetic field terms expressed in cylindrical coordinates. 


av / 



av 


m 

sin 0 + 

m 

dp \ 



9z 




9V_ 

— H =pR 

m 

cos 0 - 

m 


ae 


0r 


3z 


cos 0 


sin 9 


(7a) 


(7b) 


or letting ~ 


9V, 


m 


dp 

av. 


m 


00 


= R(B sin 0 + B„ cos 0) 

iU 


pR(B cos 0 - B_ sin 0) 

X Z 


(8a) 


(8b) 


Prior to entering the collector, the electron beam passes through a magnetic field 
which is used to confine the beam. Part of this field leaks into the collector through the 
entry hole and influences the flow of the electrons. It must be considered in the equa- 
tions of motion. 

• 2 

Also used are the velocity terms, 0 and <p , and the magnetic flux ip . These are 
expressed as 


1 

1 + (p’/p)^ 

cp = E 

2np^R^ sin^0 


2r,eVg(Vi.V) 3 3 

<p sin 0 

p'r2 


(9) 


(10) 
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B^r dr 


( 11 ) 


^ = 



which are equations (B16), (B8), and (B9) of reference 1, respectively. 

Finally, the components B^, and B^ of the magnetic field (eqs. (D9) and 
(DIO) of ref. 1) are expressed as 


and 



B (r, z) =y X 9^ . "'^ B( p, . z) ^ /r\2n-2 

" Z^[(n-1):]2 3z2n-2 W 

n=l ^ 


where 


B(0,z) = Bq 



1 z/a 

’ 1 1 (zW) 


(14) 


The approximations of B_(r, z) and B_(r, z) that were used are given in appendix E. 


The Equation of Motion 

The second- order differential eqxiation of motion (eq. (3)) was transformed into two 
first-order equations and solved by using a fourth- order Runge-Kutta method. As is 
typical in integration techniques of this type, the choice of step size is based largely on 
experience. In this problem, the most critical region is near the entry hole. For a 
large collector radius (R = 16 cm) and a small entry hole (a = 0.05 cm), an initial step 
size A0 of 0. 01° was required. For a small collector radius (R = 8 cm) and a large 
entry hole (a = 0. 10 cm), an initial A0 of 0. 05° was sufficient. As the trajectory 
moved away from this critical region, the A0 step was increased to minimize computer 
execution time. 

Initially, the solution was obtained in the p-9 spherical plane. Then when any par- 
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ticular trajectory started to arc over (p' = 0), a change was made to the r-z cylindrical 
coordinate system to avoid any possible discontinuities in the integration process. 


Poisson's Equation 

For any space charge field, the set of right-hand sides for Poisson's equation writ- 
ten at all interior grid points may be computed, and the subroutine of reference 2 may 
be applied to solve Poisson's equation for the internal potential field. The method of 
calculating the space charge field from the trajectory field is now presented. 


Calculation of Space Charge 

The source term, or space charge term, in equation (6) is 

^0 


(15) 


where 

- 14 

€q dielectric constant, 8.86x10 F/cm (farads per centimeter) 

O 

Pg volume charge density. C/cm (coulombs per cubic centimeter) 

Furthermore, 

Pe =- (16) 

® u 


where 

2 

J current density, A/cm (amperes per square centimeter) 
u velocity, cm/sec 

Figure 4 shows where a tjrpical Pg is computed; namely, at a point midway between two 
trajectories. To ensure accuracy, space charge is calculated at selected points along 
center lines midway between adjacent trajectories where the p-mesh lines are crossed. 

Since u is known at various points along a trajectory, one may interpolate to com- 
pute u at any point, say at a position midway between two trajectories. Therefore, the 
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Polar 

axis 



Figure 4. - Cross-sectional view of entry hole in spherical collector, show- 
ing typical injection pattern of electron energy groups 2 to 9 at radii r 2 , 
r 3 , . . . rg. First subscript refers to electron trajectory number; second 
subscript j refers to p-mesh line crossed by trajectory. 


space charge at each mesh point may be determined once J is known, 
how this is accomplished. 

In general, 

J.I 

A 

where 

I current, A 

A any particular element of area through which current flows 
Referring again to equation (6), we may now write 


We now show 


(17a) 


(17b) 
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where 


Iq initial cvirrent at injection 
i i*^” electron energy group 

Referring again to figure 4, it is seen that the equally spaced electron energy groups 
form eight annuli at the entrance hole . The areas dA of these annuli may be expressed 

dAj = n^r| 

^2 ■ ^^3 ~ - dAj^ = 3 dA^ 

dAg = 7Tr\ - nr^ = iriZr^f - 7r{2r^)^ = 5 dAj 

dAg = ?rrg - ffr| = ^(Srg)^ - = 15 dAj 

Since I = JA, the current at injection through each annulus may be expressed 




^ 0 , 1 " *^0 ^1 

^ 0, 2 ^*^0 *^2 ^ 3 Jg dAj 

^ 0, 3 '^0 ^3 ^ ^"^0 ^1 

^ 0 , 8 ^ *^0 ^8 ^ ^^*^0 *^1 

^0,i ^ 



(19) 
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In the example of figure 4, the point shown represents a typical space coordinate 
midway between two adjacent trajectories, in this case, trajectories 5 and 6. The 
normal distance between the trajectories at this point is d. The area of the annulus 
is then 


A = 2jrrd = 2ffPc . sin 0- z d (20) 

o,j o>j 

The general expression for the current density, J = I/A (eq. (17a)), may now be 
written as 


(21- DJjdA, 
or 

/ (21 - l)r|j„ 

Notice that current is preserved in each annular re^on, n^unely 


l0_i = (21- l)JodAi=Ji_jAj_j 

Since A. . can be computed at (p, i, 0. J by 

\]'2'Vl.j“l»\i'‘i.i 

where d- ^ is the normal thickness of an annulus at this point. 
J 


(21a) 


(21b) 


( 22 ) 


(23) 


Use of Irregular Mesh Spacing 

The successful solution to this problem, especially those phases dealing with high 
current densities, was largely due to the use of irregular mesh spacing in both the p 
and 0 directions. From figure 5 it is clearly apparent that the region in which the 
solution varies most rapidly is immediately along and adjacent to the negative polar 
axis, where the entering electron beam is most dense and of greatest velocity. Conse- 
quently, a closely spaced mesh was used in this region. Conversely, a coarse spacing 
was used where the solution was varying very slowly and quite predictably. 
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Figure 5, - Cross section of right half of spherical cone-shaped 
collector, showing typical pattern of electrodes, electron tra- 
jectories, and irregular mesh spacing in both the p and 0 
directions. 


Convergence Acceleration 

Experience has shown that in solving this problem, the iterants are well behaved 
and ultimately converge to an acceptable approximation to the solution. But because of 
the steadily decaying oscillatory nature of the iterants, a considerable savings in com- 
puter execution time was achieved by halving the newly calculated changes in the vari- 
ables at each step. Normally, a maximum of five iterations would provide an acceptable 
solution. 

. FORTRAN PROGRAM 
General Description 

The program consists of the following 16 subroutines: 

U) MAIN is the executive routine for processing multiple cases and has primary 
control of logical flow throughout the complete program. 

(2) INPUT reads. and prints initial conditions, parameters, and other data required 
for internal use by the program (printout frequency, looping indices, convergence cri- 
teria, etc.). 
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(3) INTGRN controls integration looping of equations of motion, step size, change- 
over of coordinate systems, and the number of iterations required for convergence. 

(4) DE is the spherical coordinate version of a routine for computing trajectories by 
integrating the equations of motion. 

(5) DERIV calculates 8V/9p and 3V/36 for any specified interior point. 

(6) RK is the spherical coordinate version of a Runge-Kutta integration subroutine 
used to compute trajectories. 

(7) MESH generates graduated mesh point arrays in both the r and 9 directions 
for any specified collector geometry. 

(8) RZTRAJ is the cylindrical coordinate version of a routine for computing trajec- 
tories by integrating equations of motion. 

(9) RZRK is the cylindrical coordinate version of a Rimge-Kutta integration subrou- 
tine used to compute trajectories. 

(10) RZDE contains the equations of motion in the cylindrical coordinate system. 

It is similar to subroutine DE. 

(11) OUT controls the printing of results at specified intervals along the trajec- 
tories. 

(12) VFIELD provides for the calculation of the potential field for any particular 
set of boundary values and source terms. It uses the program described in reference 2. 

(13) EQLINE calculates the spherical coordinates of any specified number of equi- 
potential lines. It is used only after convergence is achieved in the main portion of the 
program. 

(14) RHSCAL generates space charge values along center lines of the various annuli 
formed by the trajectories. It interpolates for these values at each mesh point, as re- 
quired in the solution of Poisson’s equation. 

(15) PSICAL calculates the magnetic flux p of the magnetic field at each mesh 
point. 

(16) BRBZ calculates B^ and B^, the r and z components of the magnetic field, 
at each mesh point. 

All calculations were done on the IBM 7094 11/7044 Direct Couple System computer. 
The program will perform about 210 integration steps per minute. This figure is based 
on trajectory calculations for eight electron energy groups. 

About 15 000 storages are required for the program. This does not include that re- 
quired for the solution of Poisson's equation. An additional 10 000 storage locations 
are needed for a solution to Poisson's equation similar to that described in reference 2. 

The FORTRAN program is listed in appendix B, FORTRAN symbols are defined in 
appendix C, and flow charts are presented in appendix D. 
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Using the Program 


The program always starts by reading in a nine- card input data deck. The contents 
of each of these cards is now described. 

Card 1 (Format 55H. . . ) 

Contains any desired information to identify the particular case being run 
(card columns 1 to 55 only). This will be printed out at the top of a listing. 

Card 2 (Format 4E10.3) 

R radius of collector, R, cm 

BO magnetic flux density, Bq, G 

DTHTD 6 step size used in Runge-Kutta integration (spherical coordinate), 

A6, deg 

THl cone angle, deg 

Card 3(Format 8E10. 3) 

ZIO initial current^, !«, A 

1 ^ 2 

ZJO current density , Jq, A/cm 

VO initial electric potential, Vq, V 

DELR r step size used in Runge-Kutta integration (cylindrical coordi- 
nates), Ar 

CONFLO indicates confined flow when equal to 1 and Brillouin flow when equal 
to 0 

RAD radius of collector entry hole, a, cm 

FLOOR nvunber of iterations for space charge calculations (If no space 
charge calculation is desired, this value should be set to 1. ) 

Card 4 (Format 4012) 

NT number of trajectories to be calculated 

KPOI printout frequency control, (e. g. , 4 means to print results at every 
fourth integration step of a trajectory) 

Nl"| 

N2 > integrate trajectories N1 to N2 in steps of N3 

N3j 


^These two values are used only in the calculation of the source term in Poisson’s 
equation when space charge iteration is necessary. 
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IBAR when nonzero, indicates boundary values will be supplied along the 
negative polar axis 

LCX)P counter for nxunber of space charge loops 
KPOSC print control for interim space charge calculations 

Card 5 (Format 4F5. 1) 

These values are used by the integration subroutine in the spherical coordinate 
system, which is the system used in the beginning portion of the program. 

step- size controls; increase step size DTK by a factor of THSFAC 
every THSDEL (deg) in the range from THSWHI (deg) to 
THSWLO (deg) 

Card 6 (Format 16F5. 3) 

VSMI(I), electron energy classes for each of NT trajectories to be solved 

1 = 1, NT for 

Card 7 (Format 12, 10F5. 3) 

This calculation is done only once, after convergence has been achieved in the 
main part of the program . 

NEQL number of equipotential lines to be calculated 

EQL(I), array of potential values for each equipotential line 

I = 1, NEQL 

Card 8 (Format 16F5.3) 

XT(J, 16), boundary values on surface of spherical collector from initial 
J = 1, 16 cone angle 9 ^ to 180°, V 

Card 9 (Format 16F5.3) 

XT(J, 16) continuation of card 8 
J = 17, JF 

Two sample test runs are given in appendix A. These runs were chosen to more 
clearly illustrate the makeup of the input data deck and to show the various options in 
output format available to the user. 


THSFAC 

THSDEL 

THSWHI 

THSWLO 


CONCLUDING REMARKS 

We have described the numerical technique used to obtain a solution to a problem 
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associated with the design of an axisymmetric, spherical, depressed collector. Vary- 
ing initial conditions include current density of electrons, surface potential on the col- 
lector, collector geometry, and the velocity and position of entering electron groups. 

Of particular importance to the successful solution is the scheme used near the 
entry hole in the collector. Here the entering electron stream is of highest velocity and 
greatest density, giving rise to computational problems, especially in cases of high 
initial current densities. 

Calculations by others in the past were successful in handling current densities to 
50 A/cm before numerical difficulties precluded further computation. The numerical 
approach described in this report is used to solve systems with current densities to 
lOOA/cm^. 

Computational stability was achieved primarily by using a graduated mesh in both 
the r and 0 directions, with greatest mesh point density in the hole region. In addi- 
tion, computational efficiency was improved by taking advantage of the fact that all en- 
tering electron groups were distributed uniformly both in positioning and velocity group- 
ing. That is, the fastest group entered closest to the axis, and the slowest group entered 
nearest the edge of the entry hole. 

Space charge effects were found to be negligible for initial current densities less 

O 

than 10 A/cm . Consequently, the problem was reduced to requiring only one iteration 
and essentially solving Laplace's equation to obtain the potential distribution inside the 
collector. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 2, 1972, 

111-05. 
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APPENDIX A 


SAMPlf PROBlfMS 

Short, partial, representative computer output listings of two sample runs are 
presented. Run 1, depicted graphically in figure 6, is for a 0.2 spike, no space charge 
iterations, confined flow, and a cone angle of 60°. Run 2, figure 7, is for no spike, 
space charge effects calculated, Brillouin flow, and a cone angle of 45°. 

Figure 7 shows a poor pattern. Beam spread and separation is rather nonuniform 
and erratic. In contrast, figure 6 shows well- shaped and well-behaved beams with good 
spread and uniform spacing. 



Figure 6. - Sample run 1. No iterations necessary (current 
density, 8A/cm^); 0.2 spoke; cone angle, 60^. 



hole 

Figure 7. - Sample run 2. Five iterations for space charge 
effect; no spike; cone angle, 60° 
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Sample Run 1 


SAMPIE RUN N(]« 1 

NO SfACF rHARUF EFFECTS 

r.CNEfNEn ELflto. 

fKPUT> 

PO = H.RhOE-14 KP = 6.1S8E-07 

CMFC«= 1.4l4FflO ETAC> U753F«15 
in« 4.00nF>01 

JCb 7.0‘)CF«00 R 3 8.OC0E+00 

VO- 7.500F«’03 RO - l.OOOE-06 

TAUOP- r..24 3F-17 RAO ° O.^nOF^^LO 

CUIPLT EVERY 4 STEFS. 

INITIAL VALOFS- 


I 

UFL 1 

DEL 1 

theta 

VSH 

VI 

INJEC. 

PSIO 


1 RAU) 

lOEGI 

IDEGI 



ANGLES 


7 

u.r 06 

n.358 

179.642 

0,900 

6.75E»03 

2.000 

3.925E-09 

3 

0.013 

U.716 

179.284 

C.BOO 

6.00E«’03 

3.C03 

1.568E-08 

U 

0 .f . 1 9 

1.C74 

178,926 

0,700 

S.25E«^03 

4.000 

3.517E-C8 

5 

0,075 

1.433 

17B.567 

C.6P0 

4.5CE*03 

5.C00 

6.225E-C8 

6 

0.031 

1.751 

178.209 

L.500 

3.75E«>03 

6. COO 

9.659 E-C8 

7 

0.0 38 

7.149 

177.851 

C.40C 

3.0CE4-03 

7.C00 

l,373E-07 

fl 

I). 044 

2.5C7 

177.493 

0.20C 

1*5CE«^03 

8, COO 

1.809E-07 


EAR VALUES' 

R V 

0. 'l.SOU 

C.200 '1.500 

C. 310 '0.7fl? 

0.3SI3 '0.6P5 

0.4P0 '0.503 

O.SoO '0.411 
0.633 '0.377 

C.700 '0.760 

C.76C 'O.IQR 
0.R13 '0.148 

C.B60 'O.lUH 
C.900 -0.076 

0.933 '0.04H 

0.960 -0.078 

0.980 '0.014 

0.993 'O.OflS 

1.000 O. 







POTENTIAL 

FIELD 






INPliT- 






RHO 







tefta 

C.560 

C.633 

0. 700 0. 76C 

0. 

813 

0.860 

0.90C 

0.933 

C.960 

0.980 

fl.993 

l.flflO 

6C.U0 

-l.SOE^OO 

-1.5aE«'00 -l.SGE^OO -l.SOE^OO -1. 

50E*00 

-l.50E<>C0 -1.50E<“00 

-1.50E+00 

-1.50E+00 -1.53E*00 

-1.50E*03 

-1.50E*“0 

75.U0 

0. 

0. 

0. n. 

0. 


0. 

C. 

0. 

C. 

3, 

n. 

-l.25E>flfl 

90.00 

0, 

c. 

0. 0. 

0. 


G. 

0. 

C. 

C. 

fl. 

fl. 


1C5.00 

0. 

0. 

0* 0. 

0. 


0. 

c. 

0. 

0. 

0. 

r. 

-7.50E-fll 

I 20 .cn 

0. 

c. 

0. 0. 

0. 


C. 

n. 

0. 

c. 

0. 


-5.?0E-fll 

135. CO 

0. 

c. 

0. 0. 

0. 


0. 

c. 

c. 

0. 

0. 

fl. 

-2.5CE-fll 

150.0C 

0. 

c. 

0. L. 

0. 


0. 

0. 

0. 

c. 

0. 

n. 

-fl. 

165. CO 

0. 

c. 

0. Cl. 

0. 


0. 

c. 

0, 

e. 

0. 


0. 

1 76 .( 0 

0. 

c. 

0. L. 

0. 


0. 

fl. 

0. 

0. 

fl. 

hi 

0, 

176.75 

C. 

c. 

0. n. 

0. 


0. 

c. 

0. 

0. 

0. 

0. 

3, 

176.50 

0. 

c. 

n. 0 . 

c. 


c. 

0. 

0. 

c. 

3. 

0. 


176. T5 

0. 

0. 

0. 0. 

0. 


0. 

G. 

0. 

0. 

0. 

0, 

•fl. 

177.00 

f - 

0. 

0. c. 

n. 


0. 

0, 

0, 

0. 

fl. 

0. 


177-75 

0. 

c. 

0. 0. 

0. 


0. 

0. 

fl. 

c. 

0, 

0, 


1 77 . 50 

o. 

fl. 

0 . c. 

0. 


0. 

c. 

0. 

c. 

0. 

0. 

C. 

177.75 

0. 

c. 

0. 0. 

0. 


0. 

c. 

C', 

0, 

0. 

c. 


178.00 

n. 

c. 

0. 0. 

0. 


0. 

0. 

0. 

G. 

fl. 

0. 

0, 

178.75 

0. 

c. 

0. 0. 

0. 


0. 

0. 

0, 

0. 

fl. 

r , 

0. 

178.50 

0. 

c. 

0, 0. 

0, 


c. 

0. 

0. 

0, 


0, 

0, 

178.75 

0. 

c. 

0. r>. 

0. 


fl. 

0. 

0. 

c. 

hi 

0. 

0. 

179.00 

0. 

c. 

0. n. 

0. 


c. 

c. 

c. 

0. 

?. 

c. 

fl. 

179.25 

0. 

c. 

0. 0. 

0. 


0. 

0. 

c. 

fl. 

fl. 

0, 

7, 

179.50 

c. 

c. 

0. 0. 

0. 


0. 

0. 

0. 

fl. 

0. 

fl. 

fl. 

179.75 

0, 

c. 

C. 0. 

0. 


0. 

c. 

0. ■ 

c. 

0, 

0, 

A 

180.00 

-4.11F-01 

-3.27E-01 -2.60E-01 -1.98E- 

01 -1. 

48E-01 

-1.08E 

-Cl -7.57E-02 

-4.83E-C2 

-2.8CE-02 -1.40E-C2 

-4.60B-fl3 

0, 

*l«-1.50 

















POTENTIAL 

FIELD 






f UTPUT- 






RHO 







ThFTA 

0.200 0 

.3C0 0.393 

C A ^ BO ^ 

• 63 3 

0, 700 

0.760 

0.813 0.860 

0.900 0 

,933 0 

,960 0.980 C 

.993 l.flflfl 

60.00- 

1. 5000-1. 

5000-1.5000 

-IaSOOG-IaSI^OO-Ia 

5000-1 

,5000- 

1.5003- 

1.500P-1.5030- 

1,5000-1. 

5C00-1. 

5000-1.5033-1, 

50CC-1.50 

?c 


75.00- 1.38 34-1.35 15' 1.3773-1 .3C 86-1 .2939-1 .2 824-1 .2734-1.26 65-1 .2614-1. 7575-1.2548-1, 2 528-1 .2515-1 . 2507- 1.25C2-1 .25''? 

90.00- 1.7819-1.2208-1.17 39-1.1368-1 .107 1-1 .083 2-1 .0638- 1.048 l-l . C355-1 .C254-1.? 175- 1, 0113“ I .0066-1 ,0032- 1. O'* 1 1-1 . ?0 ''C 

105.00- 1.1946-1. l0 58-l.n366-0.98O7-0.9347-O.896o-0.8648-O,8384-0.8l65-0.7984-D,783 7-0.7 72C-0.7630-0,7564-C'.7521-f*.75CC 
120 . 00-1.17 7 7-1.7 C74-C.916 1-0.840 8-0.7774-0. 7233-0 .6772-0,6380-3 . 604 7-0. 5 7 69-0. 5 5 3 9- 0 . 5 3 54- 0 . 52 1?-0 , 5l?4-0. 5C34-?. 5? '*? 
1 35. CC-l. 06 90-0 . 9 7 79- 0. 8 1 57-0 . 72 16-U. 6404-0, 569 5-0. 5074-0. 45 30-0. 4058-0. 3654-0. 3 31 5-0. 3C37-C. 2820-0. 2659-0. 2553-0. 2500 

150.00- 1.0390-0.8706-0.7399-0.6301-0.5342-0, 4487-0,3716-0.3018-3.2388-0,1823-0.1326-0.0898-0.0547-0.0276-0.0093-0, 

165. C0-1.047C-0. 83 78-0, 6917-0. 5727-u. 470 8-0,3819-0.3042-0. 2369-0. 1793-0. 1308-0. 0910-C. C591-0. 0347-0. 0170-0. 0056 0. 
176.C0-1.1 701-0. E2 51-0. 669 2-0. 5474-0. 4454-0. 3573-0.2816-C, 2 169-0. 162 5-5. 1176-0,3813-0.0 526-0. C3C 8-0. CIS l-C. CO 50 , 

176. 7 5- 1. 1 75 1-0. 8 7^6-0,6685-0. 5467-0, 4448-0.3 567-0. 281 1-0. 2 165-0. 162 2-0. 1174-0.081 l-0.C525-O,03«7-'?, 015 1-0.0050 0. 
176.50-1. 1304-0.82 <0-0.667 7-0. 5459-0.444 1-0. 3562-0. 2807-0.2 161-0. 1619-0,1 172-0.0810-0,0524-0.0307-5,0 150-0. 0C49 0, 

176. 75-1. 1362-0.8233-0.6670-0,5451-0.4435-0. 3556-0. 2802-0. 2157-3. 1616-0. 1169-0. 3808-0. 0523-0.0306-0, 0150-0. 0049 <*, 

177.00- 1. 1475-0.8226-0.6661-0.5443-0.4428-0. 3550-0.2797-0.2153-0.1613-0.1167-0.383 7-0. 0522-0. 0305-0. •0150-0, 0049 ?. 

177. 25- 1.149 3-0. 82 18-0,6653-0.5434-0.4421-0, 354 3-0. 279 2-0. 2 149-0. 1609-0. 1 165-0,3805-0.052 1-0 ,0305-0. C149-0,0<'49 

1 77.50-1.1566-11.82 10-0.6643-0. 5425-0.441 3-0. 3536-0 .2787-0.2145-0. 1606-0. 116 2-0.3833-C.052C-C.03C4-'>, 0149-0,0049 0, 

177.75- 1. 165 2-0. 82Cn-o, 6633-0. 54 14-0.4404-0. 3 529-0. 27B2-C. 2 143-0. 16C2-C. 1160-0. 380 2-0.0519-0,0303-0.0149-0.0049 0. 

178. rO-l. 1745-0.8150-0.6622-0. 5403-0.439 5-0,3521-0.2776-0,2135-0. 1598-C, 1157-0,3800-0.0517-0 .0303-0,0 148-0,0049 
17«- 75- 1.185 1-0.81 78-0. 6609-0.5350-0. 4385-0. 3512-0.2 769-0. 2129-3. 1594-0. 1154-0, 0798-0. 0516-0 .030 2-0. 0148-0. CO 49 

1 78. 50- 1,1 574-0. a 164-0. 6594-0. 537 5-0. 437 4-0. 35C 2-0 ,2762-0. 21 23-0. 1589-0. 1150-0,0 796-C. 0514-0. 0301-0. 0148-0. 0''49 


20 



17fl.7S-1.21?O-O.ai4B-O.6577-r,53S8-(J.436l-O.3<»9O-O.2754-C,2ll6-O.l503-C.ll<^7-?.?794-C.C'513-^*"3'??-?,vl47-C*C''^3 
I 79. 110- 8128-0.6557-0. 5338-0. 9345-0. 3976-0. 2 7<»4-0.2Ua-0. 1577-0. 1 142-0. 3 791-C, 05 ll“?,C299-^,^ 147-''. 

\ 79. 25- 1.25 2»-0. 8102-P. 6530-0. 53 11-0.432^-0.3959-0. 2 73 l-C. 2097-0. 1568-C.l 137-0. 0788-C.C5C 8-C .02 97-0. ?! 146-0. 48 

I 79. 50-1. 2850-0. 80^5-0.6494-0. 5274-0. 4296-0. 3434-0. 2 714-C. 208 2-0. 1557-0.1 129-0. D784-C.0 504-0. C295-:.?l45-C.Cr 43 
1 79. 75- I -338 7-0. 8C04-0.6433-C. 521 3-0. 4245-0. 3393-0 .268 5-0. 2058-0.1537-0.111 6-0. 0777-C.C498-0.C29 I-"*, ''144-C.C047 
I 80. CO- 1. 500C-0. 78 20-0.O25O-0. 5213-0. 4249-C. 3393-0. 2685-0. 2058-0. 15 37-C.1116-C.0 777-C.C498-C.C29 1-0. C144-C.CC 47 


TABLfc- OF HHO-VflLUES FOR 
EOUlPOfENTlAL LINES 


THFTA 
6G.00 
75.00 
90. OU 
IC5.C0 

-0.90 -0.80 

U. 0. 

0. 0. 

0. 0. 

0.627 0.856 

-0, 70 

0. 

0. 

c. 

0. 


-0.60 

0. 

0. 

0. 

0. 

-0.5C 

0. 

0. 

0. 

0. 

-0.4C 

3. 

0. 

0. 

0. 

-0.30 

0. 

c. 

0. 

c. 

-C.2C -0 

0 . c 

0 . c 

0 . c 

0 . c 

.10 

120.00 

135.00 

0.412 0.531 

0.323 0.408 

0,667 

C.5L1 


0.821 

0.oC2 

l.OCC 

0.7C8 

3. 

3.82C 

0. 

0.938 

0. 0 

0, P 

• - 

150.0C 

C.283 0.350 

C.425 


0.505 

0.589 

0.675 

0.762 

0.845 C 

925 

165. CC 

O.270 0.324 

C.388 


0.460 

0.537 

0.618 

0.7C4 

P.794 C 

891 

176.f.C 

0.275 0.315 

0.375 


0.443 

0.517 

0.598 

0.684 

0.777 C 

879 

176-25 

0.275 0.315 

C.374 


0.442 

0.517 

0.597 

0.683 

0.776 C 

879 

176.50 

0.275 0.314 

C.374 


0.442 

0.516 

0.597 

0.683 

0.776 0 

879 

176.75 

0.275 0.314 

0,374 


C.441 

0.516 

0.596 

C.682 

0.775 0 

879 

1 77 .00 

0.276 0.313 

C.373 


0.440 

0.515 

0.596 

C.662 

0.775 0 

879 

177-25 

0.276 0.313 

0.373 


0.440 

0,514 

0.595 

0.682 

G.775 P 

878 

177.50 

0.276 0.312 

0.372 


0,439 

0.514 

0.595 

0.661 

0.774 0 

878 

177.75 

0-277 0.312 

C.371 


0.438 

0.513 

0.594 

0.601 

0.774 C 

078 

178.00 

0.277 0.311 

0.371 


0.438 

0.512 

0.593 

0.680 

0.773 P 

878 

178-25 

0.278 0.311 

C, 370 


0,437 

0.511 

0.592 

0.679 

0.773 C 

877 

17R.50 

0.278 0.310 

0.369 


0.436 

0.510 

0.591 

0.679 

P.772 C 

077 

178-75 

0.279 0.3T9 

0.368 


0.434 

0,509 

0.590 

P.678 

0-772 0 

877 

1 79 ,0C 

0.279 0.308 

0.367 


0.433 

0.507 

0.589 

P.677 

0.771 0 

076 

179.25 

0.2ft0 0.306 

0.365 


C.431 

C.5C5 

0.587 

0-675 

0.770 C 

676 

179.50 

0.280 0.304 

0.363 


0.428 

0.502 

0.585 

C.673 

0.768 0 

875 

179.75 

U,7B1 0.300 

0.360 


0.424 

0.498 

0.581 

0.670 

0.766 0 

874 

IRO.OO 

nfLTHFTA=-( 

11.284 0.297 

1.0200 OEG. 

0.349 


G.414 

0.490 

0.581 

0.670 

0.766 0 

874 

THJ THFTA 
STFP 0 

RHO ■ VT 

OVORH 

UVDTH 

RH-PR 

RH-PPR 

THD2 

PHD2 

PSI 

2 179.64 

i.noo c. 

0. 

0. 


2.408 EfC-1 

1.1456*03 

6.3666*14 

1.4706*12 

3.9146-09 

3 179.2ft 

i.f.oo r. 

0. 

0. 


1.5316*01 

4.606E*02 

1.3966*15 

1. 5306*12 

1. 5636-08 

4 1 7ft. 93 

l.OOO -C. -0. 

C. 


1. 1226*01 

2.4616*02 

2,267E*15 

1.6676*12 

3.5076-08 

5 178-57 

l.OOO -0, -0. 

0. 


8.8426*00 

l.6206*C2 

3.1136*15 

1.8966*12 

6,206 6-08 

6 178-21 

l.nou -c, -0. 

0. 


7.290E*00 

1.023E*G2 

3.7946*15 

2.3506*12 

9.6256-08 

7 177.85 

l.OOO -0. -0. 

0. 


6.1966*00 

7.2366*01 

4.1726*15 

3.5286*12 

1. 3676-07 

8 177.49 

STFP 4 

1.000 -C. -0. 

0. 


5.3826*00 

4.7596*01 

2.7436*15 

7,7756*12 

1.7976-07 

2 179.56 

0.967 -2.4UE-02 

5.592E403 

0. 


2.317E+01 

5.3016*02 

6.6946*14 

2.7906*16 

1.0306-09 

3 179.20 

0.979 -1.554E-02 

5.671EfU3 

n. 


1.5096*01 

-2,8606*01 

1.4C9E*15 

2.2446*16 

9.3196-09 

4 17ft. £5 

0-984 -1.147F-02 

5.54664^03 

0. 


l.ll3E*0l 

-8.l48E*ri 

2.2676*15 

1.6606*16 

2.3556-08 

5 178.49 

0.988 -9.nS7E-03 

5.5736403 

u. 


a. 801 E *00 

-7.8926*01 

3.0956*15 

1.3526*16 

4,4136-08 

6 178.13 

0.990 -7.539E-03 

5.5956403 

0. 


7.278£*00 

-7.3666*01 

3.7506*15 

1.2576*16 

6. 9746-08 

7 177.77 

0.991 -6.441E-03 

5.6156403 

0. 


6.210E*00 

-6.6436*01 

4.0896*15 

1.4366*16 

9,6436-''8 

8 177-41 

STFP 8 

0-992 -5.646E-03 

5.633E4C3 

0. 


5.459E*03 

-5.1376*01 

2.593E*15 

2.3616*16 

1.1026-07 

2 179.48 

0.935 -4.8866-02 

5.8976403 

0. 


2.2C9E*01 

9.558E*C2 

7.1556*14 

3.6516*16 

7,8936-10 

3 179.12 

0.958 -3.16OE-02 

5.9496403 

0. 


1.504E*01 

1-5466*02 

1.3916*15 

4.7616*16 

4,9396-09 

4 178.77 

0.969 -2.3316-02 

5.7246403 

0. 


1.124E*01 

-4. 9316*01 

2.1B2E*15 

4.3216*16 

1. 4406-08 

5 178.41 

0.975 -1.843E-02 

5.7546403 

0. 


8.936E*P3 

-8.B28E*Cl 

2.957E*15 

3. 8276*16 

2.9296-08 

6 178.05 

0.980 -1.521E-02 

5.7796403 

c. 


7.400E*00 

-8.1996*01 

3.5746*15 

3.5526*16 

4.852E-08 

7 177.69 

0.983 -1.299E-02 

S.619E403 

0. 


6.295E*00 

-5.C47E+C1 

3.9166*15 

3.4666*16 

7.017E-C8 

8 177.33 

0.985 -1. 139E-02 

5.637E4U3 

u. 


5.4576*00 

1.168E*01 

2.5186*15 

3,5776*16 

8.9776-08 

SIFP 12 

2 179-40 

0.906 -7.391E-02 

6.293E403 

G. 


2.065E*Cl 

1.049E*C3 

7.9436*14 

2.9336*16 

4.1486-10 

3 179.04 

0-937 -4.816E-02 

5-959E403 

0. 


1.467E*01 

3,465E*02 

1.4296*15 

5.3116*16 

2.7516-09 

4 178.69 

0-953 -3,5586-02 

5.991E403 

0. 


1.122E*C1 

B.068E*01 

2.1526*15 

5.7746*16 

e. 8316-09 

5 178.33 

0.963 -2.8086-02 

5.7606403 

0. 


9.012E*0D 

-l.634E*01 

2.8606*15 

5.6186*16 

1.9336-08 

6 177.97 

0.969 -2.3246-02 

5.7846403 

0. 


7.4866*00 

-3.8946*01 

3.4366*15 

5.4256*16 

3,3586-08 

7 177-61 

0.974 -1.977E-02 

5.8046403 

0. 


6.348E*C0 

-2,49ie*Gl 

3.7366*15 

5.2606*16 

5.061E-08 

8 177.25 

0.977 -1. 7206-02 

5.8236403 

0. 


5.467E*00 

-1.826E+01 

2.4326*15 

4.5626*16 

7.3426-08 

STFP 16 

2 179-32 

O.H70 -5.R06E-02 

6.5156403 

0. 


1.921E*CI 

9.9796*02 

8.9076*14 

2.1946*16 

2.5876-10 

3 178.96 

0.917 -6.500E-02 

6.3166403 

0. 


1.4116*01 

4.3496*02 

1.5096*15 

4.9206*16 

1.6886-09 

4 178.61 

0-937 -4.806E-02 

5.9986403 

Q. 


I.IC4E*01 

1.755E*C2 

2.1026*15 

6. 1056*16 

5.6766-09 

5 178.25 

0.950 -3.8116-02 

6.0236403 

0. 


8.980E*00 

5,8C2E*C1 

2.8306*15 

6,4386*16 

1.3UE-08 

6 177-89 

0.959 -3. 1396-02 

6.C44E403 

c. 


7.506E*CO 

1.116E*01 

3.3616*15 

6.4796*16 

2.3796-08 

7 177.53 

0.965 -2.6666-02 

5. 8f 96403 

0. 


6.364E*00 

2.516E*00 

3. 7006*15 

6.4746*16 

3. 6446-08 

f 177.17 

0.969 -2.3156-02 

5.8276 403 

O. 


5.485E*00 

-3.021E*00 

2.3386*15 

6.1836*16 

5.0456-08 

STFP 20 

2 179.24 

0.852 -1.2126-01 

6.939E403 

c. 


1.7686*01 

9,2126*02 

9.9916*14 

1.6506*16 

1.8146-10 

3 17h-88 

0-898 -e. 1416-02 

6.5976403 

0. 


1.3496*01 

4.5786*02 

1.6156*15 

4,3006*16 

1.132E-P9 


1 

SKlIfh TO k-l COORCINATE SVSIEM- 
TBAJFCTOBY NO. 4 DELR».010 


TFETA 

RHO 

ft 2 

OVDR 

OVOZ 

2P 

ZPP 

V 

1 67.39 

0.388 

0.085 -0,379 

1.3466-01 

-1.5916*00 

1.8366-01 

-2.5206*02 

-6.9676-Cl 

1 66.59 

0.408 

0.095 -0.397 

1.0976-01 

-1.3886*00 

-3.8086*00 

-2.3846*02 

-6. 6846-01 

166.75 

0.457 

0.105 -0.445 

1.1536-01 

-1.3876*00 

-5,4946*00 

-1.1906*02 

-6.014E-C1 

167.20 

0.51fl 

0.115 -0.505 

9.9886-02 

-1.2846*00 

-6.4696*30 

-7,6246*01 

-5.207E-P1 

167.72 

0.586 

0.125 -C.573 

9.7606-02 

-1.2186*00 

-7.0836*00 

-5.0656*01 

-4.3426-01 

168.22 

0.660 

0.135 -0.646 

1.0156-01 

-1.1626*00 

-7.4886*00 

-3.2386*01 

-3.458E-CI 

168.67 

0,737 

0.145 -0.722 

1.1136-01 

-i.iiie*cr 

-7.7296*00 

-1.7236*91 

-2.5806-Cl 

169.06 

0.815 

0.155 -0.800 

1.1666-01 

-1.0146*00 

-7.8446*00 

-5.8956*00 

-1.730E-CI 

169.36 

0.894 

0.165 -0.879 

1.3596-01 

-9.6206-01 

-7.0406*00 

5.3326*00 

-9.392E-C2 

169.65 

0.973 

0.175 -0.957 

1.4286-01 

-8.4056-01 

-7.7546*00 

1.2056*01 

-2. 2766-02 

1 69.87 

1,050 

0.185 -1.034 

1.6276-01 

-8.0976-01 

-7. 6076*00 

1.7016*01 

4.1396-02 



SIFP 380 


? 

168.02 

0.27S 

-6.922F-01 

1.699E+04 

-2.550b»02 

-7.912E-02 

l.20lE*ri 

3.917E«>16 

2.703E+13 

4.377E-11 

.3 

167.66 

0.325 

-7.N53E-01 

l.ltiBE^CA 

4.398E^02 

1.515E-01 

2.C84E*0l 

1.513E+I6 

1.977E»14 

7.952E-11 

STEP 

2 

38A 

167.86 

0.2T5 

-f.<J15E-01 

l.652E*f.A 

-2.484E^02 

-1.124E-01 

1.163E^01 

3.953E+16 

2.554E413 

4. 502E-11 

3 

167.50 

0.325 

-7.960E-01 

1.187E*04 

4.388E«^02 

9.171E-C2 

2.2COE*01 

1.456E»16 

1. 887E«14 

8.117E-11 

STEP 

2 

388 

167.70 

0.275 

-e.9C6E-0l 

1.6£5E»C4 

-2.396E*02 

-1.452E-01 

l.lTCE^Ol 

3.9B2E^16 

2.411E^13 

4.634E-11 

3 

167.34 

0.325 

-7.964F-01 

1.186E^04 

4.383E*02 

2.838 E-02 

2.339E*Cl 

1.395E^16 

l.798E*14 

8.299E-11 


Swlir.H TO K-/ COUKiliNATE SYSTEM- 


TRAJECTORY NU. 
T6FTA RhO 

3 RELRe.OlO 
R 2 

OVDK 

ovoz 

ZP 

ZPP 

V 

167.26 

0.325 

0.072 -0.317 

1.733E-01 

-1.582E+00 

2.104E-01 

-2.410E*02 

-7.965E-CI 

166.18 

0.342 

0.082 -C.332 

1.506E-01 

-1.574E+0C 

-2.982E+00 

-1.719E*02 

-7.711E-C1 

166.03 

0.3B0 

0.092 -0.368 

1.73 7E-01 

-1.577E+0C 

-4.195E+0C 

-8.929E*ni 

-7.116E-CI 

166.2? 

0.427 

0.102 -0.414 

1.223E-01 

-l.386E*C0 

-4.956E*03 

-6.387E»C1 

-6.440E-C1 

166.55 

0.480 

0.112 -0.467 

1.CI04E-01 

-1.286E+00 

-5.507E*3& 

-4.997E^Cl 

-5.7ClE-ri 

1 66.94 

0.538 

0.122 -0.524 

1.125E-01 

-1»282E^0C 

-5.936E*0r 

-3.653E»0l 

-4.954E-CI 

167.32 

0.600 

0.132 -0.5E5 

1. 106E-01 

-1.216E *CC 

-6.26CE*JD 

-2.759E«01 

-4.186E-C1 

167.69 

0.664 

0.142 -0.649 

1.141E-C1 

-1.162E*0C 

-6.495E*D0 

-1.98 1E*0l 

-3.415E-01 

I68.l>? 

0.'731 

0.152 -0.715 

1.220E-01 

-1.112EH)C 

-6. 655 £♦ 00 

-1.270E+01 

-2.653E-C1 

168*3? 

C.799 

0.162 -0.762 

1.332E-01 

-1.C62E+OC 

-6.752E+0C 

-6.234E»0C 

-1.914E-01 

168.58 

0 . 867 

0-172 -0.850 

1.386E-01 

-9.669E-01 

-6.787E*0C 

-8.9646-01 

-1.2C9E-01 

166.80 

0.935 

0.18? -0.918 

1.477E-01 

-8.795E-01 

-6.772E*03 

3.81C-6^00 

-5.5656-02 

168.99 

1 .004 

0.192 -0.965 

1.6C1E-C1 

-8.159E-01 

-6.7156*D0 

7.4356»00 

3.C3CE-C3 


Sul Tf.H TO 

R-Z CnORCINATE SYSTEM- 





TRAJECTORY NO. 

2 RELR«.C10 






THETA 

RHQ 

R Z 

OVDR 

uvoz 

ZP 

ZPP 

V 

167.58 

0.276 

0.059 -0.269 

5.9C9E-01 

-2.I63E+00 

-3.479E-01 

-1.0756*02 

-8.8986-01 

165.95 

0.286 

0.069 -0.277 

5.171E-01 

-2.C83E*CC 

-1.0886*00 

-5.1816*01 

-8.6806-01 

164.71 

0.301 

0.079 -0.290 

8.156E-03 

-1. 6176*00 

-1.5296*30 

-4.2586*01 

-6.3716-01 

163.81 

0.320 

0.089 -0.3C8 

1.8 75E-02 

-1.6116*00 

-1.9496*00 

-4. 14 76*01 

-8.C89E-01 

163.21 

0.344 

0.099 -0.329 

2.0 8 56-02 

-1.608E*00 

-2.3596*00 

-4.0636*01 

-7.7406-01 

162.88 

0.371 

0.109 -0.355 

1.629E-02 

-l.609E*00 

-2.7646*30 

-4.0346*01 

-7.326E-C1 

162.76 

0.403 

0.119 -0.384 

-4.307E-02 

-1.4346*00 

-3.1666*00 

-4.0586*01 

-6.8646-01 

162.82 

0.438 

C.129 -0.418 

-3.963E-02 

-1.433E*00 

-3.5796*00 

-4.1576*01 

-6.3856-01 

163*01 

0.477 

C.139 -0.456 

-3.965E-02 

-1.4336*00 

-4.0016*00 

-4.2906*01 

-5.8466-01 

163.37 

0.520 

0.149 -0.498 

-6.120E-02 

-1.3406*00 

-4.4426*00 

-4.4956*01 

-5. 2036-01 

163-70 

0.568 

0.159 -0.545 

-7.246E-02 

-1.2806*00 

-4.9016*00 

-4.7236*01 

-4.6696-01 

16^-15 

0.620 

0.169 -0.596 

-6.173E-02 

-1.2756*00 

-5.3786*00 

-4.8286*01 

-4.C19E-01 

164.64 

0.677 

0-179 -0.653 

-5.701E-02 

-1.223E*00 

-5.8676*00 

-4.8666*01 

-3.3296-01 

165.14 

0.738 

0.189 -0.714 

1.827E-01 

-1.112E*0C 

-6.2706*00 

1.057E*0C 

-2.610E-C1 

165.60 

0.8fJl 

0.199 -0.776 

1.873E-01 

-1.0656*00 

-6.2526*00 

3.0056*00 

-1.914E-01 

165.99 

0.664 

0.209 -0.839 

1 .8486-01 

-9.747E-01 

-6.2156*30 

4.456E*0C 

-1.256E-C1 

166.31 

0-927 

0.219 -n.901 

1.954E-01 

-9.2676-01 

-6.1586*00 

6.4416*00 

-6.4716-02 

166.59 

C.989 

C.229 -0.962 

1.924E-C1 

-8.249E-01 

-6.0886*00 

7.4006*00 

-9.440E-03 

166.82 

1.050 

0.239 -1.022 

2.100E-01 

-8.192E-01 

.-6.0076*30 

8.7016*00 

4.229E-02 


Sample Run 2 


SAMPLE PUN NOr 2 
5 TTEPATIONS for SPACE 
CHAPCE EFFECTS 
BRILLOUIN FLOfc 
INPUT- 

EC = B^86CF-14 KP = lel60E-C7 

0MFGA= t341AE«^10 ETAO= 1,753E*15 

IG= 9o3CCF-02 

J0= R.;CCCE+Ol P = 8»00CE+0o 


V0= 

8«C0CF*O3 

RO =l.0OGF-06 

TAIIC2= 

0. 22PF-17 

RAO a 0:1006 * 00 

OUTPUT 

EVERY 6 STEPS? 



IN IT lAL 
I 

values- 

OF.LI 

DELI 

THETA 


(PAOl 

(OEOt 

(OEG) 

2 

0.,002 

0-.390 

i79a 910 

3 

P..003 

C.179 

1 79o 82i 

4 

0 UC5 

0?269 

179:731 

5 

0.,(j06 

0.358 

1 79c 642 

6 

.0.008 

0.443 

17Pt552 

7 

0oU09 

0. 53 7 

1 79.- 463 

e 

c 

0.62 7 

179,373 


•♦NO BAO*^ 


VSM 

VI 

INJEC. 

ANGLES 

PSIU 

Oc93 7 

7 20E+03 

2.0UO 

C. 

u«v8j3 

6.406*0 3 

3. 000 

* c . 

0.733 

5.606*3 3 

4.000 

0. 

Js6j j 

4*006*0 3 

5.000 

c. 

O.50J 

4., UOE+0 3 

6.000 

c. 

^i,:4Dj 

3., 20 E *0 3 

7. OUO 

0» 


1 . 60E+.J 3 . 

8. 00 0 

*<• 


22 



input- 






PCTEKTIAL field 
8HQ 






theta 

0r56C 

C r 633 

0: 700 

Q 760 

3 = 

313 0.-860 

Os 900 

0,933 

C.96G 

0.980 

0.993 

i.OOu 

4 5r CO 

-i. 00 6 *00 

-icOOE*t)0 

-i.ooe+co 

-A: 03E+C0 

-A.. 

)0Eh» 5 -t.OJE+uO 

-i.OOE^UO 

-i.LG6*C0 

-l.QLE*00 

-l,OOE*uO 

-1. OOE+OO 

-i.uOE*00 

60o CO 

0; 

0: 

0, 

0. 

9., 

Oo 

0, 

0. 

C. 

0. 

0* 

—8. oo E-oi 

7 5, CO 

0 

o. 

•J: 

0 

C: 

Jo 

Wo 

0. 

Co 

0* 

w. 

-6.0uE-vi 

9C- CO 

• : 

!h 

0.. 


0, 

Jo 

0* 

c. 

t* 

0. 

0* 

- 4 .UOE- 0 I 

I09r, CO 

0 

X 

0, 

c. 

0. 

Os 

0, 

0. 

Co 

0* 

v3. 

-2.trO£-Oi 

1 20. CG 

0, 

Or, 

0> 

0. 


1 ♦ 

0* 

0. 

Oo 

0. 

u. 

-0. 

)35 CC 

0. 


Oc 

vis 

0 . 


Os 

c. 

c« 

0. - 

u. 

U. 

150: CO 

0 

Oo 

0, 

0, 

0 • 

0, 

Js 

c. 

(io 

0. 

Oo 

U* 

16 5. CO 

0 

Ot 

0. 

0, 

0 , 

0. 

0, 

0. 

Co 

0. 

0. 

0. 

176 25 

0 

3, 

J. 

Os 

J-. 

0, 

0. 

c. 

c. 

0. 

0. 

0 . 

' 76* 50 

0 

‘J. 

0. 

0: 

0* 

0* 

0, 

0 • 

c. 

Co 

0. 

o. 

} 76r 75 

0 

0* 

D. 

0, 

0. 

0. 

0, 

0* 

0. 

0. 

0. 

0 . 

• 77, CO 

0 

f/. 

0- 

C‘> 

Os 

0. 

0, 

Co 

c. 

0. 

0. 

0. 

177,25 

c. 

0 

0. 

0. 

Os 

0, 

0, 

Vo 

c. 

0. 

0. 

0 . 

177.: 50 

0: 

Or 

0 

0~. 

0, 

0 . 

0, 

0. 

c. 

u. 

u. 

U. 

3 7 7. 75 

0 

Or 

Oi 

G. 

0: 

Oo 

0* 

G. 

1. 

0. 

0. 

Oo 

1 78: CC 

0. 

0* 

0: 

0. 

0: 

Oc 

0, 

0. 

Co 

0. 

Oo 

0 . 

178 25 

"l 

J, 

0.. 

Os 

0-. 

Ot 

Oc 

Co 

Cl. 

u. 

u. 

0. 

178.50 

0 

(5 

J: 

0: 

0* 

0. 

Oo 

0 . 

Cc 

0. 

0. 

u. 

I78a 7 5 

0: 

0, 

0. 

0: 

0* 

J. 

0, 

u. 

c. 

0. 

0. 

0. 

’ 79 CO 

0 

<3c 

(h 

Us 

0, 

Oc 

0. 

c. 

c. 

Oo 

o. 

u. 

.179,25 

0 

0 

0, 

O'. 

0. 

t>. 

u. 

0. 

c. 

c. 

J. 

Oo 

> 79*50 

0 

iK 

0: 

0, 

0. 

Oo 

Os 

0. 

0. 

0. 

Oo 

0 . 

?7Q 75 

0 

•c 

0: 

0, 

0, 

•3 « 

Co 

Co 

c« 

0. 

u. 

0. 

1 80c. 00 0 . 

Tl=-? .0) 

0. 

a 

Os 

Os 

J, 

0, 

Co 

0. 

u. 

0. 

0 . 


PHTEWTIAL field 

OUTPUT- RHC 

THETA »;c20C -J, Ov^B-J 0-560 0,633 3. 700 Or 760 0>S13 0.860 O.SOO 0*932 C.960 0.980 C.993 1.000 

A5, OC-1 OCOC-y.:OOOC-l OOOO-l- 0300-i.000o-L.05j0-i>000U-l.Jt&b-lv0u00“io0000-i^0au0-l,C0CC-i.u00v-l.0 jUO-1. JJuO-A.CiOJO 

60.00- Co 8830-C-6581-C, 8402-0 8269-0» SiTj-jr 8 j98-35 8047-0 , 8l 1 4-0. 7994-0. 7985-OV 1982-C, 7S 84- C . 7989-0. 7994-0* 7993-0, 8000 

75.00- 0. 781>-0c 7^^.3-l; 7074-0: 6813-0* 661 u-J. 6452-0 6330-0 6236-5 w 6i65-v . 61 12-0* 6C72-C. 604 5- U. 602 5-0, 6C12-0. 6OO4-0. 6000 
90-.00‘G T07 3-P. 64.*>fc-o. 5938-0. 9547-t* 522 8 -i. . 4965-3. 4747-0.,4568-0 .4422-0, 4303-0.42t 9-C, 41 3 5-C . 4079-O. 4039-tr* 4Ui 3-0. 40oo 

\0'.00-0 64.07-0^5600-0.4977-0. 4463-0.4026-.y, 3647->. 33iR-0 30 32-J , 2786-0, 2577-0> 2404-L, 226 5-0 . 2i 56-0. 2 077-u. 2 025 -0. 2000 
).?O.CO-t - 58M»-.l *934-0,4206-0^ 3 5 98-Oc3C69-Oi2597-jj 2169-0 .1777-5. i4i 9-0, 1092-0* 0 80U-0, 0546- C ,0334-0, Oi 7C-0, oo57-0. 

J35, 00- 0 545 2- 4430-0.- 3 6*1- 0:2 9Q7-C, 2443 -J, 197 i-)* t 564-u,.l 214-J ,0917-0.06 68-0, C464- C,0 3C 2- 0 4.- ol 7 7-Ue 0C87-w* Oo29 0. 

150.C0-C 516(:-ur 4083-0. 3266-0, ? 6l i -0. 2 o7S-jj i631-J. 1 263-5.-.0961-J ,07i 3-0, 0513-0.0353-0,0228-0. Oi 33-0. 0J65-U. 002 1 u. 

• 65, 00-0 498t-6* 3892-0. 3055-Ot 2406-u, , 884-0, 1463-J, I 122-J .0 847-5 ^062 5-0 . 0447-C. C30 7-C.C198-0. Oil 6-0, 0057-U. 00i9 0. 

3.7t„ 25-0 4932-0*382t-0o299i-0. 7 3 44-0*1829-5, i4i 5-3, IU82-0 . O815-U* 0600-0,0430-0* C 295- C.C19C-C.X;; 11-0*0054- 0.001 8 O. 

1.76. 50-0 49 3 i- 0.3819-0- 2991-0: 2 343 -U, * 828-0- l4i 5 -0 . 1O82-0 5 8 15-5 .0650-u. 0429-U. 0295-C,Ci 9C-0. OU 1-0. 0J54-0, OOA 8 0, 

:.,76(. 75- Jr 4931-i.r 3819-0- 299J-0, 2 348-y, 1 829-0.: 14i.4-j. lo 32-0-,0 8 14-0 ,u60o-O. 0429-0. 0294-C,ol 9C-C . .yl 1 1-0. 0054-u, Ujl8 0, 

177, 00-0,493)-0t 39': 9-0,299J- J 2 343-0. i 827-0^4 14-5, iO 8i-0 .0814-0 ,06Ju-U. t429-0. 0294- C. Ll9C-tl . Oi i i-o. 0054-0, Cv,l8 u, 

177. 25-0- 4931-u 38/ 8-0>2989-0* 23 42-0,1827-0. l4* 3-o, 1 0 Sl-C 5814-0,0600-0, o429-0*u294-t.0i9C-0. Oil 1-0.0054-0, Uul8 U. 

1.77. 50-0. 493C-(,«38' 7-C, 2989-0-. 234?-u, i 026-Jo I4i3-J , lo8i-0 . J8 14-0 ,O60u-t‘. 0429-C. C 294-C.5i 9t-i. wl I i-0. 0 J54-U, UOl 8 O. 

177. 75-0.493C-f\ 391 7-0 : 2989-0, 2 341-0. i 82 6-0o 14^ 3-3* lO 9>.-0 .0013-0 «o599-u.0429-i;* 0294-0, Li9C-0 , Ui i .-'J.0J54-0. OOl 8 o, 
•■7P..0C-0 493C-0. 38;. 7-0*2 99 9-0. 2 341 -Or. *026-0-. 1413-3.. 10 90-0 .08 1 3-0 *0599-0. 04 29-G* <J294-C,L iSC-O.Oi i i-O. 0054-0. J'Jl 8 0. 

1700 25-0-. 4P2V0 38*6-0, 2988-0, 234i -0. 1 826-0* 14i 2-3, lO 8o-0 , j 813-3 .0599-0, 0429-^,. C294-0.0 19C-0.01 1 i-i.. 0054-O, JOi 0 O. 

17fi^ 50-*jo4O25-G*39\6-0r:2g88-a,234i-U, 1825-0, 14*2-0, 1080-0.0013-5 ,0599-0,0429-0, C294-0. Ui 89-0, Ul i 1-0. 0J54-0. OOl 8 0, 

1 78, T5-0 40?5-C 3916-0^2988-0. 2 340-u, I 825-0,141 2-3* lU9d-0 . 0 8 1 3-J ,0599-0. 0429-0, 0294-0. LI 89-0, ul i i-O. 0J54-U= 00l8 O. 

179,00-0 4929-C. 33a 6-0, ? 997- -J, >343-0. 1 825-Oc 14 1 2-3 r 10 8u-C J8 13-J *0599-0, 0428-C, 0294-0.0189-0. di 1 1-%-. l J54-o, OJ l 8 0. 

\7P, 25-0^492 9-0,3916-0^2 987-0, 2 340-0. i 82 5-0. 141 2-0, iO 3J-0*j 813-3 *0599-o. 0428-0. C294-L,oi 89-0. di 1 1-0. 0L54-0. Ooi 0 u, 

179: 50-0: 2 9- c- 38 i 6-0. 2 98''- 0, 2340-0^ * 825-Or I4i 2-3.. iO 80-0 ..u 8 1 3-0 .0599-0* 0428-0. C2S4-0 189-0. ol 1 l-o. 0054-O, Oui 0 0. 

179 75-0 4929-c 3«i6-C* 2987- C, 2340-0* i 825-0. 14i 2-J., 10 Su-O : 08 13-0 ,u599-u.0428-C. C294- C.O J S9-0. 3i 1 1-0. 0054-0. uOl 3 0* 
180=00-0 4929-0, 38.i6-0. 2987-0.. 2340-0.1825-0*1412-3: 10 0U-O* U81 3-J *0 599-0.0428-0. 0294-C.0 1 89-U,Ui 1 i-O. 0j54-O. 0018 O, 

PASS ‘ PHLICWS- 

OfLThETAs-r. C2CC DEC: 


TP J 

theta 

Pnn 

VT 

0VD8H 

dvdth 

RH-PP 

STEP 


r, 






2 

no 

s9t. 

*:CCO 

Co 

0 

Os 

2«715E*01 

3 

i70 

82 

1 OLO 

Or 

0 

0. 

l,789E*-0i 

4 

7 79 

73 

.T COf. 

■\ 

-3= 

u. 

1.334E+3A 

5 

179 

.64 

l,CGO 

Oc 


J 

1* J62EtJ; 

6 

179 

s55 

1 sCrA' 


0. 

J*. 

8c82lE*JL 

7 

?79 

s 46 

\ ou 

G. 

0. 

0- 

7.538EVJO 

8 

179 

,37 

J.,( ':c 

0.: 

0: 

O' 

6c5i^6E+jC 

STEP 


6 






5 

370 

5? 

O 1.979 

-5- 829C-03 

2s2 5-^F*o3 

w. 

9,917E+iu 

STEP 


*2 






5 

379 

40 

0.,q58 

-1.= 158F-02 

2- 363F*03 

0, 

9,46oF*:u: 

STEP 


< 8 






5 

179 

>8 

0«S3o 

-7r 730F-0? 

2, 364F+U3 

c. 

9r,u35 E*JU 

STEP 








5 

179 

.1 6 

C ,S 20 

-’.299E-02 

2.. 5C7E1U3 

Os 

8c 638 E *70 

STEP 


3'? 






5 

179 

. i4 

C .St3 

-2r 854F-C2 

2.507F«-O3 

C’s 

0.'266E*JC 

STEP 


?6 






5 

178 

> 92 

0,. 886 

-3; 418F-C2 

2s 691EA03 

u. 

7.9i7E*UC 

step 


42 






5 

3 78 

80 

Or 873 

-3- 864E-02 

2,-.69Le*q3 

Os 

7.590E+00 

STEP 


A 8 






5 

178 

S» 

Or. 854 

-4. 5c6F-Or 

2* 923Ef 03 

Os 

7.283E*00 

step 


54 






5 

3 78 

. 56 

0,8 30 

-5...C53E-02 

3, 923F *03 

Os 

6o9Q4F+0U 

STEP 


60 






5 

: 78 

44 

Oa825 

-5.. 5>8E-02 

2= 923T^U3 

0* 

6o722E*00 

1 77 ., 

25 

0.-. 


. 65E*C: 

2s22F + Ol 3. 

;%5E*0l 

4c5 5E4Ji 5. 

:.77r 

50 

0 • 


A. 57E+ 31 

2e7CEEii 2. 

97E+01 

4;34E*0l 7. 

?.77^ 

75 

O'- 


I , 48E*C1 

:* 97E*o: 2. 

79E+0t 

4jJ6E+01 6 c 

• Tfl 

00 

0.: 


it 37FA0’. 

)s34g*0l 2.6wE*Ul 

3,796+01 ho 

)7P, 

75 

Jv 


■•* 26E* -7i 

J * 7C'F*Ci 2» 

4GE+01 

3=53E+01 5c 

I78r 

5C 

0. 


It 12F4 0’: 

n54e*Li. 2, 

LuE+Oi 

3,23E+0i 5, 

178 

75 

0 


9,. 77F.*00 

l:35E*Ul 1, 

96E+01 

2: 92E+U1 4. 

1 7Q« 

»0 

0, 


8» l6E*vJ 

1 14E*f.i Is 

67E*Ci 

2:55E+0i 4, 

}70 

25 

V 


6t >9 EaOO 

8: 93F + C0 1. 34F*0l 

2.J6E+J1 3c 

1 70. 

50 

c 


4, iOE*,)o 

6: JOETUO 9, 

24E+00 

»,49E+oi 2r 

1 79, 

75 

J 


2- lOE+OO 

3. ClE*Ou 4, 

62E+0J 

7.s53E+U0 1. 

l8C. 

CO 

0 ■ 


C-. 

Or G« 


s 0 c 


PH-PPR 

TFD2 

PHD2 

PSI 

3s401E+03 
1-881E+03 
is 236E+U3 
9* ii>iE+U2 
7,448 E+va 2 
6.55flE+w2 
3,852E+u2 

5,3446+14 

1*0926+15 

1.7156+15 

2.3106+15 

2.70CE+15 

3.0316+15 

1.9816+15 

*.8818+17 

1.8786+17 

i,87*e+i7 

i.867E+17 

1.855E+17 

i.831E+i‘i 

1.77CE+i7 

2.4286-10 
9*7076-10 
2,lB2E-09 
3.871E-L9 
io O2 9E-09 
0* 625E-C9 
1.154E-08 

2. 274E+02 

2.6226+15 

5,1566+14 

3.473E-10 

2, io5E+02 

2.8526+15 

1.7476+13 

9.595E-li 

1.95BE+02 

3,0956+15 

2,t646+ 32 

4. 5646-U 

l*835E+w2 

3.2516+15 

4,4E2E+ll 

2.7846-11 

1.713E+U2 

3.6226+15 

l*385E+ll 

i.946E-ll 

lo 613E+02 

3.9C7E+i5 

5,282E+1C 

i,470E-ii 

i.5o9E+o2 

4,2C76+15 

2,4526+lC 

1*1886-11 

i. 4286+0 2 

4.5226+15 

i,256E + K 

9.922E-12 

Is 339E+02 

4,8536+15 

7,C37E+l9 

fl*5J3E-12 

4.253E+02 

5.2LC6+15 

4,2326+c9 

7, 5o4E-l2 


3GF+01 

4*47E + C1 

c • 

Lc 

Co 

J. 

0. 

326+01 

i,05E+o2 

1.24E+C2 

C. 

c. 

0, 

0 . 

56F+01 

8,736+Oi 

6.596+Ci 

c. 

0. 

o. 

0. 

02E + 0i 

7* 58E+01 

7.49F+C1 

c. 

Vo 

Jo 

0. 

56F+01 

9c97E+0l 

1.576+C2 

2,2GE+02 

Co 

j* 

u. 

136+Cl 

a.90E + Ol 

1.886+02 

1.536+02 

Vo 

u. 

0 . 

65E+G1 

8* i4E+01 

i . 6 lE +G 2 

2*846+02 

Co 

0 . 

0. 

136+01 

7, 30 E +01 

1.41E+C2 

3. 1 2E +0 2 

3.246 + 02 

u. 

w. 

48E+0i 

6,266+01 

1*236+C2 

2,63E+02 

6.33E + J2 

8. 8 8E+02 

0. 

57E+Ci 

4s85E+Gl 

9.906+Ci 

2,1 66+o2 

5«0iE+O2 

1*176+03 

0. 

35E+oi 

2.716 + 01 

t.OLE+C* 

1. 43E+o2 

3. 566 + 02 

8.476+02 

o. 


Oo 

C. 

C, 

C. 

0, 

w. 


23 



PCTENT rCL F IELD 

OUTPUT- RHG 

THPTft O-20C O-3<:0 ui.393 Ct 480 0:550 0^63 3 3ir700 0.760 0*813 0. 06C C.5CC 0^93 2 Cj 960 C* 980 C*993 i.ooO 

45*CC“ 1- !tO< C- /.«• COUO- 1 5 0J0*>- 1, OOOo-i > 0000“1'- 0COO“i » COOC “l';OC(00-l oOOOO-ia OOUO- 1 •00CO-i» COt'C- 1 •oOOU- la OUOO-i a CUO'j'~i* 0000 
60«00-0- 8R5 2-G-. 96i)l-0.‘ 841 ’7-0^ 82 83-0, 8i8i -Of. Si‘J7-j, SO54-U.0ui9-O « 799 8-Og 79 87-0* 19 84-C* 79 E5-0o 7989- J. 7994-0* 7998-0. SuOU 
75< CO-0 795 7452-0. 7103-0, 6342-C, 6 634-0., 6471 -0, 6344-0 -.624 7-0 o6i73-0*6117-C* 6C77- C*6C47-0«6027-0. 6Ci3-u* 6004-O. 6000 

o-O. UC-Cr 7,.34-Oc 6487-0, 8992--X 5594-0« 5267- Jc 4996-j, 4771-C ,4596-0o4435-0.4313-0c 421 5-C*4i29-C. 4081-0. 4040-U.4Ul3-G.40ou 
)0‘.00-0 649C-C. 5687-0 505P-0, 4536-0, 4 C37 -0 , 369 7 -3; 3358-0 /3o63-0o23o9-0. 2594-Cj 24i6-C.2272-C.2i6i-0. 2079-0. 2U6-0. 2000 
) iC:. C0~O 59 7.5-0. 505 5-.5c 4 329-0, .8 7i 4-0-,3i72-u. 2684-0, 3240-0 ,18 33-0.1461-0. ii23-O.CS22-C.(;5e(-G.i/ 342-0. 01 74-0. OUSfl-u. 
l3 5«CCyO-, 55 87-C-45O7-0, 38 27-3. 3i 33-0, 2 625-J. 21 34-0^ 1703-0ci327-Ocl0O4-0.0732-0.05C9-C.c33I-0o0i94- j.0J95-o.0C3l 0. 

1 5Ct 00-0 5321-C 43C7-C-. 3546-0. 2 93 6-0. 2423-J. i 977-J. 1 58 i-0 ,1 232-0 ,,0930-C*. C676-C. 04 69- C. C3C4-0. Jl 78- j. Ou87-o* 0029 0. 

1 6 5o CO-0-, 51 6 9-*.. 41 63-0. 3 46 J - 0, 2 963-Oc 2 6J7-0v 2324 203 2-0 ai 7U-J .1 367-0. 10 28-0, 0724-C. 04 74-0. 0276- j. 0137-0. U045 0. 

176. 25-1- 5{ 2 3-0. 4126- 347C-0: 3 1 09-u, 3 fJ95-0* 356 2-.)- 3902-0. 4U72-3 » 3799-0. 3024-Cc 2C86-0 . i 36 2- 0. 0805- U. 0397-0. 0131 O. 
176,50-C - S' 22-0. 4; 76-C,94-»3-:,3i i.O-0,3099-‘Jft3570-J, 39i4-C,,4t84-0. 3815-0. 306fi-C,2le5-Cal425-C. 0852-0. 0423-0.0140 o. 

1 76, 75-Oo 5X 22-0, 4i2F-0. 3470-0, 31U-0, 3103“U« 3577-3, 3924-0 ,4095-J . 3829-0, 3*05-0. 2235-0 . 149 4- C . J9U3-0. 0451-Oc Jl49 u. 

177,. nc-C.. 5'i 2!- : 4J 2S-A, 347 j-'Ja 311 2-0. 3X0 5-0^35 84-J . 3933-0 ,4 1 l 3-3 . 3844-0. 3 142-C. 2 301- 1 . 1 563- 0.09 60-0. 0482-0. 01 60 u. 

;77« ?5-0>5l21-Cc 4l2f-0-.. 347J-^.3U3-0, 3U 0-3. 35 89-3, 3941-0 ,411 1-0 ,3856-0. 31 75-C, 2 365- 1 * 16 50-U . 102 3-0. 051 7-0. Ui72 u. 

50-0. S: 2 I-Oc 412 S-n„347i;-0, 3114-0, 51l3-Oc 3594-3: 3 947-9 = 41 16-3. 3865-0. 3206-0. 2430-0. 1 743-C.iO94-0. G558-U. 36 o. 

17 7- 7 5-n.. S: 21-0.4 - 24-i,-,3 47.J-0. 3U5-0: 3it 6-0« 359 8-J, 3953-3'>4 121-0 . 307 3-C. 32 32-0. 24 86-C . 1 3 2 1- O.i i 76- 0. 06u5-0. 02o3 U. 

J 7 8, 00- 3 5T 2C-Ct 412 4-0.34 70-3: 3* i. 6-0, 3U 8-U; 36.) 2-J 3957-0.,4 124-0.3379-0. 32 54-C. 253 5- C cl9C7-u.i2 72-0, 066 0-0, J224 o. 

178r 25-0: 51 26-0.412^-0. 3471-0: 3:17-0. 3120-u. 36^^5-J. 396i-C< 4J 27-0.3883-Oc 3<71-0. 2560-C p199l- 0 .*387-0. 0726-0.0251 u. 

1 78-- 50-Oc 5> 2C-C. 4?24-0,347X-a. 3U7-0. 3122-J; 36J7-3-. 3963-0 .:4 128-3 , 3886-uc 32 85-C. 261 5-C . 206 i- 1 « t438-0, 030 8-0. 0287 0. 

)78c 75-0 s;r 20-C.4}.24-Oc 347V-fj; 71) 8-0, 3) 24-0, 36ij8-3 , 3965-0 ,4129-0 *3888-0. 329 5-0. 264J-C .21 i 5- G . 159v>- 0. C9U-U. 054* u. 

17<;, GO-C: SI 2»:-0' 4,j 24-0o 3471 -;)l 31} 8-0. 312 S-0., 3609-1.. 3966-0 -4 129-0 . 3 9S9-0. 3302-0, 2 6 63-C . 2l 5 5-0. 1 67 l-O, 1049-Oo0425 0. 

:7<5 25-0 512C-Ce 4124-0.. 347J.-i}> 3U 8-j, 312 6-0. 36 v 0 -J ^ 3967- Jo4123-0 . 3838-C. 3306-C. 2677-C .21 8 :-0« *726- 0. ii 8i-0. 0573 0. 

179. 50-0. 5: ;s- Co4J 24-y , 347;.-.j, 3U 9-0, .1). 27-Ut 36 1 u -J 3967-0.4127-0. 3 887-t. 33u8-0. 26 84-0. 2*97- 0 .i 750-0. *-253-0. 0673 u. 

179, 75-L-. S7 iS-i, -412<‘-U: 3471-0, 311 °-L/, 312 7-0. 36* O -3.- 3967-0 ,4; 27-0 . 3886-3. 3308-0, 2687-C.22C 3-0.* 7 69-0. 1 279-u. 0693 o. 
.I8O0OO-O S) 7.9-0., 4*24- Cc 3471-.;, 31'. 9-J, 3127-0^ 36-.0-3* 3967-0.,4 126-0 . 3886-u. 3308-0, 2687-0.2202-0.1767-0. 1273-0. 0oo7 u. 


PASS 4 cdllTwS- 
STPP ^ 

S 179:6A .1: nro -2. 822F-U2 

3.- 5lJ9*04 

Ui 

.U062 6+31 

1.387E+03 

2.2C1E+I5 

1.867E+17 

3, 8716-09 

9TFP h 

S 179:52 

J 970 -lc27jF-0.l 

2. 013F+C4 

•Ja 

9, 139E+0C 

3.983E+U2 

2.454E+15 

6.182E+1* 

3.8106-lu 

step '.7 

5 179 , 40 

G>96) -U7*6E-Cl 

Zr. C8'. F+04 

Ja 

3.347 E+jO 

3. 6G4E+02 

2.65SE+15 

2.543E+12 

1,1646-10 

STEP 7P 

5 179,28 

0, C44 -2- 998E-0;. 

i, 36.’.E + 04 


7.774E+00 

2» 593E+02 

2.8£CE+15 

2.436E+12 

5, 9566-11 

STPP 24 

S lT9^i6 

0^9 29 - 2- 24SE-01 

1.1 999*04 

C'3 

7c264e+Jt 

2. 216E+02 

3.069E+15 

E.255E+IJ 

3= 8456-11 

s r FP 8(J 

5 179:(J4 

3^^S*.4 -2c 456E-C2 

lo 2148 + 04 

Or- 

6 0 8 1 9 E +0 0 

2c 04 1 E +02 

3.28*E+i5 

2.771E+11 

2.8*96-11 

STFP 36 

5 ^8,92 

3-9C0 -2 6S8E-CI 

28SF + 04 

Jc 

6. 4;>8 6+00 

1.928E+02 

3.49E6+i5 

1.1566+11 

2v236E-U 

STEP 42 

5 '.78 80 

'i,^oq7 - 2; 859E-01 

Iq 2 989 + 04 

Jo 

6*020E+j0 

1.782E+02 

3,7166+15 

5,6176+10 

1. 8706-11 

STFP 48 

S 178 6.8 

0 ,875 - 2' CSIE-C* 

1: 313F+04 

Oc 

5.660E+00 

1.654E+02 

3.941E+15 

3.C56E+1C 

ic 624S-U 

STEP 54 

S 178^56 

C 8£.8 234F-G]! 

1, 330E+04 

U;. 

5, 326E+J0 

1. 54uE+02 

4c 16 56 + 1 5 

x.ei36+lC 

Iv 450E-11 

STEP 69 

5 *78:44 

C.:R52 -3c 3 80E-0I 

J , 03 6F + 04 

Oa 

5a042 6+00 

ia 251E+02 

4.39*6+15 

1.1516+lC 

1.3226-11 

ST PD 66 

5 178„32 

OcE42 -3c 5C9E-CI 

C45F + 04 

0:. 

4a78fl 6+00 

i. 1736+02 

4.619E+15 

■(.69EE + C9 

*.2246-11 

9 T FP 7 > 

5 I78c?0 

0 832 -3c £33P-0! 

0S4F+04 

0. 

4 j 550 6 +JU 

la lOlE+02 

4.846E+1S 

5o375F*CC 

l«.i47E-U 

STFP -^8 

S '78, 'iS 

CaE?3 -3: 752E-01 

065E + C4 

J= 

4. 326 6+.>0 

1:0376+02 

5.C8CE+15 

2.851E+C5 

1.C86E-U 

STFP R4 

S 177 06 

C 8\4 -3;. R67F-01 

*,'376F + 04 

0 

4^iy.5E+3C 

9. 792E+01 

5.2126+15 

2. 904E+C5 

1. 0376-11 

9T^P O') 

5 J’7,, 8 4 

C:8C6 -3„ op F-01 

3. 703F + 03 

0. 

3.079 6 + jG 

5. 838E+ui 

5.5546+15 

2.222E+CS 

9, 9556-12 

step P6 

S 277,72 

C- -797 - 3c 9466-01 

72 3F + 03 

Ja 

3c 860E+1C 

5o 576E+01 

5, 7956+15 

1.732E+CS 

9,5976-12 

STEP JO,' 

S 177- 6 ) 

3 :790 -3. 9RUE-01 

3: 746F + 03 

J: 

3.7466+uO 

5a 332E+UI 

6=0416+15 

*,374E+tS 

9.2806-12 

STFP 108 
5 7 77c 48 

:-.,7 82 -4.. .j:. 3E-9*. 

3, 771 F + OJ 

Oj 

3. 637 E +00 

5, iC2E+0Z 

6.2936+15 

1,1056+09 

Sc 9996-12 

STFP i 1 4 

S 277.36 

rj.,774 -4. 045F-J1 

3: 700E+O3 


3. 532 E+CC 

4, fl88E+J^ 

6,5526+15 

9.u’i4E+c6 

8. 7496-12 

STFP *23 
S 5 77. 2 4 

Cm, 7 67 -4. C77F-01 

3: 827F+03 

J; 

3=412 E+jO 

4.6R6E+U1 

6a816E+l5 

7.444E+C8 

8. 5276-12 

s T FP ’.26 
S ?77oi2 

0.7 60 -4.,i'.,7F-0\ 

-2, 267E+03 

U. 

3.317 9+.70 

2. 13JE+01 

7.085E+15 

6.21*6+08 

8.3246-12 

STEP i->2 
S :77:,00 

0.753 -4 G83P-0‘, 

-2:2719 + 01 


3o293e+:)w 

2. llJE+Oi 

7.3626+15 

5c22EE+oe 

a, 1366-12 

STEP ■‘3 8 
5 ’ 7.6- FP 

0c746 -4;06i'.F-01 

-2, 2 73F + 03 

L a 

3s24OF.^0 

2o O90 6+uj. 

7.652E+15 

4,43iE+tE 

7=9556-12 

STFP ;.44 

5 *76 . 7ft 

0..,73O -4 C36F-0? 

-2i 27SF + 03 

Oa 

3.*20SE+uO 

2s C696+0- 

7=9546+15 

2.781E+08 

7.781E-12 

STEP i SO 
S .’ 76: 44 

0 ,733 -4: U2 3E-3'. 

-^c27SF + 03 

0. 

3o tS2E+:i; 

2e u49E + i:;i 

0=2686+15 

3-a249E +U6 

7=6186-12 

STFP IS 6 
S ’76. S'* 

0 756 930 F_r: 

-2: 2 7S9 + 03 


3r. 1 19 E+Ou 

2ou27F+ul 

Ba593E+*5 

2.8036+08 

7c 4546-12 

STFP If.? 

5 176 40 

t .7Z»: -3' 9659-01 

-2: 2 73E + C3 

0= 

3i J■»7E♦J^. 

2,0C6E+-jX 

9.932E+J5 

2.4326+Cf 

7, 2986-12 

STFP *a 

S *76.28 

C .7 , 3 - 3; c4,= -0'. 

-2-, 2 70E + O3 

J: 

3,.;3Sr+;,C' 

ia 985E + 01 

9.2826+15 

2.122E+C8 

7c 15C6-12 

STFP 176 
5 ?.7A., 16 

2 ,7ft7 - 3. 4 ?OF-y'. 

-S- 03qr + ^,> 

-5. 889F+0?. 

2 ■ 879E+UC 

Ps 532E+^a‘l 

1,282E«)6 

1.6676+08 

7.02*6-12 

STFP ;.<?j 

5 ns, '4 

C 70’ -3, 407E-0-. 

-S,4'ISE*02 

-5a 743E + G3 

2a697c*j^^ 

7, 9426+04 

l:462E+i6 

Ic6546+L6 

6,9ilE-12 

STFP '.86 
S .75,02 

0,696 -r, 383E-3,’ 

-2: 019" + 03 

-5.- 5n9E+u3 

2 a s 48 E + J 

6= 4996+o; 

io645E+j6 

i.48lE+08 

6.8J3E-12 

STEP 107 

5 -.75^80 

0 6 9C -3. 355F-i;‘ 

-,' v o5 8F*o:> 

-5^43 3E+..3 

2 - 421 F+.',*, 

5c 6i lE+Ol 

i. 8286+16 

la33eE+Ce 

6,7746-12 

step iqo 
S 75.68 

0_s6 8£ -3: 3 29F-.:\ 

-898E*03 

-5. 2946+33 

2.3U 6+jO 

4.9i*E + 'J* 

2.014E+ J6 

lc2i4E+oe 

6. 7346-12 

step ?w4 
.5 ,V 75 . 9 6 

C 6 a- 20 SF-OI 

-1-, 837F+03 

-5: i97E + 03 

2:2'5E + ;o 

4, 3476 + Oi 

2.2C.E+16 

i« lOdE+O 8 

6.6 966-12 

STFP 7 ?;• 

F T7«1.,4 4 

G:676, -3. 2P E-c . 

776F*ui 

-5. J25E+03 

2. 129 E + >t; 

3. 8866+jl 

2.3896+16 

1.0156+08 

6.6746-12 


i 


9TCO 

6 




5 

1 75 3*? 

>:■ ,fTZ 

-7. 258C-C. 

-• , i75P*-03 -4. 893E + 03 

?.T fP 

?2» 




5 

1 75- 20 

C- 668 

-3< 2 37F-3. 

-I, 655F*03 -4.776E«-07 

5T PO 

?2 8 




5 

? 75, j 8 

C 654 

-3: 2i6£-0i 

-:, 594E*j3 -4,657E**;3 

5TEP 

5-S9 





1 74 .Of 

0.,66C« 

-3. 1 96C-0’ 

-\.531 = *03 -4,542E«-'j3 

9TPP 

340 




5 

74-. 8 4 

0 j6 ^6 

-3 178F-0: 

-1«47.?P*03 -4,431F«-0? 

ST‘^P 

244 





‘ 74 7? 

C..657 

-- 1598 -C-. 

-i,4!2P*03 -4,323E*03 

CT CD 

752 




5 

1 74. 60 

0-64.'i 

-7r 147=-C*. 

-\:35;.E + 03 -4. 2i8F+D3 

5T CO 

755 




5 

174,48 

9 rf 4f 

-3.. 12 5P-0; 

-1.2 9JE + 03 -4. ,15E*-J3 

9TPP 

^64 




c 

174. 36 

0,64:? 

-3. iJOE-C' 

-1.23 3F + 03 -4;015E+0 3 

STCP 

7T0 




5 

174^.? 4 

0 638 

-3- 0948-91 

-1..J69C + U.3 -3.918EK-3 

c r CD 

77^ 




5 

174^ \ 7 

C ,6 35 

- 3, 079t-0i 

-1..1C3E+G3 -3-. 82ZFV03 


2,J5I 

3. 5u7E+ol 

2,579E4i6 

9.349E4C7 

6.661E-12 

i,98\ E+-J0 

3, '.79E+0* 

Z,77(.E4;e 

E.629E4C7 

6.646E-12 

l,9i3E*iO 

2.9i. ^z*j t 

2,962E4i6 

E.C:9E4C7 

6.64Se-i 2 

i,86GE+^G 

2, 66SE+01 

7,i56E4l6 

7.441E407 

9,640E-12 

1 9 8 .;6 E+C'O 

Z-, 46?,E*uI 

3.351E416 

6.947E4C7 

6.647E-12 

lr756E+>'C 

2, 28J.E4J1 

3,547E4!6 

6.49CE4G7 

6.653E-12 

1,710E*1C- 

2, l23E*ul 

3,744E4ie 

6,-w85E4C7 

6. 659E-12 

1.667E400 

:,993E401 

2.943E416 

5.72AE40 1 

6.6756-12 

627E+00 

1, 958E4C1 

4,x42E416 

5. 29EE4C7 

6,6946-12 

I , 589E4u-0 

lo 746E+01 

4.343E4U 

5.^.81E♦07 

6.7026-12 

1.554E+:-0 

1.646E40:-. 

4o544E4i6 

4,8CEE+l7 

6.723E-1 2 


PHS VALUe^r 

Tmi=179 6A? THP = 

AVIAL 

pdtential 

OISTP tPlUT inN 


RMH 

-V 

0.^633 

0 ’51 

0. •'OG 

0. 397 

V r 76C 

0 ■ 4 * 3 

t? . 8? 3 

0 389 

0, 8*C 

G 33 < 

0.9CC 

0 ?59 

0.933 

0 270 

O.. 96G 

0 177 

0^98C 

0. 177 

0 993 

0 067 

i,OCC 

-l>7 


rNITIAl VA( U?<- 


I 

OEt. I 

0611 

THE TA 

VSM 

VI 

INJ EC. 

PSIC 


( R^^)\ 

(n£9j 

<PEG) 



angles 


7 

0;. 2 

•% ‘)90 

\ 79, 9LC 

0 93j 

7c 20 E+^- 3 

2. OCO 

C, 

1 

r, ,r;: j 

0/^ 79 

J 79, 821 

J,8J3 

5,4U6+U3 

3v CiOC 

C j 

4 

0.0C5 

0 369 

179,731 

Jc 7-J'J 

5 6C’E*0 3 

4. Coo 

c. 

5 

0,-C G6 

3 5 8 

170, 442 


4, 80 E *0 3 

5, Coo 

w , 

6 

0,C08 

0.44 8 

170,852 

U.5O0 

4..0OE+0 3 

6o OOu 

c« 

7 

0,CC9 

Gi.537 

7 9,463 

J»4j) 

3.. 20 E ♦0 3 

7, OoO 

c. 

8 

OtCil 

Or 627 

179.373 

J,2D0 

1 .> 60 6*5 3 

8. OCO 

w. 


**NP 6AP** 


PCTENTfAL FietD 


INPUT- 






RHC 







THETA 

0>5 6C 

vi. 633 

Cv 7flC 

0.75 J 

0 0 8 1 3 

Oo 960 

0.900 

0.933 

C.96C 

C, 980 

0.993 

l.uoo 

4 5- nc 

-1 OOE*0C‘ 

-Ic OCE+OO 

-1, OGF+Og 

-1,UOE+0- 

-1-.0 jE+CO 

-ieOOE+00 

-i» uO E+uO 

-..COE+Cu 

-i.oc6*ocj 

-X, 006 + 00 

OOE+00 

-..wOE+GG 

50 no 

0 ; 

Co 

0.) 

C. 

0 = 

0. 

Gm 

0. 

0. 

C« 

0. 

-8.jw 6-0 1 

7 5. CO 

0. 

0. 

0 

Cro 

Or 

Go 

0. 

U. 

c. 

0. 

0. 

-6.WV E-vi 

90c 00 

0 

0« 

Oi, 

Oj ■ 

0, 

0 

C 3 

0 . 

Cj 

0. 

0* 

-4.00 E-y 1 

:o5.. cc 

0 • 

“•7 

C; 

U> 

0- 

0, 

Oj 

C. 

c. 

0. 

-J. 

-^.uo E-gi 

120 cc 

C) . 

3 

Oc 

G: 

0: 

0 0 

0. 

c. 

c. 

0. 

0. 

-0. 

13 5i CO 

0 

0,. 

Oo 

G. 

0; 

y 

03 

0. 

l. 

c. 

0. 

0. 

1 50=. CO 

C*; 

P'C 

09 

0. 

0 

0 c 

0. 

u. 

C.9 

u. 

0. 

y. 

1 6 5,cc 

0. 

'Jk 

0; 

0. 


Oo 

0, 

c , 

c« 

0. 

0. 

0. 

}7fr 25 

0. 

2v. C4C + 01 

2o 85F+Oi 

4. 49E+01 

7o '} 3E+01 

2,4iE+Cl 

0^ 

0. 

c. 

0. 

0. 

J. 

176, 5C 

«v. 

i 3 96E + 0\ 

2o 72E + 01 

4v ilE + Oi 

5.81E+01 

5. 24E+01 

0. 

0. 

0. 

u. 

0. 

u. 

176c 75 

0 , 

' 9 886 + 01 

2. 58E+01 

3. 856+Oa 

5 346+01 

lo03E+o2 


c. 

Cw 

0. 

0. 

u. 

177.; CO 

0, 

lo 79E + 01 

2: 45P + CI 

3o 6i6 + 0:. 

5.68E+0i 

8„85E+0l 

u. 

0 • 

03 

0, 

u. 

0. 

177„ 25 

0 

1. 71F + 01 

2, 33F + a 

3, 39c+0i 

5.n .i 7 E+01 

7.536+01 

l.P3E*C2 

0, 

c . 

u. 

0. 

u. 

177, 50 

J.> 

?.r 6 2E + OJ. 

.2i 2 0E + CI 

3, i9E+Ui 

4.826+01 

9. J5E+0i 

X.4BE+02 

y • 

Gv 

0. 

0. 

u. 

177^ 35 

n. ' 

i.« 51F + 01 

2. 07C+CI 

2. 99E + 01 

4 ,48E+al 

7„ 26E+01 

x,20E*u2 

0 . 

0. 

J. 

u« 

0. 

•. 7flr CC 

0 

?.r, 39E + 01 

1 92E + OX 

2, 79E+U1 

4,16E+01 

6.656+01 

l.,u4E + G2 

2,34E*C2 

0. 

0. 

0. 

0. 

179: 25 

0- 

Jr 279. *0\ 

:. , 75F + C. 

2,59E + ol 

3-. 86E+01 

5.i2E+0i 

l.o96+u2 

2,C?E*C2 

t. 

-J. 

0. 

u. 

) 78« 50 

07 

lo J4£*.at 

1. 5 8E + Oi 

2: 33E+01 

3 : 5 4E+0x 

5.626+0/ 

9..946+Oi 

2.f3E*C2 

t. 

u. 

(/• 

0. 

:78 35 

0 

Or 89E + G0 

1. 3RE+0* 

2. 05F+0i 

3. 16E+01 

5 u lu F +0 1 

8.396+Oi 

1.16E«C2 

’3 1*6*02 

u. 

Jo 

u. 

J 79. CC 

0 , 

8<. C9E + 00 

1 l6F + 0t 

r> 75E+CI 

2.71E+U1 

4t 49E+oi 

7.946+01 

1. S5F+L2 

3.446+02 

7. 51E + CZ 

Oo 

u. 

17°; 25 

J 

6- 15C + 00 

8. 93E + O0 

1 ; 38 E +0 A 

2.,21E+uI 

3. 7JE+0i 

6.78E+01 

i,34E+l2 

2. 906 *U2 

7, G2E + 02 

u. 

0. 

J 79 50 

0. 

5c lOP + OJ 

5, 99F + OU 

9, 38E+U0 

ir5 5E+01 

2« 746+01 

5. 16E+yl 

l.cTE+02 

i. 396+02 

5. 56E + 02 

1.29E+03 

u. 

179« 75 

3. 

2.:C5F+CG 

2r 99F + C0 

4,69E+00 

7o77E+o*j 

1 40 E + 0 i 

2,92E*Ui 

6.396*01 

i. 556+02 

3. 92E + y2 

9.3iE + u2 

y. 

lft.:«.CC 0- 

7 1 =- 3 ,no 

Oc 

0-. 

U. 

0. 

U. 

0. 

0. 

Cc 

0. 

o. 

W. 


pctfntial PIEUD 

OUTPUT- RHC 

THCTA Oc^OC C'<3r.O 0.393 0^480 U. 5bO 0,633 0,7JU .1, 760 O.Si3 0,660 0.900 C.932 C.96C 0.980 C.993 i.OUJ 

*^5oCC-l- COOC- '.oCUDO-l-. >70C:>“1. 000 J-l. J00ij-i.i 00 JU-i., OOOU-l.OUuO-i ,OuOO-i. JtOO-i,CCCO-l.{JCCC-l,JO-jO-l.OOOO-i,QUOU-l.L»uUu 
60. 00-0,. 885 A-Ct8603-Uj 84^1-0.. 82 8<‘-07 8l92-:-. 8i J8-0^ 8055-0 , 80 20-0 ^7999-0, 79B8-0.. 7 9 S4-L, 79 E 7969-0, 7994-J. 7998-0. 8uUU 
75r 00-0 . 7919-(.\ 7456-Ui 71i2-0. 6845-0.. 663 6-0. 647 2-0- 6345-G . 6248-0 . 6i.7 4-u. 6U8-0. 6C77-L.6C6e-0,6t27-0, 6Ji 3-0. 60o4-o. 6uwJ 
9Gc CC-C-. 7J.4 C--:s 6493-0, 5998-U; 5599-0, 527 i- J., 4999-Jv 4774-C .45 8 3-0 -.44 37-0, 6 31 6-0^ 42 16-C, 4l 39- C ,6u82-y. 4 J4U-U, 4ul 3-0. 4uUu 
>0 5^000.- 64 9 P- 0 .. 56 oft-o 5067-0. 45^^3-o» 4094-^, 37'J 3 -7 . 3 362-0 3066-0 o 2Bi 1-0, 2596- 1. 281 7- C.22 73- 0 , 2i ti-0. 207 9-0, 2 026-0. 20uo 
?20o00-C.o59e6-0^5t;6fl-0c 6342-0. 3 725-0. 3 182-0., 2693-j 2247-0 o.: 8 39-0 ..146 5-U , li26-tj„ C 624-0.056 1- C.O'343-0. Gi 74-U, 0059-w, 

*. 35 C >0- 560I-«>; 46>. 4-0. 3 846-0. 3202-0. ^ 643-u. 2i 5 I -0, 1717-0 : i 3 38-J 10 1 3-..X. 07 39-C, C 5i4- C.C3 34-C . Jl 96-0. 0U96-J. U032 O, 

1 *tO^ OC-0,. 53 3 8-0, 4379-0., 35 75- 0, 2 969-yj,. 245 9- j : Zt; I 3 -:j . 1 61 5-0 , i 26 i-O ,0954-u, C698-0. C8 61- G.C3 12- v ,Ji 83-0, o090-G. Ow29 0. 

T6 5,jC 0-0 518 7-0, 41Q1-0- ?50.l- 0:. 3 Ol 8-0. 2 67 9-Oc 24 1 3 -0^ 213i-0 . 1 81 J-0 ,1455-0. 1093-0. 0 769- t.C 5C 2- C .029 5-0. 0i45-U. 0 l48 u. 

1.7 6,2 5-0. 514 2-17, 4?. 5 6- O, 3 5, '.7- Or 3J 3 2i i -j. 3750-J , 4191-0,.4479-0 . 4 3u 6-0 . 32 5 i-C . 22 20- 0 . 1442- C, 084 9-0, 04l 8-0. Oi.38 0. 

176,- 50-0, .514 1- u- 4156- 0.3 51 7- Cl 31 83-0. 3 21 5-G, 3730 -0 . 4204-0., 449 3-0 . 43i 8-0. 3309-0. 2296- C . 1 5C4-G ,:j893-0, 0444-u. 0146 w. 

176« 75-0. 514 J-Or 4? 56-0 354 7-0„ 31 84-0, 3 2,i 9-0, 376 8-0 ; 4215-0, 4504-U , 4.328-0. 3367-G. 237 7-0. 1 5 7i- G ,j94l-y, 0471-0. Ul 56 O. 

'.7 7 CC-Cv5'i 4 3-C.4; 56-0 3 51 7-C- 3 «. 85-0. 3 22 3-0. 377 5-0 . 4226-0 ,4513-0,4337-0. 341 6-C. 2466- C . 1 644-0 ,0993-0. 0502-0. OX 67 u. 

17T, 75-G 514C-C 41^5-0,3517-0. 3i 86-0.3227-Uc 3731-0; 4234-i'. 4521-0, 4344-0. 3459-0, 2564-C.1724-C,xo5G-0,g538-0,0*79 0, 


25 



1 tl. 3«L?-C. 3 ■ *? ‘^-Cv 2 J J-j.. 3736-.) , 4^4 2-C ..•♦5 27-J 4349-jt 3497-j . 2 64 2-uc i 3 1 3-u « i* i 4- J, u579-g, ui 9<» o. 

' 7 7., 7 5-0 9' 4^-.' 4:..S6-.., 3 317-C-. 3i 83-C o 3 233- j. 377 l-j 4243-i.- .-45 3 2-J i 455 i-J-. 2520-C > 2 7v.7- C - * 5 1 2 - 1 .* i 8 7-U. 062S-w. i u. 
*» 7 3<r CO-0 .5* 35-0. 41 5 5-. r, 5 53,7- C. 3 . 37-0.- 3236-^- . 3775-:';^ 42 53-0 ^45 ?5 -m .. 4 35 5-o , 35 54-C > 2 16 2- C« 2 ^ 26-v. . i27l- u, u585-u^ j233 w. 
.78c? 5-0 5i35-.. 415?-.,^ 3513-0 31 9 j-y. 3 .23 8-j. 3773-j, 42 57-0 45 33-U-.4356-. 3575-C. 28o8-C* 2 a i9-o uA365-'j. w765-l», u2aw o. 

’73,50-0 5-’*39-»; 41.54-i .35! 3-,}, 3;. 9».-C-,324i-.).. 3300-0. 426^-0 45 37-'j ;4356-v. 3592-C- 2 8-<6-C,2l55-w, *478-0. 0363-O. u274 o. 

i7?.. 15-0 5‘ 3 5-r,4i54-c.35lH-0.. 3i9:-u,3 242-x.. 33 . 2-^.* ; 42 5.?-0..4S43-j ^ 435 5-f. 3604-0., 2 67 5-C. 2 2 5 3-y , a 61 5-0. u994-o. U34w «• 
' 79, uC-0.. 5t 3 5-, 4i 54-r , 35- 3-0 3 ( 97-Ov 3 244-0^ 33j3-J . 4263-j ,45 39-0 ..4 354-i,.'. 3fcl5-0- 2 856 -Ca 22 56-0 <. i 7i 7-'J. 1^75-0. o4j3 o. 
175r 25-«> 5V35-‘}r4154-0. 35?. 3-^, 3 a 9.'»-0, • 245-0 . 33J4-J, 4264-0 45 39-0 ,4.352-0., 36i8-w. 251o-0-. 2324-C. a 786-0. 13ww-0* w49v u. 
179* 50-0 5>.3«-C«4l54-r, 35i3-0- 319t»-,j. J246- j„ 33g 5-j , 4264- s/ . 45 33-0, 4349-0 . 36 2u-,. . 2 5 1 7-,, * 2 34 1 - 0.18 2 6-0.13 66-0-, J644 u« 
■? ?e.> I 5-0 5? 39-t:- 4’ 54-.'-. 35). 3-y. 3i93-t, 3246-Oi 38y5-J- 4265-0 4537-0 .4 343-i, . 2621-u.. 252J-C.,234?-y.i84i-0*l 338—0,0688 U. 

i30,C0-O. 51 .3 5- Co 41 5 4-l, 351 3-0. 3i 93-0, 3246- 0 , 33 j 5-j , 4264-'j ,45 37 -u . 4347-0. 362i-C . 2 52o-t v2346-C . l340-0. 1379-0. 0565 ,/• 

TABLE OP RHC-VALUES FC3 
FOUI^CTENTl AL LINES 


theta 

-C^o^ -Oc 30 

1 

0 

-j 

c 


-J-60 

-.is 50 

-0.40 

-U.3C 

— C • 2y — y. '* ‘J 

45., CC 

60.. CO 
75 CC 
90.-, CO 

1 0.5c 00 
■ 20. CO 

r Oi 

C',175 Ov 8;. J 

C,.l96 0.192 

C C7u Jcl40 

y. V. 57 t).U4 

0 0 5u Oc.*.->u 

Or, 

J 

C.43‘j 
2 22 
i 7 :. 

0 A 49 


Uc 39 3 
y 262 
u.a99 

.ic 
0 , 

•J 6 33 
y.,4t)4 
J-309 

3.-. 

0. 

1 -OjO 
0.578 
0 ,441 

0.1 

0. 

U. 

0.774 

y.537 

L c -J. 

t, 0. 

c. u. 

C. J, 

i.yoL 0. 

C.736 0, 877 

135.CC 

0f.C45 0. 091 

0,136 


0.182 

J 26a 

0,375 

U.5C9 

C-657 0* 

846 

: 5u-. oc 

C->C43 j.-C86 

0: t29 


0.a?2 

Ja>33 

J.341 

0.476 

0^,63 5 0. 

80 5 

165. CO 

Cr,04?. 0«CR3 

0:125 


j'. * 66 

0,219 

0.326 

0,484 

0.724 iJ. 

87i 

176 25 

<:,,:.41 0.082 

Ur l23 


y * 465 

0v214 

0.827 

y. 870 

y.9o9 0. 

953 

? 76o 50 

C.C41 0. 032 

0- i 2 3 


t-4, 165 

U.,214 

. 0,828 

0.872 

C . 9 i 2 J. 

955 

176., 75 

C;C4r o*o«.? 

0: :23 


Gcl65 

(i..2i^ 

0.829 

0.875 

0,916 Oc 

957 

•' 77 00 

C-t4\ O.032 

U. 123 


0.16 5 

0.214 

0,8 3j 

y«878 

C.9i9 Jv 

960 

177o 25 

0 ,(>41 Or C3.*: 

0: 123 


C. 1 65 

0^214 

y.83i 

u,eei 

L,922 0, 962 

’ 77- 50 

u..i:4l 0-. 082 

ti, i.23 


yvl65 

0:214 

y ,8 32 

U.683 

L.926 Oc 

964 

1 77 7 5 

C.C41 0. C32 

y- i23 


0.165 

214 

0*833 

y « 6 8 6 

tt:930 0, 

967 

1 78c CC 

Cc i;4i Oo 032 

0 123 


y. A 65 

0.-.214 

0.834 

0.888 

0.934 J.969 

's 73 .. 25 

(J,04:i J. 032 

Ut 123 


L. c 6 5 

0.214 

0,835 

0.S9C 

L.9 38 Jo 972 

' 78,^ 50 

0 04a 0-fi32 

U: 123 


y- 16 5 

Ur214 

0,835 

0.892 

y a 9 4 1 y. 

976 

7 fir. 75 

C.-C41 O.U32 

-1. A 23 


y. 165 

0^214 

0,835 

0.893 

C.944 y« 

98J 

l'^9 OC 

C : 04: 0.082 

0. 4 23 


U« lb5 

y..214 

0o836 

0.694 

y.947 y. 

983 

179 25 

0 C4? Or. 082 

123 


y. 465 

J.214 

0.836 

U.89S 

y.949 0. 

985 

T"^9., 50 

0- r4x U.082 

0 *23 


Of 16 5 

y.,214 

0.836 

0.895 

0.951 J. 

987 

:• 79 . 15 

-■,04i. OeOB2 

Oc J 23 


U« i. 6 5 

0o214 

0-,8.36 

0.895 

C.952 0.967 

180 CC 

0,04:. 0.032 

0, 123 


Uo 165 

0.214 

0.836 

Oc895 

0.95Z 0.987 

STEP 0 

2 ;• 19„ 97 

i CCO -3c 808E-02 

3,597r4-04 

0, 


2.715E+01 

5. 41uE+y 3 

5ilieE+14 

1.8816+17 

2.4286-10 

3 179o82 

loCCO -3- 933E-Q2 

3c55.3E + v4 

Go 


1.789E+UI 

2- 883E+U3 

1.0386+15 

1.8786+47 

9.707E-IC 

4 179.1 73 

1,CC0 -4.027E-O2 

3, 516E+04 

0. 


1. 3346+01 

4. 885E+03 

1.617E415 

I. 8746+17 

2.1826-C9 

5 119.64 

1, CCO - 3 .. 965F-02 

3.504E + O4 

Oj 


Io0 62 6+Ul 

1.41OE+0 3 

2.4576+15 

1.8676+17 

3. 8746-09 

6 179.55 

l.OCC -3.902F-02 

3,491E«-04 

0- 


8.821 E+OO 

I. 1686+03 

2.563E+15 

4.8556+17 

6.C29E-G9 

7 1 79.. 46 

I.OOU -3o 696F-0.» 

3^509E + 04 

j. 


7.538E+00 

U060E+03 

2.751E+15 

4. 8316 + n 

8.6256-09 

8 l"»9. 37 

IcCOC -3c 283E-02 

3.569E+04 

Uo 


6c576E+0G 

1. 633E+03 

1.6556*45 

i.77CE + n 

4.1546-08 

STPP 6 

2 179.79 

C .-.050 -2.033E-01 

1.516E + 04 

y. 


2.172E+01 

1. 6376+03 

6.4596*14 

7.1156+12 

7.3846-12 

3 I79i70 

Ur 966 698F-01 

i , 81 BF+04 

0. 


ic494E+Ji 

8. 678E +0 2 

1-232E*15 

6.1216+43 

4. 5646-11 

4 179.^61 

0 9 75 -ic 500E-o:i. 

lo 82 8F + 04 

Ou 


ic i24E+Jl 

5. 257E+02 

1.8946*15 

2.5056+14 

1.5866-10 

5 379,52 

G.9 80 -it 377E-Oi 

1. 837E+U4 

0. 


8.975 E+00 

3. 70UE+02 

2-4866*15 

6.423E+14 

3.8856-lU 

6 179.43 

0,9E3 -\ IT9F-01 

4,A74F+04 

Uo 


7.458 E +UG 

5. 313E+U2 

2.959E*i5 

4.2136+15 

7.5766-10 

7 170,34 

0..986 -.XrCO3E-0r 

4d 664P +04 

y. 


6.321 E+uO 

4. 836E +o2 

3.2106*15 

1.8296+15 

1.2546-09 

0 179., 25 

C.988 -6.- 127E-02 

4., 853E + U4 

0. 


4.093 E+00 

6. 908E+02 

2.0886*45 

2.3416+15 

1. 8^06-09 

STEP 12 

2 179o67 

CoOCS -2 c 787e-01. 

1.378F+04 

u. 


1-.B73E+D1 

L. 261E+03 

7.7476*14 

2.G2CF+11 

2.8IU6-12 

3 1 79,. 5 8 

O. 5 O 37 -2t28.3E-0i 

i,535E+04 

0» 


1.337 E+Oi 

6.835E+02 

1.3966+15 

4.7896+12 

1,4406-11 

4 179,49 

C ,«52 977F-0i 

:.547E+04 

G« 


1.024 E+Oi 

4. 245E+02 

2.0836*15 

e. 5 5 7E + 1 2 

4. 793E-ii 

5 179.40 

.1 ,962 -Ic 773E-0*. 

lu 08UE+O4 

u. 


8.241E+GC 

3. 337E+02 

2.6916+15 

2.707E-H 2 

1.2U36-10 

6 1 79o 3 \ 

0.968 -Ic 6 COE-O 1 

i , 918E+04 

0. 


6.786E+0U 

2. 6196+02 

3.1716+15 

6.6216+13 

2.5206-10 

T 179,s22 

0-,97 3 -It 4.55E-0i 

I 967E + 04 

G. 


5 s 60 2 E +0 C 

2. 219 E +0 2 

3.450E + 15 

1.3 '8E + 14 

4.6246-10 

8 1 79r r 3 

0,979 -1. 263E-0.1 

2a C42E + 04 

0. 


3.642 E+OL- 

2.751E+G2 

2.271E+15 

3.4756+44 

9.3826-10 

STEP 18 

2 170.55 

G.-,871 -3;, 426F-(;t 

1.405E+04 

u. 


1.633E+01 

A.034E4u3 

9.1366*14 

2.781E + K 

1.7876-12 

3 179.46 

0.,9iC -2t 743E-0l 

1.387E + 04 

0. 


i«209E+jl 

5. 642E+G2 

1.5686*15 

2.3016+11 

e.0666-12 

4 170,37 

0.932 -2. 343F-0i 

1 . 395E+U4 

u. 


9o40l 6*00 

3. 585E+02 

2.277E+15 

1.1106+12 

2.5246-11 

5 179,28 

0 ,945 -2t 093E-0I 

1., 606E + U4 

0. 


7. 627 E+00 

2. 7586+02 

2.8996*45 

3.6966*12 

6.1866-11 

6 179.19 

C,955 882E-01 

l,643F+04 

0, 


6.279E+J0 

2. 164E+02 

3.3076+15 

9.9226*12 

1.3006-10 

7 179rl0 

0 ,962 -1: 696E-01 

2.U68E+04 

y. 


5.150 E+OU 

2. 10 IE +0 2 

3-6796+15 

2.2986+13 

2.4946-10 

8 179,01 

0.972 -lc394E-0i 

2.144F + U4 

0. 


3.090E+00 

2. 5296+02 

2-5306+15 

6.7366+13 

6.0096-10 

STEP ?4 

? 170.43 

0,£39 -3r 949F-01 

J.2 52E+04 

0 , 


1.442E+01 

8. 1936+02 

1.0616+15 

7.384E+C9 

I. 3716-12 

3 179-34 

0:886 -3c 162E-01 

1.413F+04 

u> 


1.098E+01 

4. 944E+C2 

1.7486*15 

5. 6466+10 

5.660E-12 

4 179- 25 

0,013 -2c 688F-0.’. 

1.405E + 04 

0. 


8.690E+00 

3.2186+02 

2.4756+15 

2.6656+11 

1.6846-11 

5 110.16 

C .930 -2c 380F-01 

i. 418F+04 

0, 


7.081 E+00 

2. 3376+0 2 

3.1106+15 

6.9936+13 

4.0226-aI 

6 170-07 

0.942 -2.123E-01 

i, 7C0E+04 

u. 


5.840E+00 

2.U26E+G2 

3.6036+15 

2.5166+12 

8.4606-11 

7 173 9P 

0.9 52 -Ic 895E-01 

1,747E+04 

0. 


4.7696+00 

4.7146+02 

3.9CiE+i5 

6.3C2E + 12 

1.6436-iO 

8 1.73,89 

vJ 966 -lo 494E-01 

2>2 94E + 04 

0, 


2. 584 E+00 

2. 30 5E +0 2 

2.9146*15 

3.2556+12 

4.5596-10 


• STEP 30 


SWITCH Tfl R-? CnmniNATE SYSTEM- 
trajECTHRY no. 2 06LR=.010 


theta 

RHO 

P Z 

DVDR 

OVDZ 

ZP 

ZPP 

V 

1 2 ?<: 3? 

0. C9C 

0., 075 -0.04 9 

1: 3C8E+00 

-1,7416*00 

1.2556+00 

-5.1656+01 

-8.*56E-C1 

U 4.63 

0.C94 

0.0 85 -0.03 9 

t, 3336+00 

-1.6876+00 

8,5426-01 

-2.0776+01 

-6.2C6E-01 

108.50 

C* ICC 

C. 095 -Or 032. 

1,4076+00 

-i.5C?2E+0o 

6.0726-0 1 

-i. 9746+01 

-8.1836-yl 

lP4r2I 

Oc JCfi 

Or 105 -0.02 7 

1, 3806+00 

-1.614E+00 

4.40 76-01 

-1. 5066+Oi 

-€.il9E-0i 

101, ?5 

c« J 17 

0,115 -0.023 

la 399E+UJ 

-1. 5286+00 

3.1096-0 1 

-l.i20E+0i 

-8.0286-01 

9P.22 

0,127 

0, 125 -C.020 

5 41 06+0 J 

-is 471 E +00 

2s 1236-01 

— 6* 7G0 E+OC 

-1.9376-ui 

P7.83 

Oc 136 

y»‘35 -Or 019 

1,4166 +00 

-1.4326+00 

1,3436-01 

-7.0uiE+OC 

-7.8216-01 

P6.91 

Or 146 

0. 145 -0.01 8 

1,41 96 *-00 

-1.4C7E+00 

7.063E-O2 

-5.794E+CC 

-7,6936-01 

96.31 

0.156 

0,3 55 -0o0l7 

Iw 4216+00 

-1.390E+00 

1.7356-02 

-4. 90 56+00 

-1.5576-01 

95.95 

C>. 166 

Oe 165 -0.04 7 

1: 4ii26+0u 

-1, 38i 6*00 

-2r8196-0 2 

-4.23CE+CC 

-7.4146-Ci 

05.77 

0. ) 76 

0,115 ^CcGJ. 8 

1- 42^e+tio 

-1v376E+00 

-6.7766-02 

-3. 705E+UC 

-7. 2656-Oi 

95.72 

0, 186 

Ocl 85 -Or 019 

l. 422E + 0iA 

-1.374E*00 

-ir 0266-0 1 

r3o2876*G0 

-7.il*E-ci 

95o77 

0,196 

Oi I 95 -0,020 

1, 4226700 

-1, 3766+CU 

-3*3386-0* 

-2c948E+0C 

-6.9536-Ci 

95.90 

0 . 2 c 6 

V: 205 -C, 021 

5, 8y66-ul 

-1.273E+00 

-lc627E-0i 

-2.8O5E+0C 

-6.843E-CI 

9 6.10 

0.2? 6 

0,?15 -0.023 

5- 8l3E-0:, 

-1 247E*wO 

-1.898E-D1 

-2.63iE70u 

-t.762E-Ci 

96.34 

C.227 

0.225 -0.02 5 

5: 8146-01 

-1. 2236+00 

-2a54B-01 

-2.476E70C 

-6.68C6-C1 

96.62 

0.237 

0.235 -0* 027 

5,8116-01 

-i.203E*00 

-2.3946-0 1 

-2. 3316+OC 

-6.594E-C1 

9 6. 93 

0. 247 

Or 245 -C.O3 0 

5: 804E-31 

-lcl84E+0y 

-2.62U-01 

-2.211E70C 

-6.5C6E-C1 



<J7. 

“.„257 

Ck 255 -0: G33 

5,7956-n 

-l,i68Ei-00 

-2,8376-0 i 

-2.U98E+CC 

-6.4166-01 

‘’7,62 

(-•.?63 

3c 765 -0.03 5 

5 784E-01 

-l,153Ei-00 

-3,J4l6-0 1 

-1.994E+C0 

-6.3246-01 

97,90 

0,278 

0.275 -0.w;39 

5, 772E-U4 

-I. 14wE+00 

■-3c2366-01 

-1.90CE+0C 

-6.23U6-01 


C-.2P« 

3, 2 85 -0^.j42 

5, 760E-01 

: -1,1286 + 00 

-3.4226-J 1 

-1.814E+GC 

-6.1356-Ci 

9R.75 

0.. 200 

0. 295 -C.-545 

5. 746E-C/1 

-1.117E+J0 

-3^5996-Oi 

-1.736E+U0 

-6.0386-Cl 

9 0i? e. 

C. 2CC 

(;a3C9 -0c040 

4, 505E-o; 

-1,0366+00 

-3.7736-01 

-1. 7176+00 

-5.9516-Ci 

90,54 

C 2?C 

3.315 -0c053 

4.492E-ax 

-1,0756 + 00 

-3,942E-o; 

-1.6556+00 

-5.8646-Cl 

0 9,04 

U- 230 

Go 32 5 -0,0 57 

4, 473E-Ji 

-1. J65E+0G 

-4.1C4E-0 1 

-1.597E+CC 

-5. 7766-01 

I0r,35 

Cc 341 

0> 335 -0.061 

4..465E-01 

-Ic0566+C0 

-4.26 16-0 

-1. 5446+CC 

-5*6676-01 

Uii\ ->5 

0 : 351 

0. 345 -0,066 

4,45x6-01 

-KU48E+0D 

-4.,4l 36-3 1 

-1.494E+3L 

-5.5576-01 

If'! slO 

0. 262 

;i..355 -C.O'^O 

4; 43 8E-W. 

-5.-5 J40E + CC 

-4r.560E-Ji 

-1.4476+00 

-5.5C66-01 


0. 27 3 

0, 365 -0» J75 

4.425E-Ji 

-1.0336*(;0 

-4.7u3E-0i 

-J,4C4E+0C 

-5.4136-Ul 

K’UOe; 

0. 38 3 

C 375 -0, )79 

4, 4i2F-ui 

-1.0266+00 

-4c94lE-0 i 

-1. 3626+00 

-5.32C6-C2. 

lu2,35 

0 , 304 

0.- 385 -0, J54 

3 727F-J1 

-1. 0 056+Cu 

-4i9766-01 . 

-1.356E+0C 

-5.2266-Cl 

1J2,74 

0. 405 

0, 305 -CnU89 

3, 7x56-01 

* -9.9976-Oi 

-5,1106-01 ■ 

-2.322E+CC 

-5.1396-01 

is'3 4 

tc 616 

C.4C5 -Oc .305 

3: 7C4E-Ji 

-9.943E-0X 

-5t 241 6-3 X 

-l«2916*00 

-5.05C6-01 

10 3, *5? 

U.627 

0.415 — V,. *, ')0 

3,693E-Jl' 

-O„S93E-01 

-5o368E-3l 

-1.2616+00 

-4,9606-01 

]03, 9i 

0,638 

0. 425 -C,i 05 

3: 602E-U1 

-9.84RE-0i 

-5.4036-0 X 

-1.232E+OC 

-4.87CE-CI 

l( ^o.!0 

C , 449 

Co 435 -0.1 11 

3 67i.E-0i 

-9c306E-01 

-5,6156-31 

-1. 206E+CC 

-4. 7786-01 

104,67 

(?., 46u 

0v445 -0.M7 

3: 66 J6-U. 

-9,767E-0'; 

-5 .7346-0 1 

-i.ieiE+co 

-4.6666-01 

?05<,.T4 

C. 47! 

0,- 455 -:%122 

4-t i 4lF-Ji 

-7.9006-01 

-5,8466-01 

-8.3426-01 

-4.5546-Cl 

U. 5,41 

C, 682 

(‘.465 -Oci20 

3, 7246-01 

-7, 768E-01 

-5.9296-Ox 

-a. 36 26-01 

-4, 5C76-UI 

105^77 

C« ^04 

0,475 -3.134 

3. 7i5E-:>x 

-7,7546-01 

-6jUi 16-0 1 

-8.2116-01 

-4.424E-C1 

ic*e,.i2 

0. 5 j5 

0.^85 -Cn4C 

3. 7C6E-J1. 

-7o74ZE-0l 

-6=0936-31 

-8.U67E-01 

-4.34i.E-Cl 

1(16,47 

C.-516 

Oc405 -0.\46 

3, 697E-0, 

-7,73lE-0\ 

-6.1736-01 

-7.9296-01 

-4.255E-01 

J36.,f!o 

0,5 28 

U- 535 -Oo'a 53 

3. 688c -‘Jj. 

-7,7 2;. E-Ol 

-6.2516-01 

-7.7976-Cl 

-4.17CE-01 

33 7. 

0 , 5 30 

Sc 515 -O. i, 59 

3. 6806-Oi 

-7,712E-0i 

-6,3296-01 

-7.67CE-01 

-4.C656-C1 

'10 7q46 

0^5 51 

575 -C.) 65 

3:.67i6-OX 

-7, 7056-01 

-6..4G56-01. 

t7. 5486-01 

-3.5556-01 

’G7,?q 

Of. 567 ' 

0c535 -Oa?2 

3,4i6E-Ji 

' -7, 6216-01 

■-6.480 6-01 

-7.5486-01 

-3.513E-C1 

lOflclO 

0.. .574 

Or 545 -U, 4 78 

3= 40gE-ui 

-7,6286-01 

-6.55 56-01 

-7.4576-01 

-3.3256-01 

lOP,4i 

t B 5 8 5 

C. 555 -Q,185 

3, 4U0F-01 

-7.'635E-01 

-6,6296-01 

-7.3706-01 

-3.745E-01 

:0P 71 

0, 597 

Or 5 65 -0. a71 

3, 391 E-Jx 

-7.642 6-01 

-6c7C2E-01 

-7.2866-01 

-3.66CE-01 

Xr-9,0! 

0, 6C 8 

0. 5 75 -Q,:-,9 8 

3, 383E-OX 

‘-7,6506-01 

-6. 7756-31 

-7.2046-01 

-3.5756-01 

X0 0,3J. 

Ol 62M 

Oc 585 -0.205 

3. 375F-J; 

-7-6586-01 

-6,8476-31 

-7*1256-01: 

-3,4656-Cl 

:-C9.S9 

0, 63? 

Oc595 -0=212 

3: 367E-J1 

-7,6666-01 

-6.9176-31 

-7*0486-01 

-3.4C2E-C1 

!09,fto 

Cc 644 

0. 605 -0.21 9 ■ 

3= 220E-0f 

-7.6396-01 

-6. 9806-01 

-7*0576-01 

-3.3176-Cl 

i3<:.l6 

Or 655 

3. 6x5 -Qc 22 6 

3. 2liE-Ul 

-7r 6646-01 

-7:0586-01. 

-7. JC5E-G1 

-3.2316-01 

UO 44 

0,667 

Or 625 -0. 233 

3; 2ulF-Ui 

-7,6886-01 

-7o)23E-31 

-6.9626-01 ‘ 

-3.145E-C1 

71Cr;7j 

C. 6 79 

0. 6 35 -0:240 

3,1 92E-02 

.-7.7126-01 

-7.1986-01, 

-6*9156-C1 . 

-3.0576-C2 

11C,98 

0. 691 

Oc 645 -0. 247 

3. X82E-01 

-7.735E-G1 

-7, 2676-31 

-6*8706-01 

-2.9706-01 

in. 24 

Oc7C3 

0. 655 -Cc255 

3-. U9E-0. 

-7,7436-01 

-7.3356-31 

-6.8566-Cx 

-2.8626-Cl 

in. 50 

C. 715 

Cc 665 -C, 762 

3:‘07E-0'x 

-7.7856-0! 

-7.4046-31 

-6.8406-tl 

-2.753E-C1 

11 1,76 

G.727 

Cc 675 -0, 269 

3. C96E-01 

-7,8276-01 

-7r4726-0l 

-6,82^6-01 

-2.7C46-C1 

U2.01 

Oc739 

Or 685 -0.277 

3. 084E-0i 

-7^8686-Ul 

-7.5406-31 

-6.8076-Cl 

-2.6146-01 

112.26 

6)c751 

0: 695 -0,2 85 

3: 073E-0X 

-7,ou86-yi 

-7.6086-01 

-6.7906-01 

-2.524E-C1 

1- 2.: 51 

C.76 3 

C-_ 705 -0,2 92 

3-. 0766-01 

-7. 9606-01 

-7.6766-31 

-6.7756-01 

-2.4326-Cl. 

11 2^75 

0.775 

0, 715 -Cf 300 

3, C6iE-ul 

-8cJ2xE-01 

-7.744E-ui 

-6.787E-C1 

-2.34CE-C1 

112:, 09 

0; 788 

0; 725 -C,308 

3, 046E-0i 

-8,08x6-01 

-7.8126^1 

-6.7976-01 

-2.2476-Ul. 

r:3.2 3 

fl, 8CC 

Cc 735 -0. 3x5 

3. G3lF-ui 

-8, 1406-01 

-7.8606-31 

-6,8076-01 

-2.1536-Ci 

U3,46 

0cEl2 

0- 745 -0.323 

3. 01 7E-Ui 

-8= X 97 6-01 

-7,9486-01 

-6,8166-01 

-2.C58E-C1 

113.60 

o, 6 25 

0. 755 -0.331 

3 C72E-OX 

-8,3CflE-Cx 

-8.0166-01 

-6*8216-01 

-i, 9616-01 

113,02 

(.c E37 

0, 765 -0.339 

3- D53E-U1 

-8,3876-01 

-8 3846-01 

-6.8576-01 

-J.064E-O1 

U4,15 

Or E50 

0. 775 -0.348 

3: 034E-Ox 

-8c 4646-01 

-8,1536-01 

-6.8926-Oi 

-1.7656-01 

114.37 

0. F62 

0, 785 -0,356 

3,133E-0t 

-8,5986-01 

-8,2226-31 

-6. 8806-01 

-1.6656-Cl 

1.14 frC 

0- F75 

0.705 -0.364 

3, iiOE-ol 

-8,6986-01 

-8,2916-01 

-6,9426-Ci 

-i. 5626-01 

114.82 

C > ER7 

0:: 8C5 -0.372 

3= C87E-gx 

-8. 7966-01 

-8,36x6-31 

-7.0026-01 

-1.4566-01 

115.0 4 

0 : 5CC 

0. 815 -0.331 

3: 0646 -ux 

-8,9926-01 

-8.4316-31 

-7.U57E-CI 

-1.3536-01 

U'^.75 

Cc 9] 2 

0. 825 -0c389 

3,1 846-v»x 

-950796-01 

-8,5026-01 

-7.0786-01 

-1.2456-01 

U5.47 

Oc 925 

0. 835 -0,398 

3>156E-0i 

-9.1966-01 

-8,5736-01 

-7.159E-01 

-1.1356-Cx 

115,6P 

Cc 938 

0-. 845 -0,406 

3> 285E-ui 

-9,3956-01 

-8.6456-01 

-7.X7S6-CI 

-1.023E-C1 

1' 5.80 

0.951 

Oc B55 -0.415 

3= 253E-0t 

-9, 530E-G1 

-8.7176-01 

-7.28CE-01 

-5,085E-C'2 

1-16,10 

C.963 

0. 8 65 -0.424 

3, 370E-01 

-9,7426-01 

-8.7906-01 

-7.3206-01 

-7.9156-02 

116.30 

(K; 9 76 

0. 8 75 -0.43 3 

3. 334E-J1 

-9.8926-01 

-8.8646-31 

-7.4376-Ux 

-6.7146-02 

116.51 

Cr 989 

0-. 885 -0.441 

3,419E-0I 

-i.OllE+OO 

-8,9386-01 

-7.509E-C1 

-5.4626-02 

V 6. ?1 

ir CC2 

Oo 895 -n.480 

3.3 8UE-J1 

-X. 0276+00 

-9.0146-01 

-7,6346-01 

-4.2276-02 



APPENDIX B 

FORTRAN PROGRAM 


c* 

c 

c 

c* 

c* 

c* 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

r. 


PROGRAM FOR CALCULATING POTtNTlAL FIELD IN AN AXISYMMETRIC 
El ECTROSTATIC CULLECTOR, TAKING SPACE CHARGE INTO EFFECT. 
EXECUTIVE ROUTINE FOR PROCESSING MULTIPLE CASES 
ANO MAIN CONTROL OF LOGICAL FLOW THROUGHOUT PROGRAM. 


COMMON HO.CVF.DRHTI ICJ fOKHMl 16 ) • UTHTD • DTH TR t DVDR ( 10 ) t 
1 UVDTI 10 ) .UVRORtOVROZ .EOLI 1C) « ET AO • IB AR t I T t JF • JT t K POt KPO 1 • 

? KPCISC.I OOP.KEUL.Nl«N2fN3f PHD2A( ID I • P I « PI 2 . P SI A ( 10 ) « PS I 0 ( 10 ) t R t 

3 RAOtRHMI 17) •RHOPI 10) «RHOPP( 10 ) «RHT( ID) «RHVI 10 1 1 SPCH« TAU02 » 

4 T HOI 10) .THMDI 2 5) t THSDEL • THSF AC • TH SWH I t TH SWLO» 

5 THTOI lO) .THTRI 10) • THl , TH2 . V ( I 7 #25 ) ,VQ#VR# 

6 VSMII 10)# VT(IU)*Wll6) •XCTRtXT<25tl6) #Y(16) «ZJO 
EQUIVALENCE (OT H# Y I 1 ))#( TH . Y ( 2 )), I RH • Y I 3 ))» I RH2#RHP,Y(4) ) , 

« iRHlPtYl 5) ) «(kH2P« Y( 6) ) # IMI 1) #OELR) # I w(2 ) «AK) • 

* I W(3) #Z ) # I WI4) «Z2#ZP) f ( W< 51 .ZIP) , I WI6) tZ2P) 

SUHROUTINES- 

1. MAIN 

2. INPUT 
INTGRN 

4. ))E 

5. DERIV 
hm RK 

7. MESH 

8. H/TRAJ 

9. k/RKiRZOE) 

10. R/DE 

11. OUT 

12. VFIFLO 

13. EULINFS 

14 . RHSCAL (A#ZJ0#V0#THl«THF#XtV2#iaAR#RHM#TH2) 

15. PS ICAL Ik# Z«P SI «B0# ASM) 

IF. HRHZ (R#Z#BR#B/fBO#ASM) 

FORMAT I 5H0P ASS# 13# IX • 8HFOLLOWS-) 

initialize space charge FIELD TO ZERO TO SET UP 
LAPLACE ECUATION. 

on 1 i=i#i6 
on 1 U=l.25 
XTIJ# I )=0.0 

INITIALIZE MESH SIZE# BGUNOAKY CONDITIONS# CONVERGENCE 
CR ITER I A. 

pi=o.n 

LKTR=1 
CALL INPUT 
L-LOOP 
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r. 

f. COMPUTE PCTENTIAL FIELO FOR GIVEN SPACE CHARGE FIELD. 

7 CALL VFIELO 

(F( SPCh.EO.n.) GO TO 3 
C 

C CONDITIONAL OUTPUT OF POTENTIAL FIELD. 

MRITE(6.10t LKTR 
3 CALL OUT 
C 

C INTEGRATICN OF ELECTRON CLASS TRAJECTORIES THROUGH 

C THE COMPUTED POTENTIAL FIELD. 

CALL INTGRN 
L*L-1 

LOCP=L OOP-1 
JTsl 

LKTR*LKTR*l 
CALL INPUT 
IFIL.GT.II GO TO 2 
C 

Nls? 

N? = fl 
N3 = l 

SPCH=0.0 

KPOsO 

IFU.EO. II GC TO ? 

C 

STOP 

FNO 


SUBROUTINE INPUT 
C* 

C* READS AND PRINTS INITIAL CONDITIONS. PARAMETERS AND 

C» MISCELLANFOUS DATA REOUIREO FOR INTERNAL USE BV 

C* PROGRAM (OUTPUT FRFOUENCY. LOOPING INDICES. CONVERGENCE 

C« CRITERIA. FTC.I. 

C* 

COMMON BD.CVF.DRHT( 10 I . ORHMU 6 1 .OTHTD . DTH TR . DVOR ( lOI . 

1 OVDTdOl.DVROR.DVRDZ.EULdOI.ETAO.IBAR.IT.JF.JT.KPO.KPOl, 

7 KPGSC.LnOP.NE0L.Nl.N2.N3.PHD2A(10l .PI .P(2.PS1A(10I .PSIOdOl .R. 
3 RA0.RHMI17I .RHOPdOl .KHOPPdO) .KHT( 101 .RHVdO i .SPCH.TAU02. 

A THD( 10I«THMD(25I.THSDEL.THSFAC.THSHHI .THSWLO. 

5 THTOdni.THTRdOI.THl.TH2.V(17.25I.VO.VR. 

G VSM(d0I.VT(10I.Wd6I.XCTR.XT(2S.16) .YI16I.ZJ0 
EQUIVALENCE (DTH. Yd 1 1 .( TH. Y( 21 1 .(RH.Y(3I I . ( RH2.RHP. YI4) ) . 

« (RH1P.Y( 5) I .(RH2P.YI6) I, (Mdl .DELR) . (M(2 I .AR) . 

* (W(3) .ZI.(W(4I .Z2.ZP).(W(5).Z1PI.(H(6) .Z2PI 

C«»« 

IN EORMAT(2X.3Hins.LPE13.3) 

17 FORMAT(2X.3HJOs, 1PE13.3.AX.5HR s.EIO.3) 

Ifl F0RMAT(?X.3HV0*. 1PE13.3.AX.5HBU =.E10.3I 

19 FORMAT(2X.6HOMFGA=, 1PE10.3.AX.5HETAO=.E10.3I 

20 F0RMAT(2X.6FE0 = . IPE 10 . 3. AX ♦ 5HKP =.E10.3I 

21 FORMATdX.GH INPUT-1 

22 F0RMAT(4OI2 I 

23 F0RMATI16H BRILLOUIN FLDM.) 
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24 F0RHAT(/1X*12H0UTPUT E VER Y . I 2 1 IX t6HSTEPS. ) 

25 FQRMAT<2X«6HTAU02^«EIG.3»4X«5HRA0 *«EI0«3) 

26 F0RMAT(16F5.3) 

27 F0RRAT(3XtlH U8X*4HDEL it9Xf4H0ELI. 5X • 5HTHETA f9Xf 3H VSM t 
# 9X«?HV1 «8Xf EHiNJEC.f 8X«4HPSI0) 

28 FORMAT (2Xt I 2 ♦ 2F 1 3* 3 • F 10. 3 • F 12 . 3. IPE 1 4. 2 »0 PF 10. 3 . 1 PEI 5. 3 I . 

29 F0RMAT(16H0 INE TIAL VALUES-) 

30 FORMAT (12X«5FE RAO) »8Xf 5HE0EG) « 4X « 5H E DEG) « 3lX .6HANGLES) 

31 FORMATEIHI ) 

32 F0RMATI55H 

33 FORMATE8E10.3) 

34 F0RMAT(15H CCNF ENED FLOto.) 

35 FORMATE I 2* 10F5.3) 

36 FORMATE/1X»2SHCONE ANGLE IS OUT OF RANGE E=.F5.0t9H DEbREES)) 

37 FORMATE /1X,30HSTEP SIZE TOC LARGE, 0 EL THE TA= ,F5. 2 ) 

38 FORMATE /1X.32HC0NFL ECTING ROUNOARY CONOITEONS-) 

39 FORMAT(2X,18FSURFACE POTENT! AL=*1PE9. 2 ) 

40 FORMAT E 2X • 15FBE AM POTENT I AL= • IPE9. 2) 

41 F0RMATE4F5.1 I 

42 F0RMATE24H0NC SPACE CHARGE EFFECTS) 

43 FORMATE IX, E2,21H ETERATEONS FOR SPACE) 

44 F0RMATE15H CHARGE EFFECTS) 

C 

C 

C ENITEALE7F TRAJECTORY C ALCULAT E ONS. 

C FIRST TIME SETTING OF INITIAL CONOITIQNSf BOUNDARY 

C PARAMETERS, CONVERGENCE CRITERIA. 

JT = 1 

EFEPI.NE.O.) GO TO 9 
C 

r CONSTANTS. 

P 1=3.14159265 
PE2=2.«Pl 
CVF=180./PE 
F0=8.H60E-14 
FTAO=1.753F+15 
0M=1.414E+1C 
C 

C REAC/WRETE. 

WRITE! 6, 31 ) 

RFADE 5,32 ) 

WR ITEE 6, 32 ) 

REAC(5,33) R .RO.DTHTD, THl 

RFA0(S.3B) 7 EO,/JO, VO,DELR,CONFLG.KAO 

RFAC(5,22 ) NT,KPOE ,N1,N2.N3, I BAR , LOOP , KPOSC 

r 

C INCREASE STEP SIZE DTH, BY A FACTOR 

C OF ITHSFAC) EVERY (THSOEL) DEGREES 

C EN THE THETAI5) RANGE FROM ( THSwHI ) 

C TO (THSwLC) DEGREES. 

READ! 5,41 ) TFSFAC,THSOEL,THSWHl.THSwLU 
RFADE 5,26) E VSM I < I), 1=1, NT) 

READE5,35) NFOL, EEJLEE), E=1,NEUL) 

TH2=177.-TH1 /6G. 

JF=25 

KP0=0 

ZKP=Z lU/E VO^^l. 5) 

C 
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r>r> nr> nn no no nnr> 


TFSI MAGNITUDE OF CURRENT DENSITY 
TO SEE IF SPACE CHARGE CALCULATIONS 
SHnULO BF MADE. 

SPCH=0.0 

IF( 7JO.LT.10.0) GO TO 14 
SPCH=l«n 

WKITEI6.43) LOOP 
MR ITEI 6.44) 

GO TO 15 

14 MRITFI6.42) 

15 TAIJ<12=(R*R)/(2.4ETA0*V0) 

INDICATE TYPE OF FLOW, 

IFICONFl n.FC-0.0) GO TO 1 
MR ITEI 6. 34) 

GO TO 2 

1 MR1TEI6.23) 

2 CONTINUE 

READ ROUNCARY CUNOITtONS. 

READI5.26) (XTU.16). Jsl.JF) 
no 3 Jal.JF 

3 XT(J«16)s.xT(J.16) 

XCTRsXTI 1.16) 

no 4 1=1.15 

4 XT! 1. I ) = XCTR 

INPUT DATA CHECK. 

(FI THl.GE.O.C.ANO.THl.LE.174.) GO TO 5 
MRITEI6.36) THl 
CALL EXIT 

5 IF(DThTD.LT.0.0.AND.DTHT0.GE.-.05) go to 6 
MRITEI6.37) CTHTO 

CALL FXIT 

6 IF( VSNI I 1) .LE.ABSIXCTR ) ) GO TO 7 
TO=VSM II NT ) 

MR ITEI 6. 38) 

MRITEI6.3Q) XCTR 
MRITEI6.4C) TO 
CALL FXIT 

INITIAL OLTPUT. 

7 MRITFI6.21 ) 

WRITE I 6. 20) EO.ZKP 
MRITEI6.19) CM.ETAO 
MR1TEI6.16) 710 
MRITEI6.17) 7J0.R 
MR ITEI 6. 18) VO.RO 
MRITEI6.25) TAU02.RAU 
MRITEI6.24) KPOI 
DThTR=DTHTn/CVF 

DETERMINE ENTRY CONDITIONS. 

DIMENSION ALPHI 10) .DELK 10) 
TO=RAO/FLl)ATINT-l) 
no 8 (=1.NT 
T1=FI OATI I-l )*T0 
DEL II I ) = ARSIM Tl/R ) 
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on nr> onnnnn 


ALPHl I ) = 1 .0*R0AT( I 1 
T2=pi-oFL n n 
SI=SIN( T2» 

CT=cns< T2I 
APF=H*Sr 
7E0=R*( l.O+CT ) 

CALL PSrCALI ARF *2FU *PS I . BO . RAO ) 

psio( r >=ps I 

IFfCONFLQ.EO.O. ) PS[0(I)s0.0 

8 CONTINUF 

SFTTtNG OF RHO-THETA VARIABLE MESH. 

CALL MFSH 

SPREAD INITIAL TRAJECTORY VALUES. 

EVFRY TIME INITIALIZATION OF TRAJECTORIES. 

9 on 10 IT=2,NT 
RHTI IT » = .999S 
ALPHA = ALPH( IT) /CVF 
GAMMA=PI-ALPFA-OEL II IT)+DTHTR 

T2=-RHT I IT I /( DTHTR*SINIGAMMA) » 

T3=SIN( ALPHA+OEL I( IT ) )-S IN (GAMMA) 

RH2P=-T2*T3 
ORHT( IT»=RH2P 
RHOPI IT) =RH2F 

10 THO( IT ) = VSMI ( IT)/( 1. +RHOP ( I T ) **2 ) 

CALCULATE VT(U. 
no 11 1=1. NT 
THTR( I ) = PI-CCL I (I ) 

THTD( I ) = THTR ( I )*CVF 

11 VT (I )=V0*VS)<!I (I) 

PRINT INITIAL TRAJECTORY VALUES. 

URITE(6.29) 

WR ITE(6.27) 

MRITF(«.30) 

DO 12 1=2.8 
DEL ID = nEL 1 1 I )«CVF 
T1=180.-DEL IC 

IP WRITE (6. 28) I .DELK I ) .DEL ID .T1 . VSMI ( I ) . VT( I ) . ALPH (I ) , PS 10 ( I ) 

on 13 1=2. NT 
IT=I 

th=thtr( it ) 
rH=RHT( it ) 

RHP=RH0P(IT) 

CALL DE 

13 RHQPP( IT ) =RH?P 
C 

RETURN 

END 
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SURRCUTfNE IMGRN 


€♦ COMTRGtS INTEGRATION LOOPING OF EQUATIONS OF MOTICNt 

C* STEP-SIZE, CHANGEOVER OF COORDINATE SYSTEMS AND THE 

€♦ NUMBER OF ITERATICNS REQUIRED FOR CONVERGENT SOLUTION. 

C« 

COMMON RO,CVF,DRHT( 10I,DRHM( 16 ) • OT HTD • OTH TR » DVOR 1 101 t 

1 OVOTI 10) •DVROR,DVRDZ,EQLI 10) tETAO,IBAR,ITf JF, JT.KPO.KPOI t 

2 KPCSCtLOnPtKEQL#NlfN2,N3,PH02A(10) ,PI • P I 2 • P SI A( 10 ) • PS i 0 ( 10 ) • R , 

3 RAO,RHM( 17) tRHOPC 10) tRHOPPIlO ) ,RHTI 10) »RHV( 10 ) t SPCH « T AUO 2 , 

A THO( 10) .THMCI 25) • THSOEL , THSF AC, THSWH I • THSWLCt 

5 THTOI 10) ,THTR( 10) ,THl,TH2tV( 17,25) • VO, VR , 

6 VSMK 10),VT(10),W(16),XCTR,XT(25,16),Y(16) ,ZJC 
EQUIVALENCE (OT H, Y< 1 ) ) , < TH , Y ( 2 ) ) , ( RH , Y ( 3) ) , ( RH2 ,RHP, Y(A) ) , 

* (RH1P,YI 5)1 ,IRH2P,Y(6) ) , Ui ( 1 ) ,DELK) • (M(2) ,AR) , 

^ IW<3),Z),I)il(A),22,ZP),lU(5),ZlP),IW(6),22P) 

IG F0RMAT<7H0DTFTD=,1P69.2) 

EXTERNAL OE 
OTF=:DTHTR 
f. 

C CHECK FOR SPACE CHARGE ITERATION* 

IFI SPCH.EO.O.) GO TO 1 
Nls5 
N2 = 5 
N3 = l 

IFILSM.EQ.9) GO TO 1 

LS)» = 9 

TH0=THT0(5) 

C 

1 DO 2 IT=N1,N2,N3 
TH sTHTRl IT ) 

RH sRHTC IT) 

RFP=DRHT< IT) 

RCI n=RH 
CALL RK(OE) 

RHOP ( IT )=RHF 
RHOPPI IT )=RH2P 
THTRI ITi=TH 
THTCI IT)=THTRI IT)*CVF 
RHT( IT )=RH 
? DRHTI IT)=RHP 
C 

C INCREASE STEP SIZE DTH, BY A FACTOR 

C OF (THSFAC) EVERY (THSDEL) DEGREES 

C IN THE THETA(5) RANGE FRGM (THS)«lHl) 

C TO ITHSWLO) DEGREES. 

IF( THTOI 5) .GT.THSNHI ) GO TO 3 
IFI THTD(5).LF.THSttLO) THSWHI^O.O 

THSRHI=THSWH I- THSDEL 
DTH=THSFAC*DTH 
T1=CVF«0TH 
WRITEI6,10) T1 
C 

3 CONTINUE 
JT=JT^1 
CALI OUT 

IFCSPCH.EQ.O*) GO TO 4 
IF(HH.GT.RHM(7) ) GO TO 1 
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THF=THTn(5) 

CALL KHSCAL(RAO*ZJO»VO.THO.THF.XT, VSMl . IBAR* RHH» THMO ( 1 ) t KPOSC ) 
KPC=0 
RFTIJRN 

r 

f. R-7 COORDINATE TEST. 

A IF( THTR(2).DT.3.1) GO TO 5 

If <RHnP(N2).LT.O.I CALL RZTRAJ 

5 IFLRHTINl I.GT.1.0) GO TO 6 
GO TO 1 

C 

6 RETURN 
END 


SUBROUTINE CE 
C* 

C* CONTAINS FOUATICNS OF MOTION IN SPHERICAL 

C* COORDINATE STSTTEM. 

C* 

C*** 

COMMON BO.CVF.ORHTI 10 I ,ORHM< 16 I ♦ OTHTD , OTH TR . OVDR 1 10 ) t 

1 UVOTI 10) ,OVROR.UVKOZ.EUL( 10) . ET AO . 1 BAR « I T , JF , JT t K POt KPOi t 

2 KP0SC.L(j0P*NE0L.Nl.N2.N3*PH02A(10) tPI .PI2.PSIAI10) tPSIOI 10)f R« 

3 RAD»RHM( 17).RH0P( 10).RHOPP(10),RHT(10),RHV(10)#SPCH,TAUtJ2, 

4 THO( 1C) .THMCI P5).THSD£L»THSFAC.THSWHI .THSWLO. 

5 THTOI 10)»THTR( 10) .THl,TH2«V(17.25)»VO,\/Rt 

6 WSMni0).VT(10).N(16).XCTR«XT(25*16)«Y(I6) *ZJO 

EOUI VALENCE IDTHt Y ( 1 ) ) * ( TH . Y ( 2 ) ) . ( RH . V ( 3 ) ) . ( RH2.RHP* Y(4) ) t 

* (RHIP.Y(S)) .(KH2P« YI6) ) • I m( 1) «OELK) « <M(2 ) tAR) t 

* (U(3) .7) «( R(4) .22.ZP).(M<5) .ZIP). (W(6) ,Z2P) 

C*** 

SI=SIN<TH) 

ST2=ST**2 

CT=COS(TH) 

C 

C COMPUTE nv/ORHO AND DV/OTHETA. 

CALL OERIV 
ARE=«H*R#ST 
7ED=R»( I.+RH^CT ) 

C 

C COMPUTE PSI. 

CALL PS I CAL ( ARE .ZED. PS I. BO. RAD) 

C 

C COMPUTE BR AND BZ. 

CAL) BRHZ) ARE. ZED. BR.BZ. BO. RAO) 

T1 =-FTAO/( ARE«ARE*P 12) 

PHD=Tl*( PSI-PSIO) IT) ) 

PHD2=PHO»*2 
PHD2A) IT )=PHC2 
PSIA < IT)=PS I 
C 

C COMPUTES FIRST AND SECOND DERIVATIVES OF RHO WITH 

C respect TC theta (SPHERICAL COORDINATES). 

DVMDRH=R*( HR4ST+8Z4CT ) 

0VMDTH=HH*R* (BR*CT-RZ*ST ) 

Tl=RH*RHfR2*R<‘ST4PHD 
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n 


PAR1=DVDR( IT )/R + ST*PHO*OVMOTH 
PAR2=0V0TI IT )/R-Tl*OVMORH 

TO=(VSMI( IT l+VTI IT) ) / 1 T AUC2*RH*RH« ( I .♦ ( RHP/ RH M>*2 ) ) 
THn2»TO-(PHD5*ST2) /( l.*-IRHP/RH)**2) 

ThOI IT ) = ThD2 
T1=RH*ST*PHC 
T2 = ( Tl*ST*PI-0» /THD2 
T5=PTAO/R 
T3=TS/Thn2 
T1=RH+T2+T3*PAR1 
T2=-( RHP*ST=»CT*PHD2 ) /TH02 
T3=( 2.*RHP**2) /RH 
T4=-(RHP*T5 )/( RH*RH«‘THD2) 

RH2PaTl*T2tT3+T4*PAR2 

RHlP=RH2 

f, 

RETURN 

END 


c* 

.SUBROUTINE DERIV 




c* 

CALCULATES TI'F RATES OF CHANGE 

OF 

THE 

POTENTIAL 

f ♦ 

FIELD IN SPHERICAL COORDINATES 
INTERIOR POINT. 

FOR 

ANV 

SPECIFIED 


nOMMON bOtCVf .UKHTi 10) .DRHMI16) fOTHTOf OTHTR.DVDR(IO) f 

1 DVOTI 10) .OVRDRfOVRUZ tEOL ( 10 ) . ETAO . I b AR . 1 T t JF t JT t KPQ f KPO I « 

2 KPnSC.Ln0PtNEUL«Nl.N2fN3«PH02A(lG)*PI.P12f PSIA(IO) tPS10(lC)«Rf 
? RAO.RHMI 17).RH0P( 10 ) « RHUPP ( 10 ) *KHT ( 10 ) .RHV ( 10 ) » SPCH . T AU02 i 

4 TFOI lO.THMOI 2 5) , THSDEL . THSFAC, THSWHI .THSWLO. 

5 THTDI 10).THTR( 10) ,TH1 .TH2 . V( 17.25 ). VO t VR . 

6 V.SMK 10 ), VTdO ).W( 16) .XCTR.XTI25.16) .y(16) . 2 JO 

EOUl VALENCE I OT h. Y ( 1 ) ) , ( TH . Y ( 2 ) ) . ( RH , Y ( 3) ) , I RH2.RHP, Y(4) ) , 

* (RHIP.YI 5) ) .(KH2P.Y(6) ) . (R(l ) .OELR) . (M(2 ) .AR) . 

* <M(3).2).IM(4),22,ZP).(M(5).Z1P).(M(6).Z2P) 

C.*** 

T2=TH*CVF 
T3=15. 

DETERMINE RHO INTERPOLATION FACTOR (FR). 

DO 1 1=1.17 

IFIRHMI D.OT.RH) GO TO 2 

1 CONTINUE 

2 I F( I.FO. 17 ) 1=16 

FR = (RH-RHM( 1-1) )/DRHM( 1-1) 

f. 

C DFTERMINF THETA INTERPOLATION FACTOR (FT). 

DO 1 .1 = 1.25 

IF(THMD( J).GT.T2) GO TO 4 
^ CONTINUE 

4 IF(T2.GT.TH/ ) T3 = .25 
FT=(T2-THM0( J-1 ))/T3 
T3=T3/CVF 
C 

VPl=V( I-l. J-1) ♦FT*( VI I-l.J )-V( l-l. J-1) ) 

VP2 = V(l . J-1)+FT0( VI I .J )-V(I .J-D) 
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VP3 = V(l-l.J «J )) 

VP4 = V( l-lt V( I 

f. 

VT< m=VPl+FP*( VP2-VP1I 

f. FVALUATFS OV/QRHO AND OV/OTHETA AT (RH.T2) FROM CURRENT 

c potential field. 

OVORI IT I =11 VF2-VPU /ORHMI I-l ) I*VG 
OVOTI IT) = ( ( VF3-VP^»/T3H‘V0 
IFITHTDI IT».L-T.TH2) DVOT11T)=0.0 
f. 

RETURN 

END 


subroutine RK(OE) 

C* 

C INTEGRATES ONE TRAJECTORY ONE STEP USING RUNGE-KUTTA 

C FORMULAS IN SPHERICAL COORDINATES. 

f.* 

C*** 

COMMON HO.CVF, DRhTI 10) .ORHMI 16 ) . OTHTO . OTH TR . OVOR 1 10) . 

1 OVDTI 10 ) ,DVROR ,0VR02 .EOL I 10) . ETAO.l BAR .IT. JF, JT.KPO.KPOI . 

2 KPnSC.LnOP.NEaL.Nl.N2.N3.PHD2Al 10 ) . P I . P I 2 . P SI A 1 10 ) . PS 1 0 1 10 ) . R . 

3 RAO.RHMI 17) .RHOPI 10) .RHOPPI 10) .RHTI 10) .RHV 1 10 ) . SPCH , TAU02 , 

A THQIUn,THMCT25),THS0EL.THSFAC.THSWHl,THSWLQ. 

5 ThTDI 10) .THTRI 10) .TH1.TH2.VI 17.25). VO.VR. 

6 VSMK in ) . VT 1 10 ) .Ml 16 ) .XC TR. XTI25 .16 ) . V 1 16) . ZJO 
EOUIVALENCE ( DT H. Y ( 1 ) ) , I TH . Y ( 2 ) ) . I RH , Y 1 3 ) ) . ( RH2 . RH P, Y ( 4 ) ) . 

* (RHIP.YI 5) ) .(RH2P. Y( 6) ) . (Ml 1) .OELR) . I Ml2 ) .AR) . 

« (MI3) .2 ) .1 m(4) .72.ZP).(M( 5) .ZIP). (M(6) «Z2P) 

C.** * 

Y(15)=Y(3) 

Y(16)=Y(4) 

CALL DE 
Y( 7)=Y( 1 )*Y( 5) 

Y( B )=Y( 1 )4Y46 ) 

Y(2)=Y(2)+Y( l)/2. 

Y(3)=Y( 15)+.5«YI7) 

YI4)=YI16)».5*Y(8) 

CALL DE 

Y(q )=Y( 1 )4Y4 5) 

Y( 10)=Y( 1 )*Y(6) 

Y<3)=Y(15)*.5#YI9) 

Y(4) = Y( 16M-.54YI 10) 

CALL DE 

Y(1 1 ) = Y1 1 )*Y (5) 

Y( 12 )=YI 1 )*Y (6) 

Y(2)=Y(2)+Y(l)/2. 

Yf3}-Y(15t*ym ) 

YIA )=Y(16)*YI12 ) 

CALL DE 

Y( 13 )=Y( 1 )*Y(5) 

Y( 1A)=Y( 1 )«Y(6) 

C 

Y(3)=Yll5)+( YI7)«-2.*YI 9 ><-2.* Y 1 11 > Y (13 > ) /6. 

Y(A)*YU6) + ( Y(8)+2.*Y(10H-2.*Y(12)*Y(1A)) Z6. 
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CALL DE 

return 

END 


SUBROUTINE MESH 
C* 

C* FOR ANT SPEC lEIED COLLECTOR GEOMETRY, GENERATES 

C* GRADUATED MESH POINT ARRAYS IN BOTH THE RHQ AND 

C* THETA DIRECTIONS. 

C* 

C.*** 

COMMON BO.CVF.ORHTI 10 I ,DRHM( 16 ) , OTHTO, DTHTR , UVDRI 1C I , 

1 OVDTI ID ) .UVRDR.QVRDZ ,EUL( ID I , ETAO. I BAR, I T, JF, JT,KPO,KPOI , 

?. KPnSC,LnOP,NEOL,Ni,N2,N3,PHD2AI 10 ) . P I , Pi 2 , P SI A( 10 ) , PS 1 0 ( 10 ) , R . 

3 RAD,RHM( 17 > ,RH0P( 10) ,RH0PP(10) ,RHT( 10) ,RHV I 10 » , SPCH , T AU02 , 

4 thOI 10) .THMCI 25 ) ,THSOEL ,THSFAC, THSmHI ,THSRLU, 

5 THTOI lOiTHTRI 10 ) , Th 1 ,TH2 , V 1 1 7 ,25 ) , VO , VR , 

is VSMII ID), VT 110 ) ,MI 16) ,XCTR,XT 125 ,16) , YI16) , 2J0 
EOUIVALFNCE lOTH, Y I 1 ) ) , I TH , Y I 2.) ) , I RH , Y I 3) ) , ( RH2 , RHP , Y 1 4 ) ) , 

* IRHIP, YI 5) ) ,IRH2P, Y( 6) ) , (WI 1 ) ,l)ELR ) , I W ( 2 ) ,AR) , 

* lHl3),Z),IWi4),22,2P),(M(5),ZlP),lWI6),Z2P) 

C*** 

C GENERATE RHO AND DELRHO ARRAYS. 

Tl=n. 8/120. 

KHMI 17 ) = 1 .0 

GENERATE variable MESH RHO IN RHM ARRAY RANGING FROM 
RHMIl)=n TO RHMI17)sl WITH FINE STEP AT END. 

no 1 1 * 1,16 

J=17-I 

DRHMU )*FLOAT( I )*T1 

1 RHM I J)*RhM(J+l )-ORHMIJ) 

0HHMI1)=0.P 
RHM I I )*0.0 

GENERATE THETA ARRAY. 

JR*TH1/15. 

JKs12-JK 

GENERATE VARIABLE MESH THETA IN DEGREES IN THMD ARRAY 
RANGING FROM THMD(1)»TH1 TO THMDI 2 5 ) *1 80 WITH FINE STEP 
AT END. 

no 2 .1 = 1 , JX 

2 THMDI.) )*TH1*15.*FL0A7IJ-1) 

.)Kl*JKtl 
JK2*JX+2 

THMDI .)K1 )=17 7.-.2 5*FLOATI JR) 

TH?=THMDI JKl ) 
no 3 J=JK2, JF 

3 THMDI J)=THM0IJX1)*.25«FLDATIJ-JK1) 

RETURN 
END 
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SUBROllTtNF R2TRAJ 


C* 

C* CONTROLS INTEGRATION OF ECUATIONS OF MOTICN IN THE 

C* R-7 RECTILINEAR COORDINATE SYSTEM. 

C COVERT FROM SPHERICAL COORDINATES (RHO. THETA) TO CYLINDRICAL 

C COORDINATES (R.Z) AND CONTINUE INTEGRATION OF TRAJECTORY UNTIL 

C IT REACHES SPHERICAL SURFACE. 

C* 

C.*** 

COMMON BO , CV F, DRHT 1 10 ) .DRHF 1 16 ) .DTHTD. OTHTR . DVDR 1 10) , 

1 OVDTI 10) .DVRDR.0VRO7 .EULI 10) . ET AO . I BAR . I T . JF . JT .K PO . KPO I . 

;> KPOSC.LOOP. NEUL.N1.N2.N3.PH02AI 10 ) . P I . PI 2 . P SI A( 10 ) . PS I0< 10 ) • R . 

3 RAO.RHMI 17 ) .RHOPI 10) .RHOPP(IO) .RHTIIO ) .RHVI 1C ) * SPCH.TAUG2. 

A THOI 10) .THMCI 25).THSDEL,THSFAC.THSRHI.THSRL0. 

5 THTDI 10) .THTRI 10) . THl , TH2 . VI 1 7 ,25 ) . VO , VR , 

6 VSMI I 10 ). VT (1C ) ,W( 16) .XCTR.XT(25.16) . YI16) ,2J0 

EUU I VALENCE IDT F. Y I 1 ) ) . I TH . Y I 2 ) ) . ( KH . Y ( 3 ) ) , I RH2,RHP,Y(A) ) , 

« (RHIP.YI 5) ) .IRH2P.YI6) ). I W( 1 ) ,0ELK) . (M(2) .AR) , 

* (UI3) .Z ).( R(A) .Z2.ZP).( W(5) .ZIP) ,(M(6) ,Z2P) 

C*** 

10 F0RMAT(33HlShlTCH TO R-Z COORDINATE SYSTEM-) 

11 F0RMATI15H TRAJECTORY NO. , I 3, 2 X , 5HDELR = , F A. 3 ) 

12 FORMAT! /AX,5FTHETA,3X.3HRH0,5X ,1HR.7X,1HZ,6X ,AH0V0R,8X,AH0V0Z.9Xt 

♦ 2HZP.I0X.3HZPP.10X.1HV) 

13 F0RMATI2X.F7.2,3F7.3,1P5E12.3) 

EXTERNAL RZCE 
MRITE(6.10) 

WRITE(6.11) N2.DELR 
URiTE(6.12) 

CONVERT TC R-Z COORDINATES. , 

AR=*RHT(N2)#S INI TH) 

Z =RHT(N2 )*CCSI TH) 

INITIALIZE R-Z INTEGRATION. 

TO=RHOP( N2 ) *Z-RHT( N2 )*AR 
T1=RHT(N2 )»Z<-RH0P<N2)#AR 
ZP=T0/T1 
CALL RZOE 

THRZ=90.-CVF*ATAN( Z/AR) 

KHRZ = SORT( AR«AR^Z'«‘Z) 

T2=DVR0Z-ZP*CVR0R 
TO = VSMI I IT) +VR 
T1=1.*ZP**2 
T3=2.*( TO/Tl) 

Z2P=T2/T3 

WR(TE(6.13) THRZ,RHRZ,AR,Z.OVROR,UVROZ ,ZP,Z2P, VR 

INTEGRATE TO SURFACE. 

1 CALL RZRKIRZCE) 

THRZ=90.-CVF*ATAN(Z/AR) , 

RHRZ=SORT( AR*AR+Z*Z ) 

MR ITEI6, 13 ) THRZ.RHRZ. AR .Z.OVROR.OVRDZ ,ZP ,Z2P, VR 
IF! RHRZ.LT.l .0 ) GO TO 1 
IF(N1.EQ.N2) TO^TO/O.O 

N2=N2-1 
RETURN 
END 
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SUfiRCUTiNE R7KK(RZUE) 

C* 

C 

C INTEGRATES ONE TRAJECTOPy ONE STEP USING RUNGE-KUTTA 

C FORMULAS IN CYLINDRICAL COORDINATES. 

C* 

C*** 

COMMON 60.CVF .ORHTI 10 J ,ORHM ( 16 ) .OThTD . OTH TR . DVOR ( 10 ) f 
1 nVDTIlOi.OVRDR.OVROZ.EOLllOI.ETAO.IBAR.IT.JFf JT.KPO.KPOIi 
?. KPOSC.LOOP*NEUL*Nl«N2.N3.PHD2A( 10 I . P I . P I 2. P SI A ( 10 ) * PS 10 ( 10 ) , R 
3 RAO.RHMI 171 .RHOP( 10 ) . RHOPP ( 10 ) «KHT ( 10 I .RHVI 10 ) .SPCH.TAU02 . 

A THD( 10) .THMC(25)«THS0EL.THSFAC.THSMHI«THSML0. 

5 THTDI 10) .THIRI 10 ) .Thl ,TH2 . VI 17 .25 ) , VO. VR . 

6 VSMK 10 ).V1 (10) .W( 16) .XCTR.XTI25. 16) .Y(16) .ZJC 

EOU I VALENCE (OTH.Yl 1) ) . ( TH. Y ( 2 ) ) . ( RH . Y ( 3) ) . ( RH2.RHP. YIA) ) . 

* (RHlP.Y(5n .(RH2P.Y(6> ). (M(l) .OELR) . (W(2) .AR) . 

■* (R(3) .7 ) .( W14) .22.ZP).( M( 5) .ZIP) . (W(6) .Z2P) 

C*** 

W(15)=W(3) 

W( l6)aW( A) 

CALL R7DE 
M(7)sm(1|«M{5) 

W(ft)>w(l )«W(6) 

N(2 )sM(2 )fW( 1 )/2. 

M(3)=N(15H-.S«M(7) 

W(A)sU(16)4-.5«W(8) 

CALL 870E 
M(9 )sM(l )«M( 5) 

M( 10 )=M( 1 )*M (6) 

W(3 )sm(15)-».5«W(9) 

U(A)=M( 16 )'«-.5«W( 10) 

CALI R70E 

M( 11 )sM( 1 )«W (5) 

WU2)-M( 1)«W(6) 

W(2)=m(2)^w( 1)/2. 

w(a)=M(i5)«w(iii 

M(A)=W(16)4^M(12) 

CALL R70E 

W( 13 )=M( 1 )«R (5) 

M( lA )=R( 1 )«W (6) ■ 

M(3)=M(15)«( R(7)^2.«W( 9)+2.*W(11H^m(13) )/6. 

W( A )=W( 16) to( 8) 4^2. «M( 10) «^2.«U< 12) ^M(l A) ) /6. 

CALL R7DE 

RETURN 

FNO 


subroutine R7DF 
C 

C COMPUTES FIRST AND SECOND DERIVATIVES OF 7 wITH 

C RESPECT TC R (CYLINDRICAL COORDINATES). 

C* 

c*** 

COMMON HO.CVf.ORHTI 10) .DRHMI 16) .UTHTO.DTHTR.DVDRI 10) . 

1 OVDTI 10 ) .OVRDR .OVROZ .EOL ( 10) . ETAO . I BAR . I T . JF . JT . KPO. KPOI . 

2 KPnSC.LOOP,NEUL.Nl.N2.N3.PH02A(10) ,PI.PI2.PS1A(10).PSIO(10),R 

3 RAO.RHMI 17) .RHOPI 10) .RHOPPI 10 ) .RH T ( 10 ) .RHV ( 10 ) . SPCH. T AUO 2 . 



4 ThO< UH ♦ THMC( 2 5) . THSDEL . THSF AC . THSWH I . TH SWLO# 

5 THTl)( 10 ) . THTR( 10) .TH1,TH2.V( 1 7 .25 ) . VO . VR ♦ 

A VSMI ( 10 ) . V T ( IC ) .t«(16 ) «XCTR. XT<25.16 ) . Y(16) * 2JC 
ECU I V Al.ENf E lOTE.Yt 1 ) ) .( TH. Y( 2 ) ) .( RH , Y( 3) ) . ( RH2 .RHP, Y(4) ) • 

* (RHIP. Y( 6) ) .(RH2P. Y(6) ) . (R( I) .UELR) . (m(2) ,AR) , 

* (W(3).2).(t.(A).Z2.ZP).(W(5).ZlP).(M(6).Z2P) 

C*** 

RH=SORT( Z**2 ♦AR^^Z ) 

TH=ATAN2( AR.Z) 

CALL HER IV 
VR=VT< IT ) 

DVRDR = < AR*OVCR( I T ) * ( Z*UVOT ( I T ) )/RH ) / ( RH#VO) 

DVROZ =(Z»OVDR( I T ) -AR*CVDT I IT ) /RH ) / I RH* VO) 

T1=DVR0Z-ZP*DVR0R 

T2 = 2.*( IVSM H m*VK)/(l.>ZP**2) ) 

Z2P=Tl/T2 

Z1P=Z2 

C 

RETURN 

ENO 


SUBRniJTINE nuT 
C* 

C PRINTS OUT SELECTEO DATA ALONG A TRAJECTORY. 

C* 

c*** 

COMMON BO.CVF. DRHTI 10) .DRHMI 16 I.DTHTD.OTHTR.OVDRI 10) , 

1 OVOTI 10) .OVRUR .OVROZ ,EOL( 10) .ETAO, IBAR ,I T. JF, JT.KPO.KPOI , 

2 KPOSC.LnnP.AE0L.Nl.N2.N3.PH02A( 1C ) . P I , P I 2. P SI A 1 10 ) . PS i 3 < 10 ) . R . 

3 RAD.RHMI 17 ) .RHOPI 10 ) .RHOPP I 10 ) .RHTl 10 ) .RHV I 10 ) , SPCH, TAUD2 , 

4 THOI 10) .THMCI 25) . THSDEL . THSF AC . THSWH I . TH SWLO. 

5 THTDI 10) .THTRI 10 ) . THl .TH2 . V ( 1 7 .25 ) . VO . VR , 

6 VSMK 10 ) . VT 110 ).W( 16) .XCTR.XT(25.16) .YI16) ,ZJO 

ECU I VALENCE (OTR. Y ( I ) ) . I TH . Y< 2 ) ) . I RH . Y I 3) ) . I RH2.RHP.YI4) ) , 

* IRHIP.YI 5) ) .(RH2P. Y<6) ) . I W( 1) .UELR) . I WI2 ) .AR) . 

* (Wl 3) .Z ) , I W (4) .Z2.ZP ). I W( 5) .ZIP) . I WI6) .Z2P) 

C*** 

C FORMAT STATEMENTS. 

10 F0RMAT(10H0DELTHETA=.F7.4.5H DEG.) 

1 1 F0RMATI5H STEP, 15) 

12 FORMAT! /1X.3FTRJ.2X.5HTHETA.2X.3HRH0.6X, 

A 2HVT.7X. 5HCV0RH.7X.5H0V0TH.6X,5HRH-PR,6X.6HRH-PPR,6X,4HTHD2, 

8 6X.4HPH02.6X.3HPSI ) 

13 FORMAT! IX . I 2 .F8.2, F6.3 , IP 8E1 1 .3) 

C 

C print TITLES. 

IFIK9.E0.9) GO TO 1 
K9*9 

WRITEI6.10) CTHTO 
WRITEI6.12) 

C 

1 IF IKPO.GT. 1 ) GO TO 4 
JT1=JT-1 
WRITEl6.il) JTl 
no 2 I=N1.N2.N3 
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2 UR[TE(6.13) i,THTO( ( ) «RHT( II t VT( I) .OVOR( I ) t 

♦ OVDTC I I .Rt-OP( I I .RHOPPI I I .THD( 1 1 ,PMD2A ( II # PS I A ( I) 
IFIRHT(N2 I.G7. 1,(11 (iO TO 3 
C 

KPC=KPOI 

RETURN 

3 N2=N2-l 
KPC=KPOI 
RETURN 

A KPG=KPO-l 

IE(RHT(N2*.GE.1,0) GO TO 4 

RETURN 

END 


SUBROUTINE VflELO 
C* 

C 

C COMPUTES THE POTENTIAL EIELU V EOR A GIVEN SPACE CHARGE 

C FIELD 8Y SOLVING POISSONS EOUATION. 

C 

C POISSl ANC P0ISS2 ARE SUBROUTINES DESCRIBED BY THE TN 

C OF REFERENCE 2. THEY HERE DEVELOPED TO SOLVE POISSONS 

C EOUATION IN A SUBREGION OF A SPHERE IN AN EFFICIENT 

C NON-ITERATIVE FASHION. USING A BLOCK DIAGONAL MATRIX 

C TECHNIQUE. 

C* 

COMMON/BLOCK l/NR.NT.RAI 17 I .TA ( 25 I ♦ X ( 25 . 16 I ,X 1 . AVS 

C*** 

COMMON BO . CV F . OR HT ( 10 I , ORHM I 16 ) , DTHTO . OTH TR . OVOR 1 10 I . 

1 DVOTI 101. 0VR0R.0VR02.E0LI1C) .ETA0.1BAR.it. JF.JT.KPO.KPO l. 

2 KPOSC.LOOP.NEOL.N1.N2.N3.PHD2AI 10 I . P I . P I 2 . P SI AIlO 1 . PS 10 1 10 ) , R , 

3 RAD.RHMI 17 I .RHOPI 101 .RHOPP 1 10 1 .RHTI 10 I ,RHV( 10 I . SPCH. TAU02 . 

4 THO( 10) . THMC( 25) .THSDEL. THSFAC.THSWHI .THSWLO, 

5 ThTD( 10).THTR(10).TH1,TH2.VI17.25).V0.VR. 

6 VSMK in ). VTIIO ) .MI16) .XCTR.XT(25.l6) .Y116) . 2J0 
EQUIVALENCE (DTF. Y( 1 ) ) . I TH. Y I 2 ) ) .<RH.Y(3) ) . I RH2 .RHP, Y(4) ) . 

* (RHIP.YI 5) ) .(RH2P. Y( 6) ) . I Wl 1) .OELR) . (W(2) .AR) . 

« (W(3) .Z).IHI4) .Z2.ZP).(U(5) .Z1P).(W(6) .Z2P) 

C4** 

LOGICAL AVS 

DIMENSION XBa5).TD(25) 

r 

23 FORMAT! 16F5, 3) 

74 FORMAT! 1X,F6.2.1P12E10,2» 

25 F0RMAT!7H INPUT-. 52X. 3HRH0) 

76 FORMAT ! 1 1H04»N0 BAR*'*) 

27 FORMAT! 13H0 EAR VALUES-) 

2fl F0RMAT!4X,1HR.8X.1HV) 

2<5 FORMAT! 2X,F5.3.3X,F6.3) 

30 FORMATIBH OUTPUT- , 52X . 3HRH0 ) 

31 F0RMAT(6H THET A , 3X . 16 !F5 . 3 , 2X ) ) 

32 F0RMAT!6H TF ETA , 3X. 12 1 F 5. 3 . 5X ) ) 

33 F0RMAT!4H Xl=.F5.2) 

34 F0RMAT!1X,F6,2.16F7.4) 

35 FORMAT! Z54X. 15HP0TENT lAL FIELD) 



n r> r> r> n r> n n n r> r> n n n 


X1=XCTR 
NR = I7 
NR1=NR-1 
NR2=NR-2 
on i 1=1 • 17 

1 RA( n=RHH( I I 
NT = JF 
i^Tl=NT-l 

DG 2 .1=1 
00 2 K=ltNRl 

2 X(J,KI=XT( Jf K1 

If AXIAL SPIKF IS PRESENT IN THIS PARTICULAR PROBLEMf 
READ IN THESE VALUES AND DISPLAY THEM. ^ 

If I IRAR.EU.O ) GO TO 4 
00 T f=l,NR2 

3 XINT, I ) = XBI I ) 

4 on 5 .l=l*NT 
TCI J)=THM0IJ } 

5 TAI J) = THMCIJ )/CVf 
If I IPAR. EO.O ) GO TO 9 

READ 8AR VALUES. 

AVS=.TRUE« 

IFILSm.EU.91 go to 7 
LSM = 9 

READ(S«231IX(NTtIl«I=lfNRll 
00 6 f=l*NR2 
Xfil n=*XINTt II 

6 X4NT* I )=-X<N7. I ) 

PRINT EAR VALUES. 

7 MRITEI6.27I 
WRiTEI6«261 
MRITf(B«29l RAflUXI 
00 8 1=1. NRl 
J= !♦! 

8 MRITFI6.29I RAIJI.XINT.II 
GO TO 10 

9 AVS=. FALSE. 

MRITFI6.26I 

SFT-liP ANC FACTOR SUB-MATRICES IN BLOCK DIAGONAL 
COEFFICIENT MATRIX. 

10 CALL POISSl 

PRINT INPUT BLOCK. 

WRITEI6.351 
mR ITFIB. 25 I 

MRITEI6.321 IRAIIl, I=6.NR1 
00 11 J=l.NT 

11 MRITEI6.24I TDIJl* IXIJ.K-Il • K=6. NR) 

MR ITEI6.33) XI 

CARRY OUT SOLUTION OF MATRIX PROBLEM TO YIELD V-FIELD. 
CALL P0ISS2 
C 
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C ADJUST £AP VALUES* 

IFi lEAft* E0*0 ) GO TO 17 
no 16 I=4tl6 

16 X(NT« i I^XiNTU I ) 

r. 

C PRINT OUTPUT BLOCK. 

17 MRrTF46.35) 

WR ITF16.30 I 

WRITE <6, 31 I <RA( IK I =2* NR I 
DO 12 Jsl.JE 

12 WRITE (6. 34 I TO ( J ) • I X ( J . K- 1 ) t 2 • NR ) 
00 13 J=2tJE 

DO 13 K=l*15 

13 XT( J.KKXI J.KI 
C 

C DISPLAY ECUIPOTENT lAL LINE DATA. 

IFILOnp.LE.n CALL EOL INE 
no 14 J=1*JF 

14 V< ltJ)=XCTR 
DO 15 1=2.17 
K=r-i 

on 15 J=lfJF 

15 VI I « J I =XTIJ.K) 

C 

RETURN 

END 


SUBROUTINE ECLINF 

c 

C FIND ON EACH THETA MESH LINE THE RHO VALUES ASSOCIATED 

C WITH NEOL PRE-SPECIFIEO POTENTIAL VALUES. AN EOU I POTENT I AL 

C CURVE PASSES THROUGH SETS OF THESE POINTS. 

€♦ 

C«*^ 

COMMON BO.CVF.ORHTI 1C I . DRHM 1 16 I . OTHTD . UTH TR . DVDR ( 1 0 ) t 

1 OVDTI lOKOVROR^OVROZ.EDLI 10) f ETAOtIBAR.IT. JF« JTtKPCtKPOIf 

2 K pose .LOOP. KF0L.N1.N2.N3.PHD2A( 10 ) . P I . P 1 2 t P SI A ( 10 I t PS 10 ( 10 ) t R t 

3 RAOfRHMl 17).RH0PI 10) .RHUPPIlOi tRHTIlG) .RHVI 10 I » SPCHtTAU02 t 

4 THD( 10) . THMCI 2 5 ) . THSOEL . THSF AC • TH SwH I . TH SWLO. 

5 THTOI 10 ) .THTRI ID) . THl .TH2. VI 17.25 ) . VO. VR . 

6 VSMII 1D).VTI10).WI16).XCTR.XTI25 .16) .Y(16)f ZJO 

FOUl VALENCF I OTF. YI 1 ) ) . I TH . Y ( 2 ) ) . I RH . Y I 3 ) ) » I RH2tRHPt YI4) ) . 

* (RHIP.YI 51 ) .IRH2P. YI6) ) . I W( 1 ) .DELR) • I WI2) .AR) , 

« lWl3)«Z)«lW44).Z2.ZP)«IW(5I.ZlP)t(w(6) .Z2P) 

10 FORMAT(/?8X.23HTAHLE OF RHO-VALUES FOR) 

11 FORMATIBOX.ISHEUUIPOTEnTIAL LINES) 
l? FflKMATC IX.F6.2. ?X. lOF 10.3) 

13 FOBMAT(2X.5HTHETA.2X.10(F10,2) I 

f. 

C ZERO OUT RHV-ARRAY. 

on 6 1 = 1,10 

A RHV( I >=0.0 

c 

WR rTE(6.10l 
WRITF(A,11 ) 
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D r> n n 


NR = 17 

MKITE(6.13t (EOUn. I = 1,NEQL) 
f. 

DO 5 J=l.Jf 
no A l^l.MEOL 
DO 1 Kl=l,NR 
K=NR-K1 

IF(E0LU ).GT.XT<J.KJ) GO TO 2 
1 COMINUE 

r. 

? (F( Kl.GT.l ) GO TO 3 
RHV< ( i =0.0 
GO TO 5 
C 

3 F=(XT(.I.KJ-ECLI m/(XT(J,K)-XT4JtK^li) 
IF4K.GT.QI G4j TO A 
F = 4 Xf.TR-EOi 4 li ) / 4 XC Tfi-XT 4 J , Kf 1 } ) 

A RHV4 I »=RHM4K*U*F*4RHM4K*2)-RHM4K+n I 
6 WRITE46, 121 THMO( Jl, 4RHV4II. I=1.NEUL» 
f. 

RETURN 

FND 


.SOflROUT4NF R FSC AL 4 A , 2 JO. VO , THO . THF . X » V2 » IBAR #RHM, TH2 . KPOSC » 

♦ 

RIGHT-HANC SIDE CALCULATION FOR POISSON MATRIX EQUATION. 
EQUIVALENT TO CALCULATION OF SPACE CHARGE FIELD. 

* 

DIMENSION V4 111 ,C( 9) , VI 4 ID I .RH411) .ADIIL) .RH S4 1 1 . 9 ) , ARG4 16 I , 

♦ THC4 1 1.91 .THR4 11,91 ,RHSI 4 10.161 .X4 25,161 ,V24 101 ,RHM(17I 

FORMAT STATEMENTS. 

22 FORMAT41X.F5.3.F7.3I 

23 FORMAT! 16F5. 01 

PA F0RMAT4 1 IHOR FS ARRAY-1 
26 F0RMAT4 IX, 16F8. 1) 

26 FORMAT42X.16Fa.2» 

27 FORMAT410F5.0» 

2fi FORMAT! /2X,3FBH0.5X.2H-V) 

29 FORMAT! 12H0RFS VALUES.) 

30 FnRMAT4/2X,3FRHO,AX,lHC,2X,AHTUBE,AX,5HTHETA.6X,3HRHS) 

31 FaRMAT(lX,F5.3.2IA,lH-,I 1,F10.3.F9.1) 

32 FORMAT45H TH l=, F7. 3, 7X . 5HTHF =,F8.3) 

33 F0RMAT4/AX,5FAXIAL/2X,9HP0TENTIAL/1X.12HDISTRIBUTI0N) 

C 

C A=RAOUiS CF ENTRY HOLE, 

f. ZJO=JO 4CLRRENT DENSITY). 

c tho=initial value of trajectory 5. 

C X=XT425,16), SOURCE TERMS. 

c v2=energy classes. 

f. TH2=C0NF ANGLE. 

TH1=TH0 
IF = 11 
JF= 7 
JF1=JF*-1 
JF2=JF+2 
FNT=8. 
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nnnn nnn r>nnr^r>r> nn r>n nr>r^nnnnnr> 


EC sR.ti60E-l« 

ETA=l.753F+15 

Cl =2.«En*S0«T( 2.*ETAI*ENT*ENT 

PI =3.1«1592^5 

r.VF = lRO./PI 

WRtTF(6.29l 

WRITE<6.32) IHO.THF 

SECTION FCR COMPUTING SPACE CHARGE ALONG 
TUBE CENTER LINES. 

I=8IG LOOP FOR 11 RHO VALUES FROM 0.63 TO 1.00 
J=SM. LOOP FOR 7 TUBES ( 1-2 .2-3 ,. . . . 7-8 J . 

THFsRHO .63 CROSSING ANGLE OF TRAJ. NO. 5. 

READ 11 AXIAL V-VALUESIRHO .63 TO 1.0> 

AND VSMI-ARR AY(9) . 

GENERATE ARG-ARRAY. C-ARRAY, AND VI-ARRAY. 

ARG(K)=THETA MESH VALUES IN HOLE REGION. 

C ( J)-].3.5.....13 (MULTIPLIERS!. 

VI (J!-1.55.1.40,..^ .Oi^.O.C (ENERGY CLASSES IN TUBES), 
no 1 Ksl.16 

1 ARG|K)=I80.-fL0AT(K)*.25 
no 2 J=1.JF 

Vl( J)-(v2(J)4V2(J4-1! \n» 

2 C(J»s2*J-l 

GENERATE AND PRINT RH-AKRAY AND V-ARRAY. 

RH( I ) ARRAY-11 AXIAL RHO MESH VALUES 
FROM .633 TO 1.000 

V ( I)-POTENTIAL DISTRIBUTION AT THESE 
11 AXIAL MESH POINTS. 

WRITE(6.33) 

WRITE(6.2fll 
DO 3 1 = 1, IF 
RHl I )=RHM( 1^6) 

V( f)=-X( 25, I +5) 

3 WR ITE(6.22) RH( I ), V( I ) 

GENERATE THETA-ARRAY AND SPREAD THROUGHOUT 
SPACE CHARGE REGION (BASED ON TRAJECTORY NO. 5). 
DTH=TH1-THF 
nRO=l.-RHM( 7 ) 

TO=DTH/nRO 
DO 5 1=1, IF 
16=14-6 

OIFF=RHM( 16 !-RHM(7) 

TI = THF-4TO«OIFF 
T2 = Tl+( IfiO.-Tl 1/8. 

DEL=( 180.-T2)/A. 
no 4 J=1.JF2 

THOI I. J)=180.-FL0AT( J)*DEL 

4 THR(I,J)=THD(I,J)/CVF 

5 CONTINUE 

GENERATE AO-ARRAY (DISTANCE BETmEEN 
TRAJECTORIES 4 AND 5 AT EACH OF 
11 RHO MESH LINES), 
on 6 1=1 , IF 

TX1=RH( I )*TAMTHR( 1,5) ) 

TX2=RH( I )4TAN(THR( 1.4)) 
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non nn r>r^n rvr\n 


6 &D( I »*TX 2 -TX 1 


CALCULATE SOURCE TERMS ECR ALL TRAJECTORY 
CROSS INliS OF RHO MESH LINES. 

AJ=A*A*Z JO 

no <3 1=1. IF 

IFIKPnSC.NE.fll WRIT6I6.30) 
no «3 J = 1.JF 
E1 = VI( JI-VI I ) 

IF I £1 .GT.O. » GO TO 7 
RhSi i..n=n.o 
GO TO a 
7 F? = SORTIEn 

ST = SINITHR( I .J) i 

DEN=CI*S0RT1 VO**3l*E2*RHI I )«ST*AOI I) 

RHSII . J) = (C(JM‘AJ) /OEN 
a J1=.M-1 
IC=C( J» 

IFIKPOSC.NE.O) MRITEI6.31) RHI 1) . I C . J . J1 . THD 1 1 . J ) . RHSII . J ) 
9 CONTINUE 


SECTION FOR INTERPOLATING SPACE CHARGE FIELD AT MESH POINTS, 
no 15 1 = 1,10 
no 15 K= 1.16 
on 10 j=i ,7 

IFI ARGIKI.GT.THOI I. J) ) GO TO 11 

10 CONTINUE 

1 1 IFI J.NE . 1 I GO TO 13 
0 ELRHS= 0 .-RHSI I .J) 

T 1 = 180 .-THD( I. J 1 

12 T 2 = ARGIK l-TFCI I .J ) 

F=T 2 /T 1 

IFIARGIKI.LT.THOI I.JII GO TO 14 
T 1 =RHS( I ,J» +F*OELRHS 
RHSII I ,K) = IRFSI I I.KM^TU / 2 . 

GO TO 15 

13 T 1 =TH 0 I I . J -1 l-THOI I, JI 
UFLRHS=RHSI I .J-ll-RHSI I. JJ 

GO TO 12 

14 RHSII I ,K l=0.C 

15 CONTINUE 

IFIKPOSC.EO.O ) GO TO 17 

PR INT. 

MR ITEia . 24 I 

WRITEI 6 . 26 I lAKGIK). K= 1 , 16 ) 

no 16 1 = 1.10 

16 MRITEI 6 . 25 ) (RHSIII.K). K= 1 . 16 ) 

2 FRQ OUT PHS PREPARATORY TO 
LOADING SOURCE TERMS. 

17 no 18 J= 2.24 
no 18 K= 1.15 

18 XIJ.K)= 0 .n 

IFI IBAR.NE.O I GO TO 20 
no 19 1 = 1.16 

19 X 125 , I 1 = 0.0 
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f. 

20 J2=13«-TH2/15. 

no 21 4=J2.2« 
no 21 K=6,15 
Jl=25-J 
Kl=K-5 

21 X ( J ,KI =RHS I (Kl, Jl) 
C 

ftPTURN 

FNn 


SUHROIJTINE PSICAL (Rf Z .PS 1 .BC. ASM) 

C* 

C* CALCUIATFS PSI. THE MAGNETIC FLUX. OF THE MAGNETIC 

C* FIELD AT each MESH POINT IN THE REGION OF THE 

C* ENTRy HOLE, 

f ♦ 

PI=3.141592F 5 
RF=R 
fiSC=R*R 
X= Z/ASM 
X2 = X*X 

AT = ATAN2I ASM.-Z) 

r. 

T= 1.+X2 

TO = PI-AT-X/T 

TA = 3.-5.*X2 

T6= 3.-7-*X2*( 2.-X2) 

T= 21,-5.*X2 

T= 1S.-X2*T 

T8= 5.-3.*X2*T 

T= 36.-5.*X2 

T= 64.-X24T 

T= 20 .-X 2 *r 

T1C=9.4( 15.-11. *X2»T) 

T= 11. -X2 
T= 7.*X2*T 
T= 19H.-T 
T= 154.-X2*T 
T= 35.-X24T 
T12 = 21.-13.*X?*T 
C 

Xl= ASM# I 1.+X2) 

V= (RFZX1)##2 
C 

T= I21.#TlO)/S.#(ll.*Il2#V)/4. 

T= 7.#Tfl#0.5*V#T 
T=( 5.#T6)Z3. *.125#V#T 
T= T4/4.+.125#V#T 
T= l.+V#T 

T= TO-I ASM#X*RS0#T ) Z1 Xl##3) 

PSI= BO*RSa*T 

C . ■ . ' 

RETURN 

ENO 
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SUBROUTINE BRBZ (R 1 2 « BR *BZ *B0 * ASM I 
C* 

C* CALCULATES BR AND BZf THE R AND Z COMPONENTS OF THE 

C* MAGNETIC FIELD AT EACH MESH POINT IN THE REGION 

C* OF T^E ENTRY HOLE. 

C* 

PI-3. IA1592E5 
IF(Z.EO.n.O) Zc.0004 
X=ZZASM 
X2-X««2 
C 

Tl = l. 

T2 = l. 

T3=1.-5.*X2 

T4=3.-5.*X2 

T5 = 3.-7.*X2*(6.-5.«‘X2I 

T6=3.-7.*X2*(2.-X2) 

T7 = l.-X2*( 27.-7.<‘X2*(9.-3.*X2) ) 

T8 = 5.-3.*X2*(15.-X2>t‘<21.-5.*X2)) 

T9=30.-X2>M 28.-5.*X2 I 
T9*5.-11.4X2*(20.-3.4X2*T9I 
T10=54.-X2* (36.-5. ♦X2) 

T10=20.-X2*T10 
T10=9-*( 15.-11. -CTIO^XZI 
C 

T11=18.-X2*(9.-X2) 

T11=195.-143.*X2*( 10.-X2*T11I 
T11*3--X2*T11 
T12=-X2*( 10C1.-91.4X2) 

T12=-X24(2574.*T12) 

T12=-X2*(20C2.*T12) 

T12=-X2*( 455.+T12) 

T12all.*(21.*T12) 

TlBa(ASM*Z )/IASM*ASM*Z*Z) 

C 

ZK»1.4X2 

Al=ASM«ZK 

A2=ASM*ZK*X 

C 

Cl=(-2.*80)/(PI4ASM*ZK4ZK) 

C2=(-4.*C1*X )/Al 
C3=C2/A2 

C4=(-6.4C3*X)/A1 

C5=C4/A2 

C6=(-40.4C54X) /A1 
C7=(3.*C6)/A2 
C8=(-14.4C74X) /A1 
C9=C8/A2 

C10=(-4.*C94X) /A1 
Cll*(45.4C10 )/A2 

C12s(-2.4C114X)/Al 

C 

AT=ATAN2( ASM,-Z) 

BB=(80/PI)4(PI-AT-T1B) 

8 1=C 1*T 1 
8 2=C 2*T 2 
B 3=C 34T 3 
B 4sC 44T 4 
B 5=C 5*T 5 
8 6*C 64T 6 
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B 7=C 7^T 7 
B B^C 8«T 8 
8 9=C 9*T 9 
B10*C10«T10 
811=C11«T11 
B12=C12»T12 
C 

R2-R/2. 

RSsR2«R2 

Tl=B10-RS^B12/36* 

Tl=B8-RS^Tl/25« 

Tl=B6-RS^Tl/16* 

Tl=84-RS^Tl/9« 

Tl = B2-RS<»Tl/4. 
B2*B8-RS*T1 
C 

Tl=B9-Bll^RS/30. 

Tl=B7-RS^Tl/20* 

Tl=B5-RS^Tl/12* 

T1=83-RS^T1/^* 

Tl=81-RS*Tl/2* 

BR=-R2<*Tl 

C 

RETURN 

END 
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APPENDIX C 


ALPH(IO) 

ALPHA 

AR 

ASM 

AVS 

BO 

BR 

BZ 

CONFLO 

CVF 

DELI(IO) 

DELID 

DELR 

DRHM(16) 

DRHT(IO) 

DTH 

DTHTD 

DTHTR 

DVDR(IO) 

DVDT(IO) 

DVMDRH 

DVMDTH 

DVRDR 

DVRDZ 

EQL(IO) 
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FORTRAN SYMBOLS 

array of entry angles, a, rad 
entry angle, a, deg 
rectilinear coordinate, r 
radius of entry hole, a, cm 

axial boundary value logical switch TRUE indicates presence of an 
axial spike 

magnetic flux density, Bg 
r- component of magnetic field, B^. 
z- component of magnetic field, B_ 

switch for type of flow (confined flow when not equal to 0; Brillouin flow 
when equal to 0) 

conversion factor, deg to rad 

array, rad 

entry angles of electron classes, 6., deg 
step size Ar for r-z trajectory integration 
Ap of mesh 
Ap of trajectories 

A0 of trajectories (equivalent to DTHTR below), rad 

A0 of trajectories (step in 6 for p-9 trajectory integration), deg 

A0 of trajectories, rad 

9V/9p array 

9V/90 array 

9Vm/9p array 

aV^/90 

9Vj./9r 

aVj./9z 

array of values for which equipotential lines are to be calculated 



ETAO 

EO 

FLOOP 

IBAR 

IT 

Jl, J2, JF, JT 
KPO, KPOI 

KPOSC 

K9 

LKTR 

LOOP 

NEQL 

Nl, N2, N3 

NF 

NR 

NT 

PHD 

PHD2 

PHD2A(10) 
PI, PI2 
PSI 

PSIA(IO) 

PSIO 

R 

RA(17) 

RAD 
RH, RHP 
RHIP 
RH2,RH2P 


initial electron charge- to- mass ratio 
permittivity of free space, €q 

storage for number of iteration loops in space charge calculation 
indicator for axial boundary values 
iteration coimter 
looping indexes 

main output frequency control of trajectory data (e.g. , 5 means to print 
every fifth point); KPOI is constant; KPO may be varied 

output control for space charge values 

initialization switch 

coimter for number of loops 

indicator switch for space charge looping 

number of equipotential lines to be calculated 

control indexes for trajectory calculations 

number of 0-mesh lines, usually 26 

number of p-mesh points, usually 17 

number of trajectories, usually 9 

cp{d(p/dt), where cp is azimuthal angle 

. 2 

(p 

2 

cp array for output 

77 , 77 ^ 

magnetic flux, ;// 
j// array 

initial magnetic flux, t//q 
radius of collector, R, cm 
array of r-mesh values 
radius of entry hole, a, cm 
p, p'(3p/30) 

^ 2 ’ ^2 
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RHM(17) 

p-mesh values 

RHOP(IO) 

P' 

RHOPP(IO) 

p" 

RHRZ 

p in v-z coordinate system 

RHT(IO) 

p- values for each trajectory 

RHV(IO) 

p- values for equipotential lines 

SPCH 

space charge indicator switch 

ST, ST2 

2 

sin 9, sin 6 

TO, T1,T2, . . . 

temporary storages 

TA(25) 

0-mesh values 

TAU02 

’o 

TH 

variable of integration, 0 

THO 

initial cone angle, deg 

THl 

initial cone angle, deg 

TH2 

collector angle at which surface charge is 

THD(IO) 

0 array 

THF 

^f 

THMD(25) 

0-mesh values, deg 

THSDEL 

0 step- size increment, deg 

THSW 

0 switch for step- size control 

THTD(IO) 

0- coordinate of trajectory, deg 

THD2 

02 

THTR 

0- coordinate of trajectory, rad 

THRZ 

0 in r-z coordinate system 

V(17,25) 

potential inside spherical collector 

VO 

initial normalized electrical potential, Vq 

VR 

interpolated potential in r- z system 

VSMI(IO) 

v^ array; beam velocities 

VT(IO) 

array of initial voltages 

W(16) 

Runge-Kutta storage area in r-z system 
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XCTR 
XT(25, 16) 
Y(16) 
Z,ZP 
ZIP 

Z2, Z2P 


value of source term at center of sphere 
source term needed for solution to Poisson's equation 
Runge-Kutta storage area in p-6 system 
z, z'(9z/0r) . . 

: 

> miscellaneous storages for independent variables in 

„ cylindrical coordinate system 

^ 2 ’ ^2 

y 


2 

ZJO initial current density, Jq, A/cm 

ZIO initial current, Iq, A 

ZKP perveance, Kp, 

ZED rectilinear coordinate, z 
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APPENDIX D 


FLOW CHARTS 


Flow charts are presented only for those five subroutines, INTGRN, DE, RZTRAJ, 
VFIELD, RHSCAL, which might not otherwise be easily understood. Diamond- shaped 
blocks indicate decision-making points in the subroutines. Numbers appearing on blocks 
refer to the corresponding FORTRAN statement numbers in the subroutines. 


Subroutine INTGRN 
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Subroutine RHSCAL 
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APPENDIX E 


MAGNETIC FIELD CALCULATIONS 


B°(0,z) 

IT 


71 - tan' 


■ft)- 


1 + x' 


(El) 


where x = z/a. Then, 


z) = 


dHI,, 


alV, 



1 

ffa 



1 


M' 


1 - 5x2 

i?a® j 


48BqX 

,,3-5x2 




B^(0, z) = - ^ X 3 - 
^ „5 / ,\6 




B^(0, z) =- 


4 X 


480BqX 




(3 - 14x^ + 7x^) 


T /TT 12 X 480Brt o 4 A 

B^(0, z) 1 (1 - 27x2 ^ gg^4 _ 2ix®) 




(E2) 


sr(o,z)=-- 


12 X 14 X 


480BgX 




(5 - 45x^ + 63x^ - 15x®) 


B^{0, z) = - 


12 X 14X 480Bn 


9(l.x2f 


! (5 - 220x2 ^ ggQj^4 _ q24x^ + 165x®) 


Y ^ ^ 480Bj-^X n 4 A R 

BC(0,z) iL (135 - 1980x^ + 5346x^ - 3564x^ + 495x**) 

10/, 2\“ 

ffa (1 + X 1 


y, 14 X 45 X 48 X 480B. o 4 a 

B^-‘^(0j z) - IQA-U-^ L 14*>Av^ . OC'T/I^D 




' (3 - 195X-" + 1430x‘* - 2574x° + 1287x^ - 143x^°) 


ytt 28 X 45 X 48 X 480B«x 

bP(0, z) 0 




(231 - 5005x^ + 22022x‘‘ - 28314x® + llOllx® - 1001x^°) I 
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To simplify these expressions, let 




CO 

II 

1 

5x2 






T4 = 3 - 

5x2 





CO 

11 

- 42x2 + 35x'^ 

= 3 - 

7x^ 

'(6- 

5x2) 


Te-3 

- 14x2 ^ rj^4^ 

-3 - 

7x2 

(2- 

x2) 

T7 = 1 

- 27x2 -t 

■ 63x'^ - 2 lx® 

- 1 - 

3x2 

[- 

7x2(3 

Tg = 5- 

■ 45x2 + 

63x^ - 15x® = 

5 - 3x2 1 

15 - 

x2(21 

+ OOOx*^ 

- 924x® 

+ 165x® = 5 - 

11x2 

[20 

- 3x2[30 - 


’(28 - 5x2)]} 


Tjq = 135 - 1980x2 + 5346x^ - 3564x® + 495x® 


= 9(15 - 11x2 120 _ ^ 2 ^^ _ ^2(3g _ 5^2^^ 


Tjj = 3 - 195x2 ^ 1430x4 _ 2574x® + 1287x^ - 143x^° 


= 3 - 13x2^15 - 11x2 jlO _ ^2 Jig _ ^2(g _ ^2jj| 


Ti 2 = 231 - 5005x2 + 22022x'^ - 283 14x® + llOllx® - 1001x^° 

- ll|21 - 13x2^35 - x2|i 54 - x2[l98 - 7x2(11 - x2) 
Then, substituting these expressions for those in (E2), we now have 


^ (E3) 
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T 

JO, z) L 


b“(0,z) =- 


.a(l+x2)" 

8BqX 


Tra 




B™(0,z) 


B^J(0,z) =- 


8B0T3 


.3(1 . xf 

48BQXT4 


7ra 


‘(‘-f 


B>,z)= lA 


Tra 




B>,z) ,^^^30 BqxT6 


Tra 




™ 12 X 480BoT„ 

By“(0,z) = °-J 

^ \8 


Tra 


bJ“^(0,z) =. 


12 X 14 X 480BgxTg 


ira 


'(‘”f 


^ 12 X 14 X 480BnTQ 

Bf (0, z) °-9 

.a9(l.xf° 


yt 14X 45 X 48X 480BnTii 

Bf(0,z)= ° 


11 


\12 


B^(0, z) = - 


28 X 45 X 48 X 480BqxT^2 
.al2(l + x2)^^ 


(E4) 
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Since 


B. 


.(r, z) = - B^(0, z) i: + -±- b“^(0, z)k^ - -i_ B^(0, z)Uf + -±- B^(0, z)(lV - . . 

2 1:2: \2/ 2:3: \2/ 3: 4: \2/ 


(E5) 


We make the substitution t2 = r/2, and nest as follows, 


z) - -rj 


and similarly, 


bI - -1^ B™ 


2 Ib^ 


lx2\ 2x3 3 x 4 


B 


vn 


2 /gEX ^ ^2 


4X 5 V 5 x 6 


(E6) 


B2,(r, z) = B^ - r| 


B°-:2(3lV_^ (3VI 


4^ 


B 


vin 


‘2fex_:2gxn 


(E7) 
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APPENDIX F 


DERIVATION OF EQUATIONS OF MOTION 


The equations of motion as derived in reference 1, equations (BIT) and (B18), are 


dA 

p'e = p6^ + p sin^e (p^ - p"^^ + — - .2- sin e - -2. p sin 0 <p — ^ (FI) 

j^2 0p R ^ R dp 

B = - + sin 0 cos 6 + — 2 — — - -2_ cos 0 . - -2- sin 0 <p — ^ (F2) 

p p2j^2 30 pR ^ pR 30 


Solving these two equations simultaneously for p" 3 delds 


p"=p + £ji2!ll!+_>i_f, +ief 

0^ R^^ ^ 


where 


f =liy - sin e K (p 
^ R dp ^ 


fn = — - pR cos 0 A,-^ 
^30 ^ 


sin g cos 6 ^ _ p't] ^ 
0^ p^R^0^ 


(F3) 


p Sin 0 ^ 


3A 




(p 


- pR sin 0 (p 


dp 

3A 




30 


(F4a) 


or, using the relations described in appendix D of reference 1, these expressions may 
be rewritten 


^ R 3p 30 

f =1Y + p2^ 0 ^ 

^ 30 3p 


(F4b) 


Away from the entry hole, the magnetic field terms <p , A^, and are negli- 
gible. Therefore, in this region equation (F3) may be expressed ' 
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(F5) 


T] ^ 3V , 2p'^ _ p’T7 ^ 8V 
r2^2 dp p p2r2^2 30 


using equation (B16) of reference 1 as 


02 


27]Vo(e. + V) 


R2p2 



(F6) 


and, with the appropriate substitutions, we may express this equation in cylindrical 
coordinates 




0z Sr 

— 

' e. + V 

2 

1 

_1 + (z')2_ 
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APPENDIX G 


MATHEMATICAL SYMBOLS 

2 

A area of annulus, cm 

a radius of collector entry hole, cm 

Bj. r- component of potential field, 9Vj^/8r, G/cm 
z- component of potential field, G/cm 

Bq magnetic flux density, G 

b radius of electron beam at injection, cm 

E V. + V 

I current, A 

2 

J current density. A/cm 

Kp perveance, A/V^^^ 

R unreduced radius of spherical collector, cm 

r, z cylindrical coordinates 

u electron energy 

V normalized electric potential 

Vjjj magnetic field potential 

av/3r 
V, av/az 

Vp av/ap 

Vg av/ae 

XU 

Vj i normalized injection energy class 

Cq permittivity of free space, 8. 86x10”^^ F/cm 

? 7 g electron charge-mass ratio, chargeAg><10^^ 

p normalized radius vector 

p unnormalized radius vector, cm 

Q 

Pg space charge density. C/cm 

Tq R/(2r]Vo)^/^ 

0, (p polar and azimuthal spherical coordinate angles 
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